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1.  INTRODUCTION  AND  RESEARCH  OBJECTIVES 


This  is  a  final  technical  report  on  work  carried  out  at  the  Gas  Turbine  Laboratory  at  MIT,  in  a 
multi-investigator  effort  on  unsteady  flow  phenomena  in  turbomachincs.  Suppcwi  for  this 
program  was  provided  by  the  Air  Force  Office  of  Scientific  Research  under  Grant  Number 
AFOSR-90-0035,  Major  D.  Fant,  Program  Manager. 

The  present  report  gives  a  summary  of  the  worit  for  the  period  10/89  -  10/92,  but  is  not 
intended  to  be  the  primary  source  and  the  referenced  repons  and  publications  can  be  consulted  for 
more  detailed  information  and  background.  These  describe  the  research  in  considerable  depth. 

Within  the  general  area  of  turbomachinery  fluid  dynamics,  four  separate  tasks  wwe  specified 
under  this  contract.  These  are,  in  brief: 

l.  The  influence  of  inlet  temperature  non-uniformities  on  turbine  heat  transfer  and 
measurement. 

n.  Development  and  application  of  computational  techniques  for  unsteady  flows  and 
assessment  of  unsteady  losses  in  rotOT/stator  interactions. 

m.  Unsteady  phenomena,  inlet  distortion  and  flow  instabilities  in  multistage  compressors, 
including  experimental  and  analytical  investigations  of  the  structure  of  instability 
precursors,  nature  of  rotating  stall  modes,  and  detailed  blade  passage  aerodynamics  and 
stall  onset  in  advanced  blading  geometries. 

IV.  Unsteady  blade-vortex  street  interaction  in  transonic  cascades  including  models  to  describe 
the  vOTtex  wake-compressor  blade  interactions. 

In  atklition  to  these  tasks,  the  multi-investigatoT  contract  encompassed  the  Air  Foee 
Research  in  Aero  Propulsion  Technology  (AFRAPT)  Program.  The  work  carried  out  in  each  of 
the  tasks  will  be  described  in  the  next  section.  Publications  generated  arc  given  in  the  individual 
task  descriptions,  and  an  overall  list  of  presentatirais  and  publications  associated  with  die  projects 
appears  in  Section  4. 
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TASK  1:  THE  INFLUENCE  OF  INLET  TEMPERATURE  NONUNIFORMITIES 

ON  TURBINE  HEAT  TRANSFER  AND  MEASUREMENT 

(Investigators:  A.H.  Epstein,  G.R.  Gucnette,  T.  Shang,  M.  Graf) 

1.0  Executive  Summary 

This  is  a  technical  progress  rcpon  on  the  research  conducted  on  the  turbine  inlet 
temperature  distortion  cooperative  research  program  through  December  1992.  The  t^jcctivcs  of 
the  program  are  to  elucidate  the  role  of  inlet  temperature  nonuniformity  and  turbulence  on  turbine 
rotor  heat  transfer  through  detailed  experiments  and  calculations.  The  sponsors  are  SNECMA, 
Rolls-Royce,  Inc.,  US  Air  Force  Office  of  Scientific  Research,  and  NASA  Lewis  Research  Center, 
The  effort  is  focused  about  a  set  of  experiments  on  the  MIT  Blowdown  Turbine  facility  in  which 
the  effects  of  inlet  temperature  distortion  on  a  transonic  turbine  stage  will  be  directly  measured. 
This  data  will  then  be  compared  with  CFD  and  analydcal  results. 

The  activities  to  date  have  concentrated  on  the  refurbishment  of  the  test  rig,  the  addition  of 
turbulence  grids,  the  upgrading  of  heat  transfer  and  aerodynamic  instrumentation,  the  adaptation  of 
a  3-D,  inviscid,  unsteady  code  to  the  measurement  geometry,  final  assembly  of  the  aerodynamic 
insmimentadon  and  its  calibration,  initial  NGV-only  testing,  and  the  resolution  of  problems 
encountered  during  the  initial  tests.  The  test  rig  refurbishment  includes  cleaning  of  the  temperature 
distortion  generator  and  reconfiguration  of  the  generate*-  to  simultaneously  produce  different  radial 
and  circumferential  distortion  in  separate  120°  sectors. 

Instrumentation  upgrading  consists  of  reinstrumentation  of  the  rotor  blades  with  heat  flux 
transducers,  the  addition  of  circumferential  rake  translators  upstream  and  downstream  the 
turbine,  and  the  design  and  construction  of  pressure  and  temperature  rakes  for  the  translators.  The 
temperature  instrumentation  measures  time-averaged  total  temperature,  while  both  time-avmged 
and  time-resolved  rakes  were  fabricated  for  total  pressure.  Wall  static  pressure  taps  have  been 
added  at  the  inner  and  outer  annulus  wall  downstream  of  both  the  NGV's  and  the  rotor.  New 
signal  conditioning  and  data  acquisition  equipment  to  acctnnmodate  the  80  additional  aerodynamic 
measurement  transducers  has  also  been  procured  or  constructed. 


The  temperature  and  pressure  rakes  were  assembled  and  calibration.  Considerable  noise 
and  induced  signals  were  observed  on  the  exit  temperature  rake  during  traverse  shakedow  n  in 
vacuum.  This  prdslem  was  solved  with  an  improved  grounding  scheme.  The  overall  DC 
accuracy  and  stability  appear  to  be  within  0.2%,  sufficient  for  the  planned  test  program. 

The  instrumentation  performed  as  designed  during  the  first  NGV-only  blowdown  tests. 
The  inlet  rake  mapped  the  distortion  generator  exit  flowficld  and  the  turbine  exit  rake  could  readily 
resolved  the  NGV  wake  structure.  During  the  first  test,  cabling  from  the  turbine  exit  translator 
became  entangled  in  it  drive  chain.  The  exit  cabling  arrangement  was  modified  to  resolve  this 
problem. 

Turbulence  screens  have  been  designed  and  fabricated  which  match  the  turbulence 
spectrum  suggested  by  Rolls-Royce  as  being  representative  of  engine  scale.  These  screens  have 
been  tested  in  a  wind  tunnel  with  satisfactory  results. 

A  3-D,  unsteady,  Euler  code  written  by  Saxer  and  Giles  is  being  adapted  to  the  test 
geometry.  The  code  will  be  run  with  three  NGV’s  and  five  rotor  passages  to  simulate  the  36 
NGV,  61  rotor  blade  count  of  the  experiment. 

Two  problems  were  encountered  with  the  temperature  and  distortion  generator  - 
inadequate  circumferential  distortion  and  particulate  cemtamination.  The  low  level  of 
circumferential  distortion  achieved  was  traced  to  the  storage  matrix  having  15  times  nxire  thermal 
conductivity  in  the  circumferential  direction  compared  to  the  radial  one,  due  to  poor  fabrication. 
This  has  the  effect  of  quickly  washing  out  circumferential  distortions.  The  approach  adopted  was 
to  double  the  number  of  spot  heaters  and  run  them  at  very  high  current  levels  (just  below  burnout 
as  established  experimentally)  immediately  preceding  the  blowdown  tests.  The  particle 
contamination  of  the  distortion  generator  has  been  a  continuing  problem  and  efforts  to  clean  the 
unit  have  not  been  successful.  The  approach  adopted  here  was  to  fabricate  and  insert  a  particle 
filter  at  the  distortion  generator  exit. 

An  operational  problem  has  developed  due  to  the  ten-fold  increase  in  the  price  of  Freon- 12, 
to  about  $1200  per  test  -  well  beyond  the  budget  We  plan  on  using  Freon  for  the  first  tests  only 


5 


for  continuity  and  then  converting  to  an  argon-C02  gas  mixture. 

The  current  schedule  is  to  do  the  NGV-only  testing  in  January  and  start  the  fuU-siagc 
testing  in  February  1993. 

2.0  Introduction 

The  objective  of  this  work  is  to  understand  how  turbine  inlet  distortion  and  tuiixilence 
influence  turbine  heat  transfer  and  aerodynamics.  It  is  centered  on  a  set  of  detailed  spatially  and 
temporally  resolved  measurements  of  a  high  pressure  ratio  turbine  stage  in  a  short  duration  turbine 
test  facility.  Numerical  simulauon  and  analytical  modelling  will  then  be  used  to  interpret  these 
measurements  and  help  draw  conclusions  of  use  to  a  turbine  designer. 

The  experimental  work  will  be  carried  out  in  the  MIT  Blowdown  Turbine  Facility.  This  rig 
can  simulate  all  of  the  non-dimensional  fluid  mechanic  parameters  known  to  be  important  in 
turbine  aerodynamics  (Reynolds  No.,  Rossby  no.,  Mach  No.,  gas-metal  temperature  ratios,  Prandtl 
No.,  ratio  of  specific  heats,  etc.)  for  approximately  300  ms.  Turbine  corrected  speed  and  weight 
flow  are  held  constant  to  better  than  0.5%  over  this  period.  During  the  test,  approximately  20 
million  data  points  can  be  taken  from  a  variety  of  heat  transfer  and  aerodynamic  instrumentation. 

The  program  plan  is  to  setup  a  temperature  distortion  generator  to  produce  both 
circumferential  (hot  spots)  and  radial  temperature  distortions  simultaneously  in  separate  120*’ 
sectors.  The  time  resolution  of  the  instrumentation  will  permit  differentiation  amemg  the  sectors. 
The  test  plan  fcR*  this  year  is  to  first  qrerate  the  tunnel  without  the  rotor  so  as  to  permit 
aerodynamic  mapping  of  the  nozzle  guide  vane  (NGV)  exit  flov^eld.  (Since  the  nozzles  are 
choked,  the  flowfield  may  approximate  that  widi  the  rotor  in  place.)  The  romr  will  then  be  added 
and  the  heat  transfer  aixl  stage  aerodynamic  data  taken. 

To  date,  the  effent  has  ccmcentrated  on  readying  the  facility  and  its  instruiirentaticm  for  ttese 
tests,  and  adaptation  of  numerical  tools.  For  the  sake  of  this  report,  we  will  divule  this  effort  into 
categories  -  (1)  facility  refurbishment  and  addition  of  inlet  and  outlet  rake  translators;  (2)  design 
and  testing  of  turbulence  grids;  (3)  new  aerodynamic  instrumentation  and  signal  conditioning 
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equipment;  (4)  new  rotor  blade  heat  flux  instrumentation;  (5)  CFD  code  adaptation  and  checkout; 
(6)  instrumentation  preparation  and  calibration;  (7)  NGV-only  testing;  and  (8)  problems  uncovered 
and  their  solutions.  In  the  following  sections,  we  examine  each  item  in  turn. 

3.0  Facility  Refurbishment  and  Modification 

Three  major  items  arc  required  to  ready  the  facility  for  the  distortion  testing  - 
refurbishment  of  the  facility,  cleaning  and  retrofitting  of  the  distortion  generator,  aitd  the  addition 
of  new  upstream  and  downstream  circumferential  translation  capability.  TTie  tunnel  refurbishment 
consists  of  disassembly;  cleaning;  bearing,  seal,  and  lubricant  replacement;  etc.,  as  well  as 
checkout  of  the  control,  heating,  eddy  brake,  and  safety  systems. 

3.1  Distortion  Generator 

The  distortion  generator  was  cleaned  and  reconfigured.  The  unit  (0.5  m  dia.  x  1  m  long) 
was  cleaned  in  an  ultrasonic  cleaner  to  remove  the  particles  which  caused  extensive 
instrumentation  damage  during  the  previous  test  series.  These  appeared  to  be  25  pm  diameter, 
clear  particles,  probably  silica  from  grit  blasting  at  the  honeycomb  manufaemrer.  The  unit  was 
cleaned  repeatedly,  the  cleaning  fluid  filtered,  and  the  particle  number  noted.  The  number  of 
particles  was  reduced  on  each  cleaning  until  none  remained. 

The  heaters  from  one  of  three  radial  distortion  sectOTS  of  the  generator  was  removed  and 
the  unit  rethreaded  as  a  circumferential  distortion  generator  (Fig.  1).  In  this  segment,  the  heatm 
are  6  mm  in  circumferential  extent  and  spaced  2  nozzle  guide  vane  pitches  apart  to  approximate 
combustor  nozzle  spacing.  Thirty  thermocouples  measuring  metal  traiperature  were  dien  spotted 
about  the  matrix  to  aid  in  pre>test  setup.  The  heaters  and  thermocouples  performed  satisfactorily  in 
air.  Thegoierattn'isreadytogo. 

Different  temperature  patterns  will  be  generated  by  varying  the  tinoe  tenqierature  history  of 
the  wire  heaters.  High  currents  for  short  times  will  generate  “spiky”  distortions,  which  wUl  soak 
out  over  time,  lessening  the  temperature  modulation.  This  is  illustrated  ccmceptually  in  Fig.  2.  As 
an  aid  in  the  setup  of  the  patterns,  a  numerical  conduction  model  of  the  generator  is  being  written 
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and  will  be  fit  to  the  experimental  data. 

3.2  Instrument  Translators 

The  explicit  inclusion  of  circiunferential  distortions  and  the  desire  to  take  high  accuracy 
aerodynamic  data  necessitated  the  addition  of  circumferential  traversing  capability  to  the  MIT 
blowdown  tunnel.  Traversing  is  required  at  the  stage  inlet  to  document  the  turbine  inflow  for  both 
heat  transfer  and  aerodynamics.  Stage  performance  measurements  demand  that  the  rotcn'  outflow 
be  measured  as  well.  The  inlet  translator  was  a  part  of  the  baseline  program,  whereas  the 
participation  of  SNECMA  has  permitted  the  addition  of  a  unit  at  the  rotw  oudet 

The  design  goals  were  similar  for  both  units:  240°  of  annular  transverse  in  300  ms  (to 
survey  both  0TDF  and  RTDF  sectors  in  the  useful  test  time),  constant  speed  (to  aid  in  data 
reduction  and  interpretation),  one  revolution  maximum  travel  (to  permit  any  two  of  the  three  120° 
sectors  to  be  measured),  ready  installation  and  removal  of  the  rake  instrumentation  (for  calibration 
checks  and  sensor  replacement,  if  necessary),  and,  of  course,  minimum  technical  risk  and  cost. 

We  wished  to  position  the  upstream  translator  as  close  to  the  NGV  leading  edges  as  possible  and 
spent  considerable  time  wuiking  on  designs  which  would  be  accommodated  within  the  current  test 
section.  However,  there  was  very  little  space  to  accommodate  the  unit,  majOT  modifications  would 
be  required,  bearings  were  not  available  in  the  sizes  needed,  and  the  high  inertia  made  the  drive 
very  large.  Thus,  fitting  the  unit  within  the  test  s^tion  at  the  current  rake  location  was  quite 
unattractive  because  of  the  engineering  difficulty,  cost,  aral  risk. 

The  approach  adopted  instead  was  to  fabricate  an  addititmal  casing  (pressure  vessel)  for  the 
translator,  which  was  inserted  between  the  disttntion  generamr  and  the  test  sectitm  (Fig.  3).  A 
cross-section  of  the  unit  is  shown  in  Fig.  4.  In  this  implementation,  the  rakes  are  noounted  to  the 
uiner  annulus  wall  which  is  rotated  by  a  gear  arrangement.  This  lowers  the  angular  moment  of 
inertia,  compared  to  rotating  the  outer  annulus  wall,  reducing  the  drive  requirements  to  about  5 
KW.  A  computer  controlled  DC  servo  drive  runs  the  translator  (Fig.  5).  The  unit  is  water-cooled 
on  the  upstream  side  to  restrict  pre-test  heating  of  sensors  and  bearings.  The  signals  from  the 
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transducers  arc  brought  out  of  the  rotary  system  by  coiled  wires  (thus  the  one  revolution 
maximum  travel)  and  out  of  the  pressure  vessel  by  a  vacuum  feedthrough.  A  measured  time  trace 
of  the  translator  motion  is  shown  in  Fig.  6.  The  unit  operates  satisfactorily. 

The  rotor  exit  translator  was  a  much  more  difficult  task  in  that  there  was  no  choice  Ixit  to  fit 
it  within  the  existing  test  section  and  the  geometric  constraints  are  severe.  In  this  case,  the  dcswirl 
vanes  downstream  of  the  rotor  are  replaced  by  the  translator.  The  tunnel  flow  path  is  shown  in 
Fig.  7,  while  the  detail  of  the  translator  is  given  in  Fig.  8.  As  for  the  upstream  translator,  the 
probes  are  mounted  on  the  inner  annulus  flow  path  wall.  This  wall  is  mounted  on  thin  line 
bearings  and  chain-driven  (the  chain  is  in  the  flow  path).  The  ribbon  cabling  is  used  which  is 
formed  into  a  spring  coil  under  the  turbine  exit  throttle  sleeve.  A  DC  servo  drive  similar  to  the 
upstream  unit  is  used  with  the  addition  of  a  gearbox  to  increase  the  torque.  As  of  this  date,  the  unit 
has  performed  satisfactorily  mechanically.  The  instnunentadon  wiring  has  not  been  completed. 

4.0  Turbulence  Generation 

The  influence  of  turbulence  on  rotor  heat  transfer  will  be  explored  by  qjcrating  the  facility 
with  high  and  low  levels  of  turbulence.  The  blowdown  turbine  naturally  produces  a  very  quiet  in¬ 
flow,  with  turbulence  levels  measured  before  the  contraction  into  the  NGV’s  of  less  than  1%. 
Engine  levels  are  considerably  more,  and  we  must  ask  the  question  “What  are  appropriate 
intensities  and  scales?”  Rolls-Royce  has  suggested  that  the  work  of  Moss  and  Oldfield  [1]  is  a 
suitable  representation  of  engine  conditions.  We  have  been  operating  under  this  assumption  in  the 
work  described  below.  More  recently,  at  the  April  program  review  and  in  subsequent 
correspondence.  Dr.  Richard  Rivir  of  the  USAF  Wright  Laboratories  has  pointed  out  that  his  work 
shows  considerably  more  energy  at  high  spatial  flequencies  than  does  that  of  Moss  and  Oldfield. 
We  are  in  the  process  of  trying  to  understand  this  discrepancy  and  would  appreciare  the  t^inions 
of  our  partners. 

Basing  the  design  criteria  on  Moss  and  Oldfield  data,  we  considered  the  suitability  of  grids 
to  produce  the  desired  spectrum.  Meshes  placed  at  two  axial  stations  were  considered:  (1)  at  the 
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NGV  leading  edge,  and  (2)  funher  upstream  before  the  contraction  at  the  boundary  layer  bleed  lips 
(Fie.  9).  A  summary  of  the  grid  propenies  is  given  in  Table  1  and  the  resultant  spectrum  given  in 
Fig.  10.  While  the  intensity  generated  by  grids  at  the  NGV  leading  edge  is  higher,  it  is 
mechanically  difficult  to  fit  the  screens  there  in  such  a  manner  as  to  be  readily  removable.  Thus, 
boundary  layer  bleed  leading  edge  placement  was  chosen.  The  grid  geometry  selected  is  a  coarse, 
square  cross-section  grid,  stamped  from  plate  (Fig.  11). 

The  grid  was  tested  in  a  low  speed  wind  tunnel  and  produced  a  spectrum  approximately 
that  of  Moss  and  Oldfield.  The  details  of  this  testing  are  given  in  the  following  sections. 


TABLE  1 

SUMMARY  OF  TURBULENCE  GRID  MESH  PROPERTIES 


Axial 

Station 

Sector 

Mesh  Size 
M  (in) 

Wire  Dia. 
b(in) 

Solidity 

s(%) 

Tu 

(%) 

WBSSM 

iKiRiyi 

1 

0.25 

.063 

44.0 

17.8 

1.14 

At  NGV 

Inlet 

2 

0.10 

.025 

43.8 

9.18 

1.13 

3 

0.05 

.012 

42.2 

5.43 

1.04 

At  B.L. 
Bleed 

RR  E)e$ign 

0.75 

.192 

44.7 

8.10 

1.18 

4.1  Experimental  Facilities 

For  the  testing  of  the  turbulence  grids,  the  MIT  12"xl2"  low  speed  wind  tunnel  was  used 
(Fig.  12).  The  turbulence  of  this  tunnel  was  found  to  have  a  negligible  effect  on  the  turlHilence 
measurements.  To  take  the  velocity  measurements  behind  the  grid,  a  Thermo-Systems  1050 
constant  temperature  anemometer  was  used,  with  a  Thermo-Systems  1213-T1.5  hot  wire.  This 
hot  wire  has  a  diameter  of  4  ^m  and  an  overheat  ratio  of  1.8.  The  anemometer  voltages  were 
recorded  using  a  Nicolet  NIC-310  recording  oscilloscope  with  16  bit  resolution.  This  provides  a 
velocity  resolution  of  0.1  m/s.  Taking  into  account  the  errors  caused  by  the  hot  wire,  pitot  tube, 
and  calibration  curve,  the  accuracy  of  the  measurements  was  calculated  to  be  approximately  ±0.5 
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m/s. 

The  calibration  of  the  tunnel  was  done  using  a  pitot  tube  attached  to  a  pressure  transducer. 
The  pitot  tube  also  acted  as  a  check  on  the  accuracy  of  the  hot  wire  during  the  tests.  The  calibration 
curve  of  anemometer  output  voltage  versus  velocity  was  then  used  to  convert  the  test  voltages  to 
velocities.  Although  the  tunnel  was  calibrated  at  the  beginning  of  each  test  sequence,  the  drift  of 
the  calibration  curve  was  checked  to  see  if  it  would  pnxiuce  an  appreciable  error  as  the  tests 
continued.  Figure  13  shows  that  the  curve  remains  approxima’ely  constant  with  time. 

Two  groups  of  tests  were  carried  out  in  the  tunnel.  The  first  group  measured  the 
turbulence  characteristics  at  a  position  half  way  between  the  turbulence  grid  and  nozzle  guide 
vanes,  and  at  the  position  of  the  nozzle  guide  vanes.  At  these  two  streamwise  locations,  six 
measurements  were  taken  at  different  transverse  positions  to  the  flow  to  check  for  changes  in 
turbulent  characteristics  in  the  transverse  direction.  The  object  of  this  first  group  of  tests  was  to 
compare  the  results  with  similar  tests  run  in  the  past  [2],  to  confirm  the  grid  was  producing  the 
correct  turbulent  characteristics. 

For  the  second  half  of  the  tests,  a  plexiglass  nozzle  section  was  added  to  the  tunnel.  This 
nozzle  is  the  same  dimensions  as  the  nozzle  section  in  the  blowdown  turbine  facility.  The 
locations  of  the  test  positions  was  the  same  as  for  the  first  half  of  the  tests.  With  these  results,  the 
fteestream  turbulence  characteristics  entering  the  nozzle  guide  vanes  is  known  for  the  blowdown 
turbine  facility.  These  turbulence  characteristics  can  then  be  compared  with  measurements  from  a 
typical  gas  turbine  engine  burner  exit. 

4.2  Data  Analysis 

When  looking  at  the  effect  of  turbulence  on  rotor  heat  transfer,  it  is  found  that  two  of  the 
rixrst  important  parameters  which  characterize  turbulence  are  turbulence  intensity,  and  the  power 
spectral  density.  The  turbulence  intensity  is  defined  as 
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The  power  spectral  density,  relating  turbulent  energy  with  wave  nunnber.  can  be  found  by 
squaring  the  finite  fouricr  transform  of  the  data  (sec  Table  2),  From  the  data,  the  one  sided  power 
spectral  density  can  be  found  as  follows 


G(k)  » 


1 


jiUNAt 


lX(f)l- 


where 


N-l 

X(f)  =  At  X  u'{n)e 

n-O 


*i2}tfnAt 


n  *= 


TABLE  2 

NOMENCLATURE 


f 

Frequency  (Hz) 

u' 

Velocity  fluctuation  (nVs) 

U 

Mean  axial  velocity  (m/s) 

k 

Wavenumber  (1/m)  =  2Jtf/U 

G(k,U) 

One-sided  power  spectral  density  =  d(u'/U)/dk 

Tu 

Turbulence  Intensity 

N 

Number  of  samples  =  4(X)0 

At 

Time  between  sampling  points  (s) 

X(0 

Fouricr  transfonn 

A 

Integral  length  scale 

Fw  the  first  half  of  the  experiments,  these  calculated  values  for  turbuletre  intensity  and 
power  spectral  density  can  be  compared  with  similar  tests  performed  in  the  f»st  [2].  Frtnn  tftese 
past  experiments,  the  turbulence  intensity  and  power  spectral  density  were  found  to  be 
^fnoximated  by 
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Tu=...3jip 

where  X  and  d  are  the  grid  dimensions  defined  in  Fig.  14,  and 
1  +  (kA)2 

where 


A  =  0.2d 


4.3  Results 

Due  to  the  difference  in  Reynolds  Number  between  the  wind  tunnel  and  blowdown 
turbine,  the  effect  of  Reynolds  Number  on  turlHiient  characreristics  was  investigated.  Figure  15 
shows  the  effect  of  Reynolds  Number  on  power  spectral  density.  From  this  graph,  it  can  be  seen 
that  the  turbulent  characteristics  are  independent  of  Reynolds  Number,  except  for  a  small  influence 
at  high  wave  number. 

The  first  turbulent  characteristic  of  interest  is  the  turbulence  intensity.  Both  the  measured 
and  predicted  values  [1],  [2]  of  turbulence  intensity  over  the  frequency  range  12.5Hz  <  f  <  25,000 
Hz  can  be  found  in  Table  3,  where  25,000  Hz  is  the  Nyquist  frequency. 


TABLE  3 

TURBULENT  INTENSITIES 


Streamwise  Podtion 

No  Nozzle 

Nozzle  1 

1/2  NGV 

NGV 

1/2  NGV 

NGV 

Tu 

23.6 

13.6 

23.0 

6.4 

18.6 

12.2 

18.6 

8.7 

li 


From  Tabic  3,  it  can  be  seen  that  the  measured  values  of  turbulence  intensity  were  higher  than  the 
predicted  values,  except  at  the  NGV  streamwise  location  with  the  nozzle  installed.  This  is  the 
position  of  importance,  because  this  position  conesponds  to  the  positiem  of  interest  in  the 
blowdown  turbine  facility.  At  this  position,  the  measured  value  of  turbulence  intensity  was  lower 
than  the  turbulence  intensity  that  would  occur  at  this  position  in  a  typical  gas  turbine  engine. 

The  other  turbulence  characteristic  of  importance  is  the  power  spectral  density.  Rgurcs  16 
and  17  show  the  graphs  of  the  power  spectral  density  for  the  no  nozzle  case  at  ihc  1/2  NGV  and 
NGV  streamwise  locations  respectively. 

These  two  figures  show  the  measured  values  of  turbulent  energy  at  high  wavenumber 
correspond  nicely  with  theory,  while  the  measured  energy  at  low  wave  number  is  noticeably 
higher  than  predicted.  A  series  of  tests  were  run  in  which  a  low  pass  filter  was  used  to  filter  out 
the  high  frequency  contribution  of  energy  to  the  power  spectral  density  to  show  that  this  high 
energy  level  at  low  wavenumber  was  not  caused  by  aliasing  of  the  signal  when  the  data  was  taken. 
In  Fig.  18,  the  graph  of  the  power  spectral  density  at  the  NGV  streamwise  location  is  plotted  with 
the  nozzle  installed. 

In  Fig.  18,  the  theory  line  [1]  is  measurements  taken  at  the  burner  exit  of  a  simulated  gas 
turbine  combustor.  The  turbulence  energy  measured  in  the  wind  tunnel  matches  the  turbulence 
energy  measured  in  the  gas  turbine  rather  well  at  high  and  low  wavenumber,  however  the  wind 
tunnel  measurements  of  turbulence  energy  are  noticeably  lower  in  the  intermediate  wavenumber 
region.  This  lower  energy  region  agrees  with  the  lower  value  of  tiubulence  intensity  measuml  in 
the  wind  tunnel. 

In  order  to  check  the  effect  of  transverse  positioning  of  the  hot  wire  on  turbulence 
characteristics,  a  number  of  measurements  were  carried  out  at  various  transverse  positions,  while 
keeping  the  streamwise  position  constant.  The  measurements  were  taken  at  the  NGV  streamwise 
posidtHt,  with  nozzle  installed.  Figure  19  shows  the  effect  of  transverse  positioning  in  the  direction 
parallel  to  die  nozzle  by  looking  at  two  positions  1/2  of  one  grid  mesh  apart  in  the  direction  parallel 
to  the  nozzle.  Figure  20  shows  the  effect  of  transverse  positioning  in  the  directicMi  of  nozzle 
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contraction  by  looking  at  two  positions  1/2  of  one  grid  tnesh  apart  in  the  direction  of  nozzle 
contraction. 

Figure  19  shows  that  the  parallel  positioning  of  the  hot  wire  in  the  nozzle  has  little  effect  on 
turbulence  characteristics,  while  Fig.  20  shows  that  hot  wire  position  has  a  considerable  effect  in 
the  direction  of  nozzle  contraction.  When  looking  at  the  turbulence  intensities  at  the  NGV 
streamwise  position,  it  was  found  that  the  turbulence  intensities  varied  in  tlw  dir»^on  of  nozzle 
contraction  due  to  changes  in  the  mean  velocity  across  the  nozzle.  The  toot  mean  square  values  of 
the  velocity  fluctuations  remained  constant  across  the  nozzle.  Figure  21  shows  the  value  of 
turbulence  intensity  in  the  direction  of  nozzle  contraction  at  the  NGV  streamwise  position. 

4.4  Cgndusigns 

The  turbulence  characteristics  of  the  turbulence  grid  to  be  used  in  the  blowdown  turtle 
facility  have  been  presented.  A  number  of  interesting  points  have  been  presented. 

1 .  The  turbulent  characteristics  of  the  flow  produced  by  the  grid  are  independent  of  Reynolds 
Number  in  rite  range  Re  =  5xl(P  to  1.5x10^,  except  far  a  small  effect  at  high  wavenumber. 

2  For  all  cases,  the  measured  value  of  turbulent  intensity  was  larger  than  theory,  except  at  the 
NGV  streamwise  position  with  the  nozzle  installed.  In  this  case,  the  turbulence  intensity  was 
found  to  be  lower  than  that  measured  by  Moss  and  Oldfield. 

3  At  high  and  low  wavenumber,  the  power  spectral  density  was  found  to  match  rather  well  with 
the  results  obtained  by  Moss  and  Oldfield.  At  the  middle  wavenumbers  however,  the  power 
spectral  density  measured  behind  the  turbulence  grids  was  found  to  be  noticeably  less.  This 
result  points  out  the  lower  amount  of  turbulent  energy  contained  in  the  turbulent^  generated 
by  the  grid. 

Assuming  no  change  in  the  velocity  perturbation  through  the  NGV,  the  turbulence  intensi^ 
at  rotCT  entrance  can  be  calculated.  In  the  absolute  finame  of  referaice,  the  turbulence  intensity  was 
calculated  to  be  0.58%,  while  in  the  relative  firame  of  reference,  the  turbuloice  intensity  was 
calculated  to  be  0.97%. 


15 


5.0  Aerodynamic  Instrumentation 

Prior  to  the  stan  of  this  program,  the  aerodynamic  instrumentation  available  on  the 
blowdown  turbine  facility  was  limited  to  the  steady  state  instrumentation  in  Table  4  and  the  time- 
resolved  instrumentation  in  Table  5.  In  the  context  of  the  blowdown  turbine,  steady  state  is  relative 
to  the  rotor  blade  passing  frequency  of  approximately  6  kHz.  For  aero  performance  quality 
measurements,  the  steady  state  instrumentation  requires  frequency  response  of  several  hundred 
hertz  to  follow  the  tunnel  transient  (step  input  startup  followed  by  a  20  seccMid  dnM  constant 
exponential  decay),  and  the  temporal  fluctuations  produced  by  moving  the  measurement  point 
through  a  spatially  nonuniform  flow  field  such  as  a  turbine  exit. 

The  baseline  research  program  added  an  additional  upstream  circumferentially  traversing 
total  temperature  rake  between  the  exit  of  the  distortion  generator  and  the  turbine  inlet,  and  wall 
static  taps  on  the  NGV  inner  and  outer  annulus  at  the  trailing  edge.  After  SNECMA  joined  the 
program,  we  were  able  to  add  turbine  exit  translating  time-averaged  total  pressure  and  temperature 
rakes,  and  a  time-resolved  exit  total  pressure  rake.  After  the  April  program  review  meeting,  an 
upstream  inlet  total  pressure  rake  and  exit  inner  annulus  static  tap  were  added.  Hie  new 
instrumentation  is  summarized  in  Table  6.  In  the  following  sections,  we  discuss  the  design  of  the 
rakes,  the  signal  conditioning  equipment,  and  instrument  calibration  plans. 

5.1  Temperature  Rakes 

The  measurement  of  total  gas  temperature  in  the  blowdown  turbine  represents  unusual 
problems  because  of  the  relatively  fast  time  response  required  (several  hundred  l^rtz)  even  for 
time-averaged  measurenKnts.  The  most  important  design  constraint  is  the  reduction  of  transient 
tiKimal  conduction  error  due  to  heat  conduction  fixim  the  thermocouple  to  the  supporting  structure. 
The  design  of  suitable  probes  was  described  in  some  detail  in  L.  Cattafesta’s  thesis,  in  which 
sevoal  alternate  approaches  were  evaluated. 

The  design  philosophy  was  to  adapt  the  basic  probe  flow  getMnetry  used  at  Rolls-Royce  to 
the  transient  environment  This  has  the  principal  advantage  of  maintaining  continuity  with 


TABLE  4 

BASELINE  BLOWDOWN  TURBINE  STEADY  STATE  AERODYNAMIC  INSTRUMENTATION 


TABLES 

BASELINE  BLOWDOWN  TURBINE  TIME-RESOLVED  INSTRUMENTATION 
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measurements  made  on  the  same  stage  hardware  at  Derby.  To  that  end,  the  basic  probe  is  a 
vented,  shield  Kiel-type  design.  A  “micro-disk”  thermocouple  (20  jim  diameter  x  5  ^im  thick)  is 
used  for  its  high  frequency  response  (300~5{)0  Hz).  The  thermocouple  is  mounted  on  a  75  pm 
diameter,  4  mm  long  quartz  standoff  to  minimize  the  thermal  conduction  from  the  T/C  to  probe 
body.  The  arrangement  is  illustrated  in  Fig.  22.  The  vent  hole  size  is  a  tradeoff  between  a  low 
velocity  need  to  minimize  recovery  factOT  corrections  and  a  high  velocity  which  increases  the  heat 
transfer  to  the  T/C  suppon  stem  and  thus  reduces  the  transient  thermal  error.  Hole  sizes  are  varied 
to  account  for  the  mean  Mach  number  of  the  external  flow  field  and  are  1  mm  and  0.6  mm  on  the 
turbine  inlet  and  exit  probes  respectively. 

These  measurements  reduce  the  transient  thermal  conduction  error  due  to  the  tunnel  startup 
to  1~2°K,  a  level  sufficient  for  the  support  of  heat  transfer  measurements  (150®K  gas-to-wall 
temperature  difference)  but  inadequate  for  aerodynamic  performance  determination.  An  additional 
transient  error  source  is  introduced  when  a  probe  is  travers>>i  through  thermal  gradients  such  as 
wakes.  This  effect  was  simulated  by  driving  a  time  accurate  conduction  model  of  the  temperature 
probe  with  the  output  of  a  2-D  viscous  multi-blade  row  calculaticMi  of  the  flow  in  the  turbine  (Fig. 
23).  A  measured  blowdown  tunnel  transient  is  shown  in  the  temperature  time  history  of  Fig.  24, 
in  which  the  top  trace  (labeled  “gas”)  is  the  output  of  the  micro  disk  thermocouple.  A  second 
thermocouple  was  mounted  at  the  bottom  of  the  75  pm  diameter  quartz  stem  Oabeled  “stem 
bottom”).  The  gas  trace  illustrates  the  rapid  temperature  rise  as  the  tunnel  starts  at  80  ms  and  the 
slow  decay  as  the  supply  tank  blows  down.  The  stem  bottom  trace  reflects  the  teii^ierature  rise  of 
probe  support  structure,  which  is  non-negligible  over  the  first  500  ms  of  interest.  The  differences 
between  the  two  traces  is  a  measure  of  the  driving  temperature  difference  responsible  for  die 
transient  conduction  error.  Given  measurement  of  this  temperature  difference,  the  crniducdon 
error  can  be  estimated  and  a  correction  made.  The  accuracy  to  which  this  procedure  can  be 
implemented  must  be  established  by  calibration  and  will  be  discussed  in  a  later  section.  The 
principal  disadvantage  of  this  approach  is  that  twice  as  many  sensors,  signal  conditioners,  and  data 
recording  channels  are  required  compared  to  just  single  thermocouple  measurements. 
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5.1.1  Inlet  Rake 

The  inlet  rake  consists  of  1 1  Kiel  head  measurement  points  (2  sensors  each)  arranged  to 
sample  equal  areas  in  the  10  cm  high  annulus  (Fig.  25).  The  unit  is  mounted  on  the  inlet 
circumferential  translator,  angled  to  account  for  the  probe  rotation  velocity.  Thermocouple  (Type 
K)  wires  are  derotated  by  a  coil  and  brought  out  through  the  pressure  vessel  wall.  A  “cold” 
junction  ccmnection  to  copper  leads  is  then  made  in  a  dewar.  Rather  than  keep  the  dewar  at  the  ice 
point,  the  cold  block  temperature  is  measured  with  two  4-wire  platinum  resistance  thermometer 
devices  (RTD’s). 

5.1.2  Turbine  Exit  Rake 

The  turbine  exit  rake  consists  of  8  Kiel  head  measuientent  points  (2  sensors  each)  arranged 
to  sample  equal  areas  in  the  3.7  cm  high  annulus  (Fig.  26).  The  unit  mounts  to  the  turbine  exit 
circumferential  translator,  which  presented  two  severe  design  constraints.  The  first  is  that  the  rake 
must  be  readily  removable  from  the  blowdown  tunnel  with  the  rotor  in  place  for  calibration.  The 
only  access  is  through  the  three  small  instrumentation  ports.  This  suggested  placing  the  cold 
junction  within  the  probe  body  itself  so  that  all  thermocouple  wires  and  joints  are  included  in  the 
calibrations.  The  second  constraint  is  space,  since  the  turbine  exit  translator  is  retrofit  to  an  area  in 
which  no  provision  had  been  made  during  the  initial  facility  design.  The  final  probe  design  houses 
a  copper  cold  plate  in  a  cylindrical  housing  below  the  flow  path.  The  therroocoiqile  wires 
transition  to  copper  on  solder  pads  on  the  cold  plate.  The  cold  plate  temperature  is  nwnitored  by 
two  RTD’s.  There  is  sufficient  thermocouple  wire  within  the  rake  body  to  permit  separation  of  the 
rake  head  and  cold  plate  fcH-  calibradon.  The  bottom  of  the  rake  assembly  is  a  41  pin  electrical 
connector  so  that  the  entire  unit  simply  plugs  into  the  exit  translator. 

5.13  Temperature  Rake  Signal  Conditioning 

Each  of  the  rake  sensors  (11  upstream,  8  downstream)  has  two  thomocouples  and  thus 
requires  two  signal  conditioners  and  two  A/D  channels.  The  principal  requirements  for  the 
conditioners  are  low  noise  (compared  to  the  thermocouple  output)  and  long-term  stalnlity  (since 
the  amplifiers  must  remain  stable  between  calibrations,  hopefully  several  weeks).  Analog  devices 
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2B31H  instrumentation  signal  conditioners  were  selected.  They  are  operated  at  a  gain  of 
approximately  800  and  the  internal  two  pole  filters  are  set  at  500  Hz.  The  same  units  are  also  used 
as  RTD  signal  conditioners  for  the  cold  junction  references.  Fony-two  channels  have  been  built. 

The  data  acquisition  system  has  been  expanded  with  the  addition  of  an  Analogic  64- 
channel,  16-bit,  PC-based  A/D  subsystem.  The  maximum  aggregate  sampling  rate  for  the  unit  is 
200  kHz.  Software  has  been  written  so  that  this  unit  mimics  the  functionality  of  the  existing 
blowdown  turbine  data  system  (41  channels  at  2(X)  kHz  per  channel  plus  128  channels  at  16  kHz 
may  rate  per  channel).  The  new  unit  is  required  because  of  inadequate  stability  for  aero 
performance  of  the  older  low  speed  channels  as  well  as  lack  of  capacity. 

5.1.4  Temperature  Rake  Calibration 

Three  calibration  steps  are  required  to  establish  the  rake  accuracy.  The  first  is  a  steady  state 
voltage  output  versus  temperature  calibration.  This  will  be  performed  in  a  stirred  silicon  oil  bath 
over  the  SOO-SOO^K  range  of  interest.  A  secondary  standard  (0.0 1®K)  platinum  RTD  is  the 
temperature  reference.  The  second  step  is  to  establish  the  long-term  stability  of  the  sensors,  signal 
cemditioners,  and  A/D,  which  will  be  established  by  multiple  calibrations.  The  third  and  final  step 
is  determining  the  overall  accuracy  of  the  rake  in  measuring  gas  temperature,  which  will  include 
the  transient  conduction  and  recovery  factor  corrections. 

The  final  aerodynamic  calibration  of  the  rakes  requires  that:  (a)  the  freestream  ten^rature 
be  uniform  and  known  to  the  necessary  accuracy  (0.1~0.3®K);  (b)  the  freestream  Mach  number  be 
variable  to  match  the  turbine  inlet  and  exit  conditions;  and  (c)  the  flow  turn  on  with  approximately 
the  same  time  constant  as  the  blowdown  tunnel  (30~50  ms).  To  accomplish  this,  we  have  built  a 
75  mm  diameter  blowdown  test  tunnel  using  blowdown  turbine  components  (Fig.  27).  The 
heated  gas  supply  is  the  blowdown  turbine  supply  tank  (-1 1  volume).  The  flow  is  turned  on 
and  off  using  the  75  mm  diameter,  fast-acting  (15  ms)  coolant  supply  valve.  An  orifice  at  the  test 
section  discharge  is  used  to  establish  the  Mach  number.  The  test  tunnel  is  closely  coupled  to  the 
supply  tank  to  keep  the  thermal  boundary  layers  thin  at  the  probe  location.  The  supply  tank  gas  is 
homogenized  prior  to  the  test  by  three  fans.  The  temperature  distribution  is  established  by  1 1 


platinum  RTD’s  spotted  through  the  tank  (Fig.  28).  The  RTD’s  are  bath-calibrated  against  a 
secondary  standard  before  installation.  Figure  29  is  data  taken  from  the  tank  sensors  and  shows 
that  the  temperature  is  uniform  to  about  0.2°K. 

The  tank  RTD’s  and  fans  have  been  removed  from  the  supply  tank  for  the  turbine  testing. 
The  plan  is  to  bath  calibrate  the  temperature  rakes,  run  the  turbine  tests,  and  then,  in  the  late 
summer-early  fall  time  frame,  set  up  the  test  tunnel  and  run  the  calibratitm  checks. 

5.2  Total  Pressure  Rahes 

Three  translating  total  pressure  rakes  are  planned  for  these  tests:  two  exit  rakes  described  at 
the  April  review  meeting,  and  a  new  inlet  rake  as  suggested  at  the  meeting.  The  inlet  rake  and  one 
exit  rake  are  designed  to  measure  time-averaged  total  pressure  and  look  externally  the  same  as  the 
total  temperature  rakes  in  Figs.  25  and  26.  They  are  vented  Kiel  head  designs  with  short  (75~100 
mm)  pressure  tubes  leading  to  transducers  mounted  in  the  rake  housing  under  the  flow  path.  This 
is  to  insulate  the  pressure  transducers  from  the  temperature  transient  of  the  tunnel  (^ration  and 
maintain  stability  over  the  first  SOO  ms  to  the  end  of  the  test  time.  One  Kulite  X(I!Q-063 
transducer  is  used  per  rake  position  (1 1  at  inlet,  8  at  outlet).  The  pressure  tubes  extend  to  the 
transducer  faces,  so  that  there  is  no  “dead”  volume  and  thus  frequency  response  is  maximized. 

The  differential  transducers  are  used  with  the  reference  tube  brought  out  to  the  laboratory.  During 
the  test  time,  the  transducers  are  referenced  to  vacuum.  Just  prior  to  and  after  each,  the  transducers 
are  calibrated  by  changing  the  reference  pressure.  Thus,  the  transducers  need  not  remain  staUe  for 
more  than  a  minute  or  two.  This  results  in  stable  and  accurate  DC  pressure  readings.  The 
pressure  rakes  are  interchangeable  with  the  temperature  rakes. 

The  third  pressure  take  is  designed  to  measure  time-resolved  total  pressure  so  the 
transducers  (Kulite  XCQ-063)  are  mounted  in  the  eight  rake  Kiel  heads,  protected  by  a  screen. 
Because  the  transducers  are  now  exposed  to  the  flow,  their  temperature  change  and  thermal  drift 
are  expected  to  be  greater  than  for  the  takes  with  the  protected  and  insulated  transducers.  If  drift  is 
a  problem,  then  the  DC  values  for  these  transducers  will  be  corrected  with  the  readings  from  the 
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Steady-State  rake. 

5.3  Static  Pressure  Taos 

Wall  static  pressure  taps  have  been  added  to  the  inner  and  outer  annulus  walls  at  the  NG  V 
trailing  edge  (mid-passage),  and  at  the  inner  annulus  at  the  turbine  exit  translator  at  the  rake  axial 
station.  All  three  consist  of  short  tubes  (50~1(X)  mm)  fed  through  Kulite  XCXJ-093  transducers. 

6.0  Rotor  Blade  Heat  ¥\ux  Instrumentation 

Since  a  principal  output  of  this  research  program  is  a  set  of  detailed  rotor  blade  heat  transfer 
measurements,  the  rotor  blade  instrumentation  is  the  single  most  important  pan  of  the  facility 
preparation.  It  is  also  the  most  time  consuming  and  expensive  by  far,  with  more  than  one  half  the 
program  resources  to  date  being  expended  in  this  area.  For  this  program,  the  instrumentation  has 
been  entirely  redone,  including  the  manufacture  of  transducers,  their  mounting  and  wiring,  and 
blade  calibration. 

Three  blades  have  been  produced,  each  instrumented  about  their  chtxd  with  the  spanwise 
location  differing  among  blades.  A  composite  of  the  three  blades  is  shown  in  Fig.  30,  which 
illustraes  the  heat  flux  gauges’  locations.  Thirty-nine  gauges  are  fully  functioning.  An  additional 
five  have  just  the  t(^  sensor  q>erational,  in  which  case  essentially  all  information  can  be  retrieved 
with  some  additicHial  data  reduction  effort  On  three  additional  sensex^,  only  the  bottom  sensors 
are  working,  in  which  case  little  information  can  be  extracted.  Two  additional  gauges  are  mounted 
on  the  blade  platform  near  midpassage  at  the  leading  edge  plane.  Resistance  thermometers  are 
mounted  to  the  underside  of  the  blade  platforms  to  serve  as  pre-  and  post-test  temperature 
references  which  facilitatt  scale  factor  and  drift  determination. 

The  wiring  for  each  of  the  three  blades  is  shown  schematically  in  Fig.  31.  The  signals 
from  each  transducer  is  brought  out  by  slip  rings.  The  signals  are  then  fed  through  34  bridge 
an^lifiers  (DC-SO  kHz  bandwidth)  and  into  individual  data  acquisiticxi  channels  (200  kHz 
sampling  rate  for  the  top  sensors,  16  kHz  for  the  bottom  sensors). 

The  blades  have  been  completed  and  the  94  amplifiers  required  checked  out.  Blade 
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calibration  will  be  done  as  the  NGV-only  testing  is  ongoing  in  order  to  free  up  manpower  for 
readying  the  facility. 

7.0  CFD  Activities 

The  principal  CFD  tool  to  be  used  at  MIT  on  this  program  is  the  thiee-dimensicmal, 
unsteady,  multi-blade  row  Euler  code  written  by  Saxer  and  Giles,  UNSFL03.  This  is  an  explicit 
time-marching,  node-based,  Ni-Lax-Wendroff  algorithm  implemented  on  an  unstructured  grid. 
Non-reflecting  boundary  conditions  are  applied  at  the  inlet  and  outlet  and  at  the  statesr-rotor 
interface. 

This  code  had  been  run  on  the  experiment’s  NGV  and  rotor  blade  geometries  by  Saxer  but 
at  equal  NGV-rotor  pitches  rather  than  the  36  NGV’s,  61  rotor  blades  of  the  rig.  The  grid 
geometry  has  been  modified  to  calculate  3  NGV  passages  and  5  rotor  passages  (i.e.,  36  NGV,  60 
rotOT  blade  count)  as  illustrated  in  Fig.  32.  This  grid  is  currently  set  up  with  about  256,000  node 
points  (56  x  36  X  21  in  each  of  3  NGV  passages  and  56  x  22  x  21  in  each  of  5  rotor  passages).  A 
coarser  version  of  the  grid  (only  9  rather  than  21  radial  planes)  has  been  run  to  convergence  as  an 
initial  checkout.  Rotor  blade  temperature  contours  are  shown  in  Fig.  33.  The  finer  grid  is 
currently  running. 

The  inlet  boundary  conditions  have  been  set  up  so  that  quite  general  radial  temperature  and 
pressure  distortions  can  be  specified.  Periodicity,  however,  constrains  hot  spot  placement  to  one 
every  3  NGV  pitches  where  the  experiment  is  set  up  as  one  hot  spot  for  every  2  NGV  pitches 
(Fig.  34).  After  we  gain  some  «cperience  with  the  code,  the  blade  count  can  be  enlarged  to  6 
NGV’s  and  10  rotor  blades  to  permit  the  experiment’s  poiodicity  to  be  modelled.  This  is  at  the 
expense  of  doubling  the  run  time,  of  course.  The  grids  have  been  interfaced  to  die  MIT  3-D 
visualization  codes. 

In  addition  to  UNSFL03,  the  two-dimensional  Euler  version  aind  a  2-D  thin  shear  layer 
viscous  version  are  available  at  MIT  and  have  been  checked  out  with  the  experimental  geometry. 
Also  available  is  John  Denton’s  MULTISTAGE  3-D,  steady,  multi-blade  row,  Euler  code. 


25 


Discussions  arc  ongoing  with  NASA  Lewis  Research  Center  about  using  Jdin  Adainczyk's  3-D 
viscous  code. 

8.0  Rig  Preparations 

Rig  preparations  included  completion  and  insialladon  of  the  new  rake  transkicai^ 
lefurbishment  and  repair  of  tunnel  subsystems,  installation  and  programming  of  new  dau 
acquisition  equipment,  and  instrumentation  calibrations  and  cables.  Tla:  inlet  and  ouitet  rake 
translators  were  <x>mp!ctcd  and  installed  in  the  rig.  The  units  performed  as  designed  during 
vacuum  testing. 

Several  of  the  rig  subsystems  required  refurbishment  and  repair,  which  was  not  surprising 
considering  that  the  rig  is  10  years  old  and  had  not  been  operated  for  about  3  years.  The  Freon 
boiler  in  the  gas  mixing  system  failed  pressure  testing  and  was  replaced  as  were  numerous  leaking 
valves.  The  vacuum  system  valves  and  vacuum  pumps  were  overhauled.  When  the  test  section 
was  disassembled  to  mount  the  rotor  exit  translator,  it  was  discovered  that  a  wire  on  one  of  the 
eddy  t»ake  magnet  coils  was  broken.  This  explains  the  sudden  shift  in  brake  calibration  noted 
during  the  previous  test  series.  The  eddy  brake  was  removed  for  repair  and  will  be  re-installed 
with  the  rotor. 

8.1  Data  Acquisition 

Two  new  data  acquisition  systems  have  been  added  to  supplement  the  current  unit  and 
serve  as  backup  (the  current  system  is  1 1  years  old).  Sixty-four  low  speed  channels  have  been 
added.  These  units  have  16-bit  resolution  (about  14  of  which  are  usable  given  the  noise  levels  at 
the  rig),  and  3  kHz  per  channel  maximum  data  rate.  The  unit  is  mounted  in  a  ‘486  PC  These  low 
iKHse,  high  resoludrai  channels  have  been  alk)cated  i»imarily  to  the  ^poatuit  rake  raiqnits. 
Additional  high  bandwidth  channels  have  been  added  as  well.  Thirty-two  channel  (24  arrived.  8 
CHI  CHder),  12-bit,  330  kHz  per  channel  A/D’s  have  been  added  in  a  second  ‘486  PC.  This  unit  was 
purchased  by  the  laboratory  and  will  be  shared  with  other  experiments.  This  brings  the  total 
channel  capacity  of  the  blowdown  turbine  facility  to  approximately  68  high  speed  analog  channels 
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(200+  kHz  per  channel)  and  192  low  speed  analog  channels  (3+  kHz  per  channel). 

Both  the  new  low  and  high  speed  data  acquisition  systems  (A/D’s  plus  computers) 
required  considerable  custom  cabling  and  programming.  The  units  were  int^;rainm6d  to  interface 
to  the  current  data  handling  systems  and  are  mostly  transparent  to  the  user. 

8.2  ln.strumgntation  Calihratinn 

Upon  completion  of  the  A^D  systems,  the  temperature  and  pressure  rakes  were  DC 
calibrated.  The  temperature  rakes  are  of  a  pitot-type  configuration  with  two  type  K  thermocouples 
at  each  head  position,  one  to  measure  gas  temperature,  the  other  to  measure  the  temperature  of  the 
new  stem  on  which  the  first  is  mounted,  to  permit  corrections  for  conduction  errors.  The  rakes  arc 
designed  with  local  reference  hot  junctions  rather  than  standard  remote  cold  junctions  (Fig.  35a). 
This  arrangement  facilitates  rapid  recalibration  of  the  rakes  without  necessitating  teardown  of  the 
test  section,  as  would  be  the  case  if  thermocouple  leads  were  run  out  to  external  cold  junctions. 
(The  rakes  can  be  removed  through  windows  in  the  test  secticHi  while  leaving  the  cabling  in  place.) 
The  hot  junctions  consists  of  an  insulated  copp^  block  on  which  the  thermocouple  to  copper  wire 
connections  are  made.  The  block  is  instrumented  with  two  resistance  thermometers  (RTD’s). 
Calibration  consisted  then  of  two  steps  -  calibration  of  the  reference  junction  sensors  followed  by 
calilnation  of  the  thermocouples.  Both  were  done  in  a  stirred  bath  of  silicon  oil  against  a 
secondary  standard  platinum  RTD  stable  to  O.OPC.  The  calibrations  were  performed  through  the 
entire  data  chain  -  sensor,  cabling,  signal  conditioners,  data  acquisition  systenL  Transducers 
agreed  with  each  other  to  O.l^C.  Repeatability  of  calibration  was  within  0.2°C. 

The  pressure  transducers  in  the  rakes  (Fig.  3Sb)  are  readily  calibrated  since  each  is  a  differ¬ 
ential  transducer  with  a  reference  pressure  which  is  brought  out  of  the  rake  and  through  the  transla¬ 
tor.  In  practice,  the  transducers  are  calibrated  in  place  immediately  before  and  after  each  test  The 
pressure  transducers  were  checked  for  ten^ierature  sensitivity  by  calibrating  them  in  an  oven  at 
temperatures  up  to  100°C.  All  transducers  are  stable  to  better  than  1%  over  the  temperature  range. 
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8.3  Instrumentation  Performance  on  Translators 

The  rake  instrumemation  was  installed  in  the  translators  (Fig.  35c)  which  were  then  run  in 
vacuum  as  an  operational  check.  The  pressure  rakes’  operation  was  satisfactory  in  that  no 
extraneous  signal  or  noise  was  observed  on  the  pressure  transducer  outputs  during  the  translatcw 
motion  (Fig.  6).  The  temperature  outputs,  however,  did  show  artifacts  due  to  both  translator 
motion  and  operation  of  the  DC  servo  systems  driving  the  two  translators.  A  signal  related  to  the 
motion  can  be  seen  on  the  output  from  an  inlet  rake  thermocouple  in  Fig.  36,  starting  at  about  200 
msec  (same  time  scale  as  in  Fig.  6).  The  signal  is  similar  in  shape  and  magnitude  on  all  22 
thermocouples  on  the  inlet  rake.  The  signals  are  clearly  related  to  the  probe  nxstion  in  that  they  do 
not  appear  on  the  reference  RTD  signals  (Fig.  37)  which  arc  mounted  in  the  stationary  frame  for 
the  inlet  rake.  The  magnitude  of  this  unwanted  signal,  about  0.2®C,  is  sufficiently  low  that  it  does 
not  serve  as  a  hindrance  to  this  research  program  but,  at  this  level,  it  is  a  major  component  of  the 
upstream  gas  temperature  measurement  uncertainty. 

The  level  of  unwanted  signal  is  much  greater  on  the  outlet  temperature  rake,  as  can  be  seen 
in  Fig.  38,  in  which  the  ‘noise’  amplitude  is  about  3°C.  As  with  the  inlet  rake,  the  signal  is  similar 
on  all  the  thermocouples  on  the  rake  to  a  remarkable  degree  (Fig.  39).  Unlike  die  inlet  rake, 
however,  the  outlet  temperature  reference  RTD’s  are  on  the  rotating  system  and  these  exhibit  a 
signal  amplitude  similar  to  that  of  the  thermocouples.  Fig.  40.  (The  plotted  ‘reduced’ 
thermocouple  temperature  traces  incorporate  the  RTD  signals  since  these  arc  reference  sources 
which  are  differenced  with  the  raw  thermocouple  temperature.) 

The  level  of  temperature  uncertainty  represented  by  the  level  of  interfering  signal  on  the 
outlet  ttmiperature  rake  (3°C)  is  clearly  unacceptable  for  aerqierformance  measurements  and 
marginal  for  reducing  heat  transfer  measurements.  Therefore,  considerable  effort  was  put  into 
understanding  this  problem.  Sources  of  the  interference  considered  included  some  combination  of 
magnetic  pickup  on  the  rotating  system,  cross-talk  between  layers  of  cabling,  grounding  problems, 
etc.  The  presence  of  this  signal  on  the  outlet  RTD’s  (which  were  operated  in  a  4-wirc 
configuration  to  eliminate  common  signals  induced  in  the  leads)  and  its  absence  on  the  pressure 
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rakes  (4-wirc  bridge  configuratior>)  which  share  the  same  cabling  was  particularly  puzzling. 
Comparison  of  these  signals,  as  in  Fig.  41,  was  interesting  but  unenlightening.  Circuit  modelling 
and  bench  testing  were  unsuccessful  in  reproducing  or  explaining  the  problem.  Mu  metal 
shielding  was  added  to  the  cabling  to  exclude  magnetic  fields  with  no  avail.  Finally,  the  grounding 
scheme  of  the  rakes  was  modified  in  an  od  hoc  manner.  This  succeeded  in  reducing  the  induced 
signal  on  the  exit  rake  by  more  than  an  order  of  magnitude  (Figs.  42  and  43),  to  less  than  0. 15®C. 
Only  minor  improvement  was  observed  in  the  upstream  rake  ou^ut  (Fig.  44).  At  its  current  level, 
these  unwanted  signals  will  not  adversely  affect  the  current  research  program.  We  fixed  things 
without  understanding  why  and  will  now  retire  from  the  field  of  battle  clothed  in  graceful  victory 
and  ignorance. 

9.0  Initial  NGV-Onlv  Blowdown  Testing 

Two  blowdown  tests  were  run  within  the  NGV-only  configuration  -  with  the  rotor 
removed,  and  with  the  outlet  rakes  moved  upstream  to  the  NGV  exit  (Fig.  45).  Several  problems 
were  encountered  during  the  runs  but  data  was  nevertheless  taken.  In  the  following  sections,  we 
will  first  examine  the  data  and  then  discuss  the  problems  which  arose  and  their  resolution. 

9.1  Initial  NGV  Test  Data 

The  primary  purpose  of  these  tests  is  ch^kout  and  evaluation  of  the  facility  and  its 
instrumentation.  Two  tests  were  run  -  test  150  with  both  radial  and  circumferential  distortions, 
and  test  151  with  radial  distortions  only  (the  circumferential  distortion  heat  burned  out,  as  will  be 
discussed  later).  A  pressure  history  from  the  inlet  total  pressure  rake  (Fig.  46)  illustrates  tunnel 
operation.  The  mainflow  valve  is  commanded  to  open  at  0  milliseconds  (ms).  The  pressure 
begins  to  rise  in  the  test  section  at  60  ms  and  peaks  at  95  ms.  The  20  second  e}q)aimental  decay 
time  constant  of  the  facility  then  can  be  seen  in  the  decay  of  the  mean  pressure  tevel  from  100  to 
600  ms.  Superimposed  on  this  is  a  ringing  of  the  cavity  formed  by  the  inlet  valve  to  the  NGV 
throat,  driven  by  the  startup  fast  rise.  This  ringing  has  substantially  dried  out  by  250  ms,  the  start 
of  the  useful  test  time  which  runs  until  550  ms  into  the  test.  (The  data  acquisition  systems  operate 
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at  their  highest  rate,  200  kHz  per  channel,  during  this  250  to  550  test  time.)  This  time  traux  also 
illustrates  the  radial  flow  uniformity  at  the  tunnel  inlet  and  thermal  drift  limits  of  the  transducer. 

All  six  of  the  radial  stations  agree  to  within  one  percent,  and  five  of  the  six  to  0.5%.  We  believe 
that  the  sixth  transducer’s  performance  can  be  improved  with  additional  calibration. 

Performance  of  the  low  fi'cquency  response  turbine  exit  pressure  rake  and  the  wall  static 
taps  is  shown  in  Fig.  47.  The  sawtooth  signals  are  the  NGV  wakes  through  which  the  rake  is 
being  translated.  The  sudden  steps  in  two  signals  are  due  to  a  wiring  failure  which  will  be 
discussed  later.  An  example  of  the  pressure  rake  signal  during  the  250-550  ms  test  time  can  Le 
seen  in  Fig.  48  in  which  the  NGV  exit  flow  structure  can  be  clearly  seen  superimposed  on  the  DC 
pressure  decay  of  the  facility.  A  similar  plot  is  shown  in  Fig.  49  for  a  high  frequency  response 
total  pressure  rake  midspan  position.  We  infer  from  an  examination  of  this  data  that  the  pressure 
rake  performance  is  satisfactory  and  as  planned. 

The  response  of  the  upstream  temperature  rake  during  the  circumferential  traverse  in  a 
blowdown  test  is  illustrated  in  Fig.  50  for  the  midspan  location.  The  initial  upspike  at  90-100  ms 
is  compressional  heating  during  the  startup  transient.  The  three  downspikes  at  250  ms,  390  ms, 
and  550  ms  are  the  thermal  wakes  of  the  stmts  separating  the  three  sectors  of  the  temperature 
distortion  generator.  The  data  from  all  eleven  sensors  along  the  span  was  used  to  construct  the 
contour  plot  of  gas  temperature  in  Fig.  51.  The  circumfoentia^  translator  is  accelerating  until  about 
120  ms  (resulting  in  the  close  contour  lines),  but  is  travelling  at  constant  velocity  during  the  250- 
550  ms  test  time.  Thus,  the  time  axis  can  be  considered  to  be  circumferential  position  during  the 
test  time.  That  period  replotted  in  Fig.  52  shows  the  traverse  of  two  120°  segments  of  the  annulus, 
a  0TDF  sector  from  250  to  390  ms,  and  an  RTDF  sector  from  390  to  550  ms.  The  stmt  wake  can 
be  seen  at  390  ms.  The  radial  temperature  distention  is  seen  to  be  uniform  circumferentially 
between  450  and  510  ms  (about  a  45  degree  sectw  or  4-5  NGV  spacings). 

The  output  from  the  midspan  exit  total  temperature  rake  sensor  as  it  is  traversed  about  the 
circumference  at  the  NGV  exit  is  illustrated  in  Fig.  53.  The  compressional  startup  spike  is  evident 
at  90  ms.  The  cold  wakes  of  the  NGV’s  can  be  seen  between  200  and  300  ms.  From  300  to  500 
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ms,  the  rake  is  traversing  through  the  6TDF  sector,  whose  heater  was  burned  out  during  this  test, 
so  that  the  temperature  is  quite  low.  After  examining  the  temperature  data,  we  conclude  that  the 
temperature  instrumentation  is  operating  as  designed  and  that  its  performance  is  satisfacmry  for 
this  research  project. 

9.2  Exit  Rake  Translator  Wiring  Problems 

During  the  Hrst  NGV-only  blowdown  tests,  the  wiring  which  runs  from  the  translator  to 
the  stationary  frame  was  seriously  damaged.  The  wiring  consists  of  three  SO-conductw  ribbon 
cables  sandwiched  between  four  layers  of  75  pm  stainless  steel  shim  spring  stock.  The  multi¬ 
layer  cable  is  3/4  of  a  mt.ter  long  and  arranged  in  a  hairpin  with  one  end  fastened  to  the  rotating 
translatca*  and  the  other  fixed  to  the  stationary  tost  section.  The  cable  is  located  at  the  rear  of  the 
turbine  exit  flow  path  in  a  protective  shroud.  This  loop  arrangement  (Figs.  8  and  54)  facilitates 
360  degree  travel  during  the  test  The  shims  both  protect  the  ribbon  conductors  aixi  provide  the 
restoring  force  needed  when  the  translator  is  reset  between  tests.  Unfortunately,  the  ribbon  blew 
out  under  its  containment  shroud  and  became  entangled  in  the  chain  drive  (Fig.  55)  approximately 
250  ms  into  the  first  test 

Several  solutions  were  considered,  including  a  rotary  seal  on  the  protective  shroud  and 
various  modificatipns  to  the  laminated  ribbon  assembly.  The  geometric  constraints  of  the  test 
section  proved  very  restricting.  A  full-scale  mockup  was  constructed  to  evaluate  various  schemes. 
The  approach  adt^ted  was  (a)  to  stiffen  up  the  ribbon  assembly  by  replacing  the  outer  two  75  pm 
thick  shims  with  150  pm  shims,  and  (b)  the  addition  of  a  rotating  containment  ring  at  the  inner 
radius  of  the  stationary  prot^ve  shroud.  This  ring  increases  the  moment  of  inertia  of  die  rotating 
system,  slowing  the  translator  acceleration. 

10.0  Temperature  Distortion  Generator 

The  temperature  distortion  generator  used  for  these  tests  is  the  honeycomb  storage  matrix 
unit  originally  designed  for  RTDF  tests.  It  uses  electric  heating  wires  at  midspan  as  a  heat  source 
and  oil  cooling  jackets  on  the  inner  and  outer  flow  path  walls  to  serve  as  heat  sinks,  thus  generating 
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a  parabolic  temperature  profile  aaoss  the  annulus  during  the  pretest  heatup  period.  The  unit  was 
modified  for  6TDF  testing  by  spacing  the  heating  wires  along  the  circumference  in  raie  120® 
sector  to  generate  hot  spots,  one  hot  spot  per  two  NGV  pitches  (Fig.  1). 

Previous  testing  had  shown  that  the  generator  was  contaminated  with  fine  particles  during 
manufacturing  which  destroyed  the  rotor  blade  heat  flux  instrumentation.  TTiereforc,  a  1(XX)  liter 
ultrasonic  cleaner  was  constructed  so  that  the  generator  could  be  immersion  cleaned.  Afta*  each 
cleaning,  the  water  was  drained  and  filtered  and  the  particles  counted.  The  cleaning  was  repeated 
until  no  particles  could  be  discerned. 

10.1  Girciimfergnfial  Distortion  Heatup  Difficnlfigs 

The  design  of  the  distortion  generator  was  based  on  the  original  design  estimates  of  the 
thermal  diffusivity  of  the  stainless  steel  honeycomb  matrix.  Variations  in  hot  spot  amplitude  and 
extent  were  to  be  achieved  by  adjusting  the  current  time  history  into  the  electric  heaters  as 
illustrated  in  Fig,  2.  During  the  first  tests,  it  was  observed  that  circumferential  temperature 
gradients  were  small  relative  to  those  in  the  radial  direction.  Thermocouples  spotted  about  the 
matrix  permitted  explicit  measurement  of  relative  thermal  diffusivities  (Fig.  56).  Figure  57  shows 
the  difference  between  the  radial  and  circumferential  temperature  gr^ents  to  step  changes  in  heat 
input.  The  matrix  equilibrates  in  the  radial  direction  in  about  an  hour,  while  only  five  minutes  is 
needed  in  the  circumferential  direction.  This  implies  that  circumferential  distortions  will  quickly 
“wash  out”  compared  to  radial  ones. 

To  quantify  this  difference,  a  two-dimensional  transient  conduction  model  of  the  distortion 
generator  was  written.  The  diomal  diffusivities  in  the  radial  and  circumferential  directkms  were 
estimated  by  fitting  the  model  to  time-temperature  histories  of  thermocouples  spotted  about  the 
matrix.  The  model  fits  the  data  well  (Fig.  58),  giving  ctMifidence  in  the  results.  The  conclusion  is 
that  the  thermal  conductivity  is  15  times  greater  in  the  circumferential  direction  than  in  the  radial 
one.  The  cause  of  this  large  a  nonuniformity  is  not  clear  (poor  brazing?)  Ixit  the  implication  is  that 
0TDF  can  be  held  for  only  short  periods  of  time  and  titat  large  amplitude  hot  spots  will  be  difficult 
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to  achieve.  This  is  illustrated  in  Fig.  59,  which  shows  the  result  of  a  transient  h^t  calculation  of  a 
sector  with  two  hot  spots.  At  heater  current  levels  compatible  with  past  practice,  a  60  to  80°C  hot 
spot  is  achieved  in  two  minutes.  After  an  additional  two  hours  of  heating,  tftesc  relatively  small 
hot  spots  are  now  superimposed  on  a  large  radial  temperature  profile.  This  pattern  is  quite  similar 
to  that  in  an  engine  (good)  but  not  so  convenient  as  pure  6TDF  for  simple  analysis  (bad). 

One  solution  to  this  problem  is  very  fast,  high  current  spikes  applied  immediately  (l*-3 
minutes)  before  the  start  of  a  test.  The  constraint  is  the  current  level  which  can  be  put  into  the 
heater  before  mechanical  failure  (burnout)  of  either  the  heater  wire  or  its  termination.  This  was 
investigated  experimentally  on  a  test  rig  in  vacuum,  and  it  appears  that  6TDF’s  of  20%  can  be 
generated  before  either  the  heater  wire  bums  out  or  the  brazed  matrix  is  damaged. 

The  solution  adopted  was  to  double  the  number  of  heater  wires  at  each  spot  from  2  to  4, 
thus  doubling  the  heat  input.  The  6TDF  wires  are  operated  at  currents  of  30  amps  (as  opposed  to 
10  amps  in  the  RTDF  sector),  increasing  the  heating  rate  by  nine  times.  This  high  heat  rate  is 
applied  five  minutes  immediately  prior  to  the  blowdown  test  and  should  pioduce  hot  spots  of 
approximately  50®C.  Vacuum  tests  with  thermocouples  in  the  matrix  indicated  this  to  be  a  viable 
approach.  As  of  this  time,  one  half  pressure  test  has  been  run  at  low  current  (20  amps)  as  a 
mechanical  checkout  Its  results  are  inconclusive  insofar  as  6TDF  levels  achievable. 

10.2  Distortion  Generator  Contamination 

Although  the  generator  was  thoroughly  cleaned  before  testing,  aluminum  witness  cylinders 
were  placed  in  the  flow  path  at  the  NGV  exit  Microscopic  examination  after  the  first  test  showed 
significant  impact,  equivalent  to  about  10^  particles  per  test  if  the  sampled  area  is  extrapolated  over 
the  entire  flow  annulus.  This  level  of  contamination  is  unacceptable. 

Particles  were  collected  and  examined.  Two  types  were  noted  -  clear  spheres  of  about  50- 
100  pm  diameter,  and  brown  cubes  of  similar  size.  Since  clear  spheres  did  not  dissolve  in  any 
acid  other  than  hydrofluoric,  they  are  thought  to  be  glass  sandblast  beads.  The  brown  particles  did 
n(K  dissolve  in  any  acid.  They  closely  resemble  and  are  themght  to  be  aluminum  oxide  abrasive 
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particles,  also  used  for  sandblasting. 

The  question  is,  of  course,  how  to  remove  the  particles.  The  ultrasonic  cleaning  removed 
all  the  particles  that  ultrasonic  cleaning  will  remove.  What  to  do  about  the  rest?  Mechanical 
cleaning  techniques  considered  included  shaking  the  unit  on  a  100  g  shaker  used  for  satellite 
prelaunch  testing,  repeated  heating  and  blowing  cycles,  and  steam  cleaning.  Chemical  means  of 
dissolving  the  particles  were  also  considered.  Tests  indicated  that  the  procedure  outlined  in  Table  7 
might  be  effective.  The  question  with  all  of  these  alternatives  is  how  to  remove  the  panicles  but 
not  damage  the  brazed  stainless  steel  matrix.  Since  we  know  that  the  matrix  biazing  is 
substandard,  all  of  these  options  entail  considerable  risk  of  damaging  the  matrix. 

TABLE  7 

ONE  CHEMICAL  CLEANING  OPTION 

High  pressure  steam  clean,  followed  by  H2O 
8-hour  immersion  in  ammonium  bifluoride  and  HCL 
8-hour  immersion  in  sulfuric  acid  at  140®F 
Neutralize  in  NaOH 
High  pressure  steam,  followed  by  H2G 
_ Clean  gas  flow _ 


The  solution  of  replacing  the  distortion  generator  with  one  01  a  different  design  or  from  a 
more  reliable  vendor  is  not  a  serious  option  due  to  cost  and  schedule  constraints.  The  solution 
adopted  was  to  fabricate  a  particle  filter  for  the  distcfftion  generator  option.  The  engineering  design 
criteria  for  the  filter  includes  filter  efficiency  (99%+  of  particles  50  iim  diameter  and  larger),  low 
pressure  drop  at  25  kg/sec  mass  flow,  low  thermal  inertia  (to  prevent  attenuation  of  the  distortion 
pattern),  mechanical  integrity,  and  rapid,  relatively  low  cost  fabrication.  The  final  design 
incorporates  a  sintered  stainless  steel  filter  mesh  rated  as  a  nominal  20  pm  particle  filter,  backed  up 
by  a  45%  porosity,  three  millimeter  thick  poforated  stainless  steel  plate  for  mechanical  support 
The  filter  in  the  rig  is  illustrated  in  Fig.  7.  Because  the  filter  extwids  the  rig  length  by  about  3  cm, 
realignimnt  of  the  supply  tank  and  modification  of  several  auxiliary  systems  was  required. 
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At  this  time,  one  blowdown  test  has  been  run  at  one-half  pressure  with  the  filter  in  place. 
No  mechanical  problems  were  encountered.  The  trapping  of  particulates  was  evident  when  the 
filter  was  disassembled  for  inspection  after  the  test.  Thus,  this  problem  appears  solved. 

11.0  Gas  AYaUabilitx.and  Test  Costs 

Since  the  blowdown  turi>ine  facility  was  last  operated,  the  cost  of  Freon- 12  has  risen 
tenfold,  fi’om  $l/lb  to  $10/lb.  With  allowance  for  leakage  and  aborted  tests,  this  translates  to 
$1200  per  test  for  Freon,  a  cost  level  which  is  not  in  the  current  budget  There  are  several  possible 
alternative  approaches  to  this  problem,  including: 

a)  Running  fewer  tests 

b)  Testing  with  the  Freon- 1 2  replacement,  HC- 1 34a 

c)  Recycling  the  gas  mixture 

d)  Testing  in  air  at  y  =  1  -4 

e)  Testing  with  an  argon  and  CXD2  mixture  at  the  proper  y 

f)  Some  combination  of  the  above. 

Each  of  these  approaches  has  some  advantages  and  serious  disadvantages  which  are  discussed 
below. 

(a)  Obviously,  running  fewer  tests  is  the  most  unattractive  since  the  entire  point  of  the 
project  is  to  take  data;  the  data,  however,  must  be  relevant,  i.e.  properly  scaled  conditions. 

(b)  The  Freon- 12  replacement,  He- 134a,  is  20%  more  expensive  so  it  does  not  represent  a 
useful  alternative. 

(c)  Recycling  the  gas  mixture  is  technically  feasible  but  somewhat  con^licated  since  the 
gas  mixture  is  used  at  a  pressure  above  tiiat  for  condensation  of  Freon-12  at  room  temperature. 
Also,  a  relatively  large  amount  of  gas  (20~50  kg)  is  used.  Initial  estimates  are  that  such  a 
recycling  system  would  cost  between  $50,000  and  $100,000  and  take  4-6  months  to  get  working. 
Thus,  this  solution  is  not  compatible  with  the  current  project  budget  and  schedule. 

(d)  Running  tests  in  air  is  inexpensive  but  requires  a  25%  increase  in  rotor  speed  to 
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maintain  constant  corrected  speed.  This  represents  a  50%  increase  in  rotor  (no  problem)  and 
instrumentation  (may  be  a  problem)  stress  levels.  This  speed  increase  is  in  addition  to  the  20% 
increase  over  the  original  tunnel  design  speed  needed  to  operate  at  the  -10®  rotor  incidence  angle 
that  is  now  the  standard  operating  condition.  This  places  the  rotor  speed  at  8,000+  RPM,  very  near 
the  dynamic  limit  of  the  current  rotating  system.  Also,  of  course,  y  is  1.4  rather  than  1.28. 

(e)  Gas  mixtures  other  than  argon-Freon  have  been  considered.  One  of  the  most  benign  is 
argon  and  CXD2  mixed  to  the  proper  y  of  1 .28.  '•  is  has  a  molecular  weight  of  42  and  a  speed  of 
sound  1 14%  of  that  of  argon  and  Freon.  Thus,  it  represents  a  low  cost  alternative  with  a  less 
extreme  increase  in  rotor  speed  and  stress  levels. 

None  of  these  options  are  all  that  attractive  but  we  are  unaware  of  any  other  alternatives. 

The  current  plan  for  the  NGV-only  testing  is  to  (a)  run  most  tests  in  air,  (b)  run  a  set  with  argon- 
CO2  for  comparison,  and  (c)  run  two  Freon  tests  for  reference  and  comparison  with  previous 
running.  For  the  rotor  tests,  the  first  few  tests  will  be  with  the  standard  argon-Freon  mixture.  The 
principal  concern  here  is  to  avoid  stressing  of  and  damage  to  the  rotor  instrumentation  until  a 
baseline  set  of  tests  have  been  completed  (or  the  budget  is  exhausted).  Testing  will  then  switch  to 
argon-CC)2. 

12.0  Near-Term  Plans 

The  plan  for  this  winter  is  to  conduct  the  NGV-only  tests  at  the  end  of  December  and  in 
January.  The  rotor  instrumentation  will  be  calibrated  in  January  and  the  rotor  installed.  The  full- 
stage  testing  should  begin  before  February  and  continue  through  the  spring. 
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Fig.  3 :  Blowdown  turbine  facility  with  distortion  generator  and  upstream  rake  translator  installed 
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Fig.  8 :  The  exit  rake  translator  installation  in  the  lest  section 
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Fig.  9 ;  Alternate  locations  considered  for  turbulence  grid  placement 
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Hgure  13:  Hot  wire  calibration  graphs 


Figure  14:  Turbulence  grid  dimensions 
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Figure  15;  Reynolds  number  effect 
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Figure  16:  PSD;  1/2  NGV  location 
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Figure  17:  PSD;  NGV  location 
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Figxire  18:  PSD  nozzle  case  NGV  location 


PSD  (dB) 


Hgure  19:  Positioning  effect  parallel  to  nozzle 


PSD  (dB) 


Position  2 


Wavenumber  k(l/m) 


F.'^iac  20:  Positioning  effect  in  direction  of  nozzle  contraction 
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Figure  21 :  Turbulence  intensity  variation 


Conduction 
Reference  Junction 


57 


Fig.  22 ;  The  general  arrangement  of  the  fixed  MIT  fast  reponse  total  temperature  rake.  The  new  units  are  similar 
the  insulators  are  3.8  mm  nuher  than  7.6  mm  long. 
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Fig.  24 ;  Time  hist<^  of  the  ex 
Test  time  is  from  250 


Fig.  25 :  Inlet  temperature  rake.  Each  of  the  1 1  sensors  is  a  vented  design  similar  to  Fig.  22 . 
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Fig.  27:  Rake  calibration  tunnel 
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28 :  Temperature  sensor  distribution  in  turbine  supply  tank 
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Fig.  29 :  Influence  of  fans  on  supply  tank  temperature  distribution 


Fig.  30 :  Composite  of  rotor  blade  heat  flux  instrumentation  projected  onto  a  plane. 
Each  chordwise  row  is  on  a  separate  blade. 
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Fig.  32 :  3  NGV,  5  rotor  blade  grid  geometry  illustrated  for  the  hub  plane 


Fig.  33 :  Calculational  domain  with  total  temperature  contours  illustrated  on 
blade  surfaces.  There  arc  56  x  36  x  21  nodes  in  each  NGV  passage 
and  56  X  22  X  21  in  each  rotor  passage. 
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Fig.  34:  The  circumferential  period,  X,  of  the  inlet  temperature  distortion 

must  match  the  3  NGV-5  rotor  grid  geometry 


Rake  Head 


i^ig.  35a;  Turbine  exit  traversing  rake.  The  thermocouple  reference  junctions 
are  mounted  in  the  rake  body. 


Fig.  35c:  Total  temperature,  high  frequency  total  pressure,  and  low  frequency  total  pres¬ 
sure  rakes  mounted  at  NGV  exit  station  on  rotating  drum  (looking  downstream). 
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SCALED  TT5_01T.  TT5_05T,  AND  TT5_06T 
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Fig.  3  9  :Time  histories  of  three  exit  rake  thermocouples  during  traverse  in  vacuum. 
The  three  outputs  closely  line  up. 
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Fig.41 :  Inlet  and  outlet  T/C’s  and  RTD  plotted  with  angular  velocity  squared  (Vg) 
during  a  vacuum  traverse. 
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Fig.  45 :  Test  section  arrangement  for  NGV-only  testing 
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Fig.  46 :  Time  history  of  the  six  transducers  from  stationary  inlet  total  pressure  rake  during  a  blowdown  test. 
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Fig.  48 :  The  periodicity  of  the  NGV  exit  flow  is  evident  as  the  exit  rake  is  traversed  l)ehind  the  NGV’s. 
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Fig.  49 :  The  high  frequency  response  exit  total  pressure  rake  shows  the  NGV  exit  flow  field  during  a  test  traverse. 
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Fig.  50 :  A  circumferential  traverse  of  the  inlet  rjike  midspmi  station.  The  downspikes  at  250, 390,  and  550 
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Fig.  52:  Acont 
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Looking  Downstream 


Fig.  54; 


Rotating  multilayer  ribbon  arrangement  on  exit  translator. 
Shown  in  staning  position,  travel  is  350°. 
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Fig.  56:  Upstream  end  of  the  0TDF  sector  of  the  distortion  generator  showing  matrix,  heater  (thick  wire),  and 
thermocouples  (numerous  smaller  wires). 
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Fig.  57 :  0TDF  sector  response  to  step  heat  inputs. 


Rg.  53 ;  Rt  of  transient  conduction  parameters  (lines)  to  data  (symbols)  during  vacuum  heatup  of  distortion  generator. 
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Fig.  59 :  Model  predictions  of  eTDF  generator  performance  showing  two  hot  spots.  1  Icight  is  temperature. 


TASK  II:  ASSESSMENT  OF  UNSTEADY  LOSSES  IN  STATOR/ROTOR 
INTERACTIONS 

(Investigators:  M.B.  Giles,  G.  Friisch) 

Gerd  Fritsch,  the  graduate  student  who  was  working  on  this  project,  completed  his  Ph  D. 
thesis  in  May  1992.  The  thesis  was  subsequently  prepared  as  an  MIT  Gas  Turbine  Labcsatory 
Report  (GTL  Repeat  No.  210,  “An  Analytical  and  Numerical  Study  of  the  Sc<x>nd*Ordcr  Effects 
of  Unsteadiness  on  the  Performance  of  Turbomachines”)  and  is  available  on  request.  In  June,  a 
paper  was  presented  at  the  ASME  Gas  Turbine  Conference  in  Cologne.  This  focussed  on  the  loss 
mechanism  due  to  the  unsteady  Stokes  sublayer,  as  discussed  in  previous  annual  reports.  A  copy 
of  the  paper  is  given  as  an  Appendix.  A  second  paper  based  on  Fritsch's  wwk  is  now  in 
preparation  for  the  1993  ASME  Gas  Turbine  Conference  in  Cincinnati.  This  addresses  the 
analysis  of  time-averaged  and  space-averaged  mixing  loss  which  was  discussed  in  last  year's 
annual  report. 

Overall,  Fritsch's  research  has  made  an  important  ccwitribution  towards  the  growing 
understanding  of  the  impact  of  unsteadiness  in  turbomachinery  flow.  He  has  shown  that  the  loss 
due  to  unsteady  circulation  is  very  small  and  can  be  disregarded  in  fuUire  research.  He  has 
explained  and  analyzed  the  m^hanism  of  the  unsteady  Stokes  layer  loss.  Unfortunately,  he  has 
not  been  able  to  accurately  assess  its  magnitude  under  typical  engine  conditions  but  even  in  this  he 
has  made  the  important  contribution  of  showing  the  limitations  of  existing  CFD  methods  in 
addressing  this  important  question.  Finally,  his  analytical  work  on  the  mixing  loss  ties  together 
and  expands  upon  the  work  done  by  others  and  shows  the  separation  of  the  loss  due  to  different 
forms  of  unsteady  wave.  It  also  highlights  the  loss  that  is  introduced  in  steady-state  multistage 
CFD  analyses  using  mixing  plane  coupling. 

Also  presented  in  last  year's  report  woe  ideas  for  continuing  this  work  in  the  future  using  a 
combined  asynq)totic  analysis  of  unsteady  multistage  flow  fields.  In  the  analysis,  tiie  base  level  is 
a  steady-state  multi-row  calculatitHi  using  mixing  plane  coupling  at  all  statorAotor  interfaces.  The 
next  part  is  a  calculation  of  the  unsteady  flow  field  using  a  linear  perturbation  assumption.  The 


third  stage  calculates  the  change  in  the  time-averaged  flow  Held  due  to  quadratic  terms  asscx  latcd 
with  the  linear  unsteady  flow  field.  A  paper  was  presented  on  this  idea  at  the  ASME  Gas  Turbine 
Conference  in  Cologne  (also  included  in  the  Appendix).  A  number  of  related  issues  were  also 
discussed  at  a  Workshop  on  Unsteady  Losses  held  at  NASA  Lewis  in  May.  This  brought 
together  researchers  in  NASA,  academia  and  industry  to  discuss  the  key  aspects  of  unsteady  flow 
in  turbomachinery,  to  review  our  current  understanding  and  the  issues  that  need  to  be  resolved. 
Discussion  of  this  is  given  in  the  section  on  Interactions. 

Work  on  the  multi-stage  analysis  will  be  continued  by  Dr.  Giles  at  Oxford  University,  to 
which  he  moved  in  June. 
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Second-Order  Effects  of  Unsteadiness  on  the 
Performance  of  Turbomachines 

GERO  FHITSCH 
MICHAEL  GILES 

Oepanment  of  Aeronautics  and  Astronautics 
Massachusetts  institute  of  Technology 
Cam&fidge.  MA  02139 


Abstract 


A  2D,  lmeari2ed  approach  is  used  to  iovestigate  second-order 
effects  of  unsteadiness  on  the  time-mean  efficiency  of  turboma¬ 
chinery.  The  objective  is  to  quantify  unsteady  losses  with  a 
nonzero  time-mean  and  to  examine  numerical  simulations  with 
respect  to  the  modeling  of  unsteady  flow  fields  and  loss  mech¬ 
anisms.  Results  of  simulations  constitute  the  input  to  the  ana¬ 
lytical  models  employed. 

Two  unsteady  I'ou  mechanisms,  oneof  invisdd  and  the  other 
of  viscoua  nature,  are  conndered.  The  uiuteady  dxculation 
losses,  i.e.,  the  transfer  of  kinetic  energy  into  the  unsteady  part 
of  the  flow  field  through  vortidty  shed  at  the  trailing  edge  of  a 
blade,  was  first  considered  by  Keller  (1935)  and  later  by  Kemp 
and  Sears  (1956).  The  vortidty  is  shed  in  response  to  an  un¬ 
steady  blade- circulation  and  determined  from  Kelvin’s  dreuia- 
tion  theorem  which  is  valid  in  compressible  bomentropic  flow. 
Use  of  a  numerical  simulation  to  obtain  drculation  amplitudes 
avoids  the  limitations  of  thin>airfoil  theory  and  yields  results 
mote  realistic  for  modern  turbomachinery.  For  the  unsteady  vis¬ 
cous  lost  mechanism,  i.e.,  the  dissipation  in  an  unsteady  bound¬ 
ary  layer  on  the  blade  surface,  Ligh thill’s  high- reduced- frequeucy 
limit  (1955)  is  used  to  obtain  the  local  velodty  disinbution  in 
the  laminar  sublayer  and  the  corresponding  time-mean  unsteady 
dissipation-  The  input,  to  the  model  is  the  time-haimonics  of  the 
pressure  gradient  along  a  blade  surface  obtained- from  a  sinuila- 
tion.  A  numerical^ study  of  theertors  introduced  by- a. departure 
from,  the* hi^nduced-freqnency  limit  is  presented.  Losses  from 
both  sources  are  found  to  be  smalL 


Nomenclature 


Latin  Symbols 
a 
c 

h, 

k,  k,  and  hn 


speed  of  sound 

blade-  axial  chord 

specific  heat  at  constant  pressure 

total  specific  enthalpy  of  a  fluid 

wave  niunber,  reduced  frequencies 


n 

average  kinetic  energy  in  a  convetled  frame 

n 

coordinate  normal  to  a  blade  surface 

P 

pressure 

s 

entropy,  surface  coordinate,  surface  length 

t 

time 

tt/v 

velocity  in  the  direction  *,  i,  or  s  /  y,  y  ot  n 

X.  * 

axial/ftreamwise  coordinate,  convected  z-coord. 

y.  y 

coordinate  normal  to  z,  x 

L.L 

incident  and  reflected  acoustic  intensity 

1^,W. 

mean  Mach  numbers  (U/'S),  (U.I'S.) 

P,  P 

pitch,  pitch  in  a  convected  frame  (/  =  Pcosu] 

«() 

real  part  of  a  quantity 

Re 

Reynolds  number  based  on  chord. 

Re.,  Re, 

,  Re  based  on  surface  length,  momentum  thickess 

A? 

average  nondimeusionai  entropy  rise 

v,v. 

mean  convection  velocity,  velocity  at  ‘.he  bJ.  edge 

Vn 

rotor  speed 

Greek  Symbols 

a 

angle  between  mean  flow  and  axial  direction 

inter-blade  phase  angle,  attenuation  factor 

4 

inter-wake  phase  angle  ia  convected  frame 

7 

ratio  of  ipediic  heat,  vortex  street  strength 

Is, 

steady  and  unsteady  boundary  layer  thickness 

s, 

laminar,  sufalayex  thickness  in  turbulent  flow 

li,  V,  K 

viscosity,  kinematic  visconiy,  conductivity 

total  pressure-  ratio 

« 

smoothing  coefficient 

V. 

uentiopic  efficiency 

P 

density 

<r,  r 

norma]  stress,  shear  stress 

$ 

perturbation  potential 

e 

angle  between  wave  vector  ard  surface  normal 

a  reduced  frequency  A„  =  jkm4i  wavelength 

W,  W/ 

angular  frequency,  blade  passing  frequency 

r 

circulation 

A 

difference,  grid  spacing 

* 

dissipation  function 

PrstantM-st  tti*  Intsmallanal  Q*t  Tufbine  «n<)  Asnicnqin*  Congran  snil  Emosiuon 
Cotogn*.  G«rmviy  4ufi«:l-4.  19^ 
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Subscripts 

e 

at  the  boundary  layer  edge 

t 

integrated 

in 

steady  uniform  fluid  state  far-upstream 

mixed 

mixed-out  state  far-downstream  of  the  rotor 

out 

unsteady  nonuniform  state  at  the  rotor  outflow 

s 

isentropic 

w 

at  the  wall 

D,  T 

design,  turbine 

HFL 

high-frequency  limit 

1,2 

first/second-order  (perturbation) 

Superscripts 

unsteady  quantity 

A 

complex  amplitude 

- 

time-mean,  time-mean  and  space-mean 

Introduction 

Flow  fields  in  turbomachinery  are  inherently  unsteady.  In  a 
turbine,  the  first  blade  row  will  be  subject  to  a  nonuniform  dis¬ 
tribution  of  velocity  and  total  temperature  due  to  incomplete 
mixing  in  the  combustion  chamber.  The  relative  motion  of  two 
adjacent  blade  rows,  in  conjunction  with  the  spatially  nonuni¬ 
form  pressure  field  locked  to  a  loaded  blade,  leads  to  an  unsteady 
pressure  distribution  through  potential  interaction.  A  blade  row 
may  move  through  and  interact  with  the  shock  system  originat¬ 
ing  at  the  trailing  edges  of  the  upstream  blade  row.  Stator 
wake/rotor  interaction  causes  unsteadiness  in  the  rotor  frame  of 
reference  because  the  rotor  is  moving  through  and  chopping  the 
wakes.  Similarly,  secondary  flow  effects  like  horse  shoe  vortices, 
passage  vortices,  and  tip  clearance  vortices  contribute  to  flow 
unsteadiness.  The  viscous  flow  past  a  blunt  turbine  trailing 
edge  results  in  vortex  shedding;  trailing  edge  vortex  shedding 
has  also  been  found  in  compressors  operating  in  the  transonic 
or  supersonic  regime.  Finally,  there  is  unsteadiness  induced  by 
the  motion  of  the  blades  themselves,  i.e.,  blade  flutter. 

Nevertheless,  the  standard  aerodynamic  design  tools  for  tur¬ 
bomachinery  are  steady  state  codes  both  inviscid  and  viscous, 
ud  steady  cascade  experiments.  Designing  a  single  stage  or 
blade  row  with  steady  state  tools  amounts  to  placing  the  stator 
and  rotor  infinitely  far  apart  thus  eliminating  the  effects  of  blade 
row  interaction  and  the  interaction  of  convected  disturbances 
with  the  downstream  blade  rows.  Unsteadiness,  however,  con¬ 
tributes  additional  losses  with  nonsero  time-mean.  First,  the 
energy  associated  with  the  niuteady  part  of  the  flow  field  is 
usually  not  recovered  in  the  blade  pass^e  or  in  downstream 
blade  rows;  it  will  eventually  be  lost  to  viscous  dissipation.  Sec¬ 
ond,  the  interaction  of  the  unsteadiness  with  boundary  layers 
and  shock  structures  can  trigger  additional  dissipation  (Hodson, 
1983;  Ng  and  Epstein,  1985).  In  a  steady  viscous  simulation,  the 
effects  of  unsteadiness  and  unsteady  viscous  dissipation  on  the 
eflSciency  are  not  captured.  An  unsteady  nonlinear  simulation 
is  still  prohibitively  expensive  for  routine  design  purposes  and 
will  likely  remain  so  for  the  foreseeable  future,  in  particular  for 
multi-stage  turbomachinery  Furthermore,  the  weakest  point  of 
any  numerical  simulation,  steady  or  unsteady,  is  the  accurate 
prediction  of  losses  and  heat  loads  due  to  the  unavailability  of 
adequate  turbulence  and  transition  models.  Testing  of  a  stage  or 


a  whole  turbine  under  unsteady  operating  condUior.s,  reproduc 
mg  the  whole  range  of  unsteady  effects,  will  remam  irnpracticaJ 
for  routine  design  purposes. 

One  of  the  niajn  impulses  for  this  investigation  was  the  devel¬ 
opment  of  the  CFD-code  UNSFLO  by  Giles  (1987,  1990,  19yla) 
and  the  visualisation  package  V1SUAL2  by  Giles  and  Baimes 
{199ib) ,  which  provides  a  detailed  picture  of  unsteady  flow  phe 
nomena.  Those  include  wakc/blade  row  interaction,  unsteady 
shock  wave  structures,  unsteady  boundary  layers,  unsteady  heat 
transfer,  and  unsteady  vortex  shedding.  There  is  a  need  to  un¬ 
derstand  the  physical  nature  of  these  phenomena,  to  quantify 
their  effect  on  the  time-mean  efficiency  of  turbomachinery,  and 
to  examine  the  modeling  accuracy  of  CFD-codes  with  respect  to 
these  phenomena 

An  alternative  to  nonlinear  unsteady  simulations,  which  is 
receiving  increased  attention,  are  linesir  perturbation  methods 
A  nonlinear  steady  state  is  found  and  the  unsteady  poitioi;  ai c 
flow  field  superimposed  upon  it  as  a  small  perturbation.  Second- 
order  terms,  i.e.,  terms  quadratic  in  unsteady  quantities,  are  ne¬ 
glected  since  the  perturbations  are  assumed  to  be  small.  Linear 
perturbation  codes  have  been  found  to  give  accurate  results  op 
to  a  surprisingly  high  level  of  unsteadiness  (Hall  and  Crawley, 
1987;  Lindquist  and  Giles,  1991)  and  will  be  more  widely  used 
in  the  future.  Linear  codes,  bke  steady  codes,  cannot  capture 
unsteady  losses.  The  time-mean  of  the  first-order  unsteady  dis¬ 
sipation  is  zero;  only  terms  second-order  in  unsteadiness  have  a 
nonzero  time-mean  dissipation. 

The  turbomachinery  community  it  moving  towards  the  con¬ 
sensus  that  increased  losses  under  unsteady  operating  conditions 
are  primarily  due  to  nonlinear  effects  like  the  alteration  of  the 
boundary  layer  characteristics  through  their  effect  on  irassitioa 
(Hodson,  1983),  the  variation  of  secondary  flow  generation  in 
downstream  blade  passages  and  their  effect  on  separation  or 
reattachment.  Nevertheless,  the  magnitude  of  effects  that  can 
be  treated  in  a  linear  framework,  remains  undetermined.  While 
the  corresponding  unsteady  losses  are  small  by  definition,  their 
exact  magnitude  is  of  importance  to  the  design  process.  Small 
changes  in  the  aerodynamic  performance  figure  significantly  over 
the  lifetime  of  an  engine.  In  a  turbine,  for  example,  eflBaencies 
often  exceed  90%  and  improvements  will  be  incremental. 

This  paper  treats  two  aspects  in  a  2D,  linearized  framework. 
In  the  first  part,  it  revisits  the  Unsteady  Circulation  Loises  first 
examined  by  Keller  (1935)  and  later  by  Kemp  and  Sears  (1956). 
Keller  estimated  the  unsteady  dxcnlation  amplitude  to  arrive  at 
the  average  kinetic  energy  in  the  unsteady  part  of  the  flow  field 
induced  by  the  vorticity  shed  at  the  trailing  edge.  Sears  and 
Kemp  used  thm-airfoil  theory  to  obtain  the  unsteady  circula¬ 
tion  amplitude.  This  approach,  while  enabling  them  to  calcu¬ 
late  the  drcnlation  amplitudes,  limited  them  to  incompressible 
flow,  blades  of  zero  thickness  and  camber  with  the  mean  flow 
nearly  in  the  blade  direction,  i.e.,  lightly  loaded  blades.  Thus, 
the  airfoils  are  more  representative  of  compressor  blades  than 
turbine  blades.  Furthermore,  they  bad  to  specify  the  steady  Hft 
distribution  on  the  blade.  In  this  paper,  the  fluctuation  am¬ 
plitudes  are  the  result  of  CFD-simulations,  which  allows  one  to 
obtain  fluctuation  ampUtudes  for  arbitrary  blade  and  stage  ge¬ 
ometries,  steady  lift  distributions,  and  Mach  numbers.  Kelvin’s 
circulation  theorem,  upon  which  Keller’s  work  rests,  is  valid 
even  in  compressible  homentropic  flows.  Eliminating  the  need 
to  estimate  the  fluctuation  amplitudes  or  to  deduce  them  from 
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thin-airfoil  theory,  results  in  a  mote  realistic  measure  for  the  sec¬ 
ondary  kinetic  energy  in  modern  turbomachines  and  the  losses 
associated  with  its  dissipation.  The  dissipated  energy  is  related 
to  the  efficiency  through  a  control  volume  argument  and  the 
linearization  of  the  unsteady  isentropic  efficiency.  Exsimples  are 
given  for  several  blade  rows  and  stages. 

The  second  part  considers  Unsteady  Viscous  Losses,  that 
is,  dissipation  in  unsteady  boundary  layers  on  the  blade  sur¬ 
faces.  Using  linearization,  Fourier-decomposition  in  time,  and 
Lighthill’s  high-reduced-frequency  limit,  the  momentum  equa¬ 
tion,  driven  by  the  harmonics  of  the  pressure  gradient  obtained 
from  a  simulation,  yields  the  local  velocity  distribution  in  the 
laminar  sublayer.  While  a  numerical  simulation  has  difficulties 
in  accurately  predicting  the  losses,  the  pressure  distribution  is 
predicted  accurately  in  attached  flows.  The  second-order  dissi¬ 
pation  function  with  nonzero  time-mean  is  calculated  from  the 
induced  velocity  distribution  and  the  corresponding  entropy  in¬ 
crease  integrated  over  the  blade  surface.  The  analysis  can  be 
regarded  as  a  viscous  complement  to  Johnson  et  al.  '1988), 
who  studied  the  surface  heat  transfer  fluctuations  in  response  to 
shock  wave  passing.  The  results  of  a  numerical  study  to  check 
the  errors  introduced  by  a  departure  from  the  high-frequency 
limit  are  discussed.  The  model  is  applied  to  a  transonic  tur¬ 
bine  stage.  Comparing  the  entropy  rise  due  to  unsteady  vis¬ 
cous  losses  to  the  total  entropy  rise  in  the  simulation,  it  is  seen 
that  unsteady  viscous  losses  account  only  for  a  small  part  of 
the  increased  losses  under  unsteady  operating  conditions.  A 
simple  analytical  model  for  the  propagation  of  pressure  waves 
is  a  CFD'SiffiiiJation  is  used  to  show  that  the  principal  loci  of 
unsteady  dissipation  in  a  simulation  may  not  be  in  unsteady 
boundary  layers  as  in  real  flows.  Due  to  the  action  of  numerical 
smoothing,  dissipation  within  the  blade  passage  can  dominate. 


Unsteady  Circulation  Losses 

Theory 

Figure  1  shows  a  cascade  in  which  every  blade  is  experiencing  an 
unsteady  drculation  due  to  the  presence  of  and  interaction  with 
the  upstream  blade  row  in  an  elementary  stage.  The  derivation 
proceeds  along  the  lines  of  Keller  (1935).  The  coordinate  system 
has  been  chosen  such  that  the  mean  velocity  17  far  downstream 
is  in  the  x-direction;  the  vortex  sheet  shed  from  the  nth  blade 
is  located  at  y  =  nP.  P  is  the  blade  pitch,  P  the  rotated  pitch 
defined  as  (Fcota),  and  a  the  angle  between  the  direction  of 
the  convected  vortex  sheet  and  the  axial  direction. 

The  drculation  of  the  nth  blade  is  given  by 


t 


V, 


Figure  1:  A  cascade  with  unsteady  lift  and  circulation 


7„(z,()  convected  with  the  mean  flow  is  related  to  the  rate  of 
change  of  the  drculation  by 


r  —  nP  sin  a'^) 

V  )\ 


(2) 


With  equation  (1)  and  the  convected  coordinate  z  =  z~Vt,  this 
can  be  expressed  in  the  form 

in  the  convected  coordinate  system  (x,y).  The  reduced  fre¬ 
quency  kn  =  Wmc/l^  it  defined  with  the  blade  axial  chord  c; 
0„,  denotes  the  inter- wake  phase  angle  of  the  mth  mode.  In  the 
stationary  frame,  the  linearized  equation  for  the  perturbation 
potential  ^  describing  the  vortex  sheet  far  downstream  of  the 
trailing  edge,  is 


where  a  denotes  the  speed  of  sound.  It  does  not  describe  per¬ 
turbations  convected  through  the  blade  row.  The  unsteady  and 
convective  contributions  are  eliminated  by  a  shift  to  the  con- 
vected  frame  of  reference;  all  that  remains  is 


di^  5y* 


=  0. 


(5) 


In  the  convected  frame,  the  perturbation  potential  satisfies  the 
equation  for  steady  irzotational  flow.  Due  to  the  linear  nature 
of  equation  (5),  one  temporal  mode  can  be  treated  at  a  time. 
The  periodicity  condition  for  the  mth  mode  is 


r»  (<)  =  r  4-  53  «P  {-*  ( * ) 

m 

with  0„  being  the  inter-blade  phase  angle  for  the  mth  mode. 
The  frequenriei  ««  correspond  to  the  blade  passing  frequency 
w/  and  its  harmonics.  The  real  part  of  any  physical  quantity  is 
implied  by  complex  notation.  The  strength  of  the  unsteady  vor¬ 
tex  sheet  shed  at  the  trailing  edge  is  determined  by  the  condition 
that  the  sum  of  bound  and  shed  vorticity  remain  constant.  This 
is  a  consequence  of  Kelvin’s  circulation  theorem,  which  is  valid 
for  any  compressible,  inviscid,  and  barotropic  fluid  with  conser¬ 
vative  body  forces  sicting  on  it.  The  strength  of  a  vortex  sheet 


+ P)  =$mii,Ti)exp{0n,}  ■  (6) 

The  potential  for  0  <  y  <  .^,  is  of  the  form 

cosh  + 

J&,„smh^hm^-A„^jexp|j(h«^-(-^)|  (7) 

where  is  defined  by  A™  =  The  values  of  the  com¬ 

plex  amplitudes  A™  and  B„  are  determined  from  the  matching 
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conditions  at  1/  =  0.  First,  the  jump  in  the  streaxnwise  velocity 
perturbation  at  p  =  nP  must  equal  the  strength  of  the  nth 
vortex  sheet.  Second,  there  must  be  no  jump  in  the  normal  ve¬ 
locity  perturbation  u„  across  it.  Using  the  periodicity  condition 
(6)  to  relate  the  velocities  at  (i,0")  to  those  at  the 

perturbation  potential  becomes 


sin’  ^  +  sinh’  A„ 


i  cosh  A„  sin  cosh 


4  sinh  A„  cos  sinh  ~  j  ^  2  ‘ 


(8) 


The  subscript  ‘in’  refers  to  a  uniform  and  steady  slate  far  up¬ 
stream  of  the  stator  while  ‘out,  s’  and  ‘out’  refer  to  the  un¬ 
steady  spatially  nonuniform  rotor  outfiow  in  an  (hyp(>t>ie!:<:iih 
isentropic  expansion  and  in  the  actual  expansion,  respectively. 
Linearization  of  (12),  calculus  of  small  variables,  and  definition 
of  an  isentropic  ‘design  point  efficiency’  between  the  steady  and 
uniform  states  ‘in’  and  'mixed’, 


Vi.T.D  = 


(iJ) 


lead  to  the  result 


Efficiency  Drop 


The  average  kinetic  energy  per  unit  mass  associated  with  the 
convected  vortex  sheet  is 

’jl'k  /.*'  I  !<*"-»’  ^  <"''-''’1  ■'»  ““ 

_  1  A„  sinh  Xn  cosh  A„ 

^  -P’  sin’  ^  4  sinh’  A„ 

which  coincides  with  the  result  of  Keller  (1935).  3?  denotes  the 
real  part  of  a  complex  quantity.  Questions  relating  to  cross- 
induced  contributions  in  the  kinetic  energy  from  disturbances 
convected  through  the  blade  tow,  which  were  brought  up  by 
Kemp  and  Sears  (1956),  are  addressed  in  Fritsch  (1992).  They 
will  not  be  discussed  here  because  they  do  not  bear  heavily  on 
the  conclusions  drawn. 


It  is  assumed  that  the  kinetic  energy  associated  vrith  the 
unsteady  vortex  sheet  cannot  be  recovered  downstream.  If  this 
kinetic  energy  is  loti  to  viscous  dissipation  in  a  constant  area 
duct,  the  average  entropy  rise  per  unit  mass  for  an  ideal  gas 
with  constant  specific  heats  7  is 

The  subscript  ‘mixed’  in  equation  (10)  denotes  the  uniform  state 
after  mixing;  the  oveibar  implies  a  mass-average  taken  before 
the  mixing  process  that  dissipates  the  kinetic  energy.  Elquation 
(10)  holds  regardless  of  the  mean-flow  Mach  number  "M  and  in 
the  presence  of  unsteady  waves  other  than  vorticity  waves,  like 
pressure  waves  or  entropy  waves.  The  entropy  rise  associated 
with  dissipation  of  the  kinetic  energy  (9)  is 


ATf  ('T  f Stttfi  cosh  An, 

4  W  ri„*^  +  rinh’A„' 


(11) 


In  spatially  nonuniform  unsteady  flow,  it  is  the  mass-averages  of 
total  pressure  and  total  temperature  at  inflow  and  outflow  that 
are  needed  to  determine  the  isentropic  efiSdency.  As  mixing  oc¬ 
curs  downstream  of  the  rotor  (or  stator)  exit,  mote  losses  are 
generated  and  the  mixed-out  uniform  and  steady  state  leads  to 
a  lower  total  pressure.  If  the  energy  assodated  with  unsteadi¬ 
ness  is  not  recovered  then  it  is  the  appropriate  state  to  use  in 
the  isentropic  effidency.  The  isentropic  efficiency  for  unsteady 
nonuniform  flow  in  an  elementary  turbine  stage  with  adiabatic 
blades  is 


(M.)m  - d(j() 


A7.,r 

V>.T.D 


1  —  X 
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AS. 


(14) 


The  'design  total  pressure’  ratio  x,,n  is  defined  as  (pi,mw.^/pt,«i)- 
Equation  (14),  in  conjunction  with  (11),  tells  how  much  of  the 
departure  of  the  isentropic  efficiency  from  unity  can  be  ascribed 
to  the  dissipation  of  nnsteady  vorticity  downstream  of  the  rotor 
outflow.  An  analogous  result  can  be  derived  for  a  compressor 
stage.  Details  of  tbe  derivation  for  both  turbine  stages  and 
compressor  stages  are  found  in  Fritsch  (1992). 


Results  -  Unsteady  Circulation  Losses 

The  model  described  in  the  previous  section  was  applied  to  a 
number  of  blade  rows  and  stages.  The  amplitudes  of  the  un¬ 
steady  drculation  were  calculated  with  the  CFD-code  UNSFLO 
by  Giles  (Giles,  1987,  1990;  Giles  and  Baimes,  1991a).  It  uses  an 
explicit  Lax-Wendroff  scheme  to  solve  the  unsteady  compress¬ 
ible  Euler  equations  in  2D  on  an  unstructured  grid  composed 
of  triangular  and/or  quadrilateral  cells.  Its  quasi-3D  capability 
accounts  for  a  varying  stieamtnbe  thickness;  arbitrary  pitch  ra¬ 
tios  are  handled  by  a  time-inclined  computational  domain.  The 
simulations  were  run  until  the  flow  was  periodic;  then  50  instan¬ 
taneous  flow  solutions  equally  spaced  in  time  were  stored,  from 
which  the  circulation  ampUtudes  have  been  determined. 


Single  blade  row  simulations 

This  section  combines  the  results  of  simulations  for  three  differ¬ 
ent  rotor  blade  rows  with  incoming  wake  disturbances  modeled 
at  the  inflow  boundary.  The  wake  parameters  were  chosen  from 
the  experimental  data  for  the  respective  stages;  no  effort  was 
made  to  account  for  the  viscous  decay  of  the  wakes  in  the  exper¬ 
iment.  Potential  disturbances  were  also  simulated,  bnt  their  ef¬ 
fect  was  found  to  be  comparably  small  due  to  the  Mach  numbers 
and  axial  gaps  in  the  corresponding  experiments.  The  numer¬ 
ical  implementation  of  wake  models  and  potential  disturbance 
models  is  treated  in  Giles  ( 1988). 

Hodson  (1984b)  conducted  measurements  of  wake-generated 
unsteadiness  in  the  rotor  of  a  large-scale,  low-speed,  single-stage, 
air-driven  turbine.  The  51  rotor  blades  are  separated  from  36 
stator  blades  by  an  axial  gap  of  about  one  rotor  axial  chord.  A 
corresponding  inviscid  simulation  of  the  2D  wake/rotor  interac¬ 
tion  has  been  presented  by  the  same  author  (Hodson,  1984a); 
the  incoming  wake  was  modeled  after  the  experiment  (Hodson, 
1984b).  The  wake  model  in  the  simulation  for  this  case,  case 
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A  of  table  1,  was  identical  to  his  and  the  maximum  velocity 
defect  and  the  wake  width  were  chosen  such  as  to  replicate  the 
maximum  velocity  deficit  and  the  momentum  deficit  in  the  ex¬ 
periment. 

Case  B  is  a  simulation  of  the  flow  through  the  rotor  of  a 
large-scale,  axial-ilow  turbine  model  investigated  experimentally 
by  Joslyn  et  al.  (1983).  It  is  a  one  and  one-half  stage  turbine 
model  consisting  of  22  first  vanes,  28  first  blades,  and  28  sec¬ 
ond  vanes,  representing  the  first  three  blade  rows  of  the  high 
pressure  turbine  of  a  high-bypass-ratio  aircraft  turbine.  The 
Reynolds  number  as  well  as  the  absolute  and  relative  flow  an¬ 
gles  are  typical  for  current  high-pressure-turbine  airfoils.  The 
wake  profile  at  the  rotor  inflow  boundary  was  Gaussian  (Giles, 
1988).  Wake  width  and  maximum  velocity  defect  were  chosen  to 
match  the  maximum  velocity  defect  and  the  momentum  defect 
measured  at  midspan  downstream  of  the  first  vane. 

Binder  (1985)  used  the  L2F-technique  to  measure  the  flow 
in  a  cold-air  turbine.  The  stage  has  20  stator  blades,  31  rotor 
blades,  and  an  axial  gap  of  113%  rotor  axial  chord;  it  is  a  scaled 
version  of  the  first  stage  of  a  helicopter  gas  turbine.  The  flow 
in  the  low-aspect  ratio,  high  hub-to-tip  ratio  stage  is  weakly 
transonic  and  the  experiments  showed  a  wake  with  a  velocity 
defect  of  6.5%  at  midspan  50%  rotor  axial  chord  upstream  of  the 
rotor  leading  edge.  Based  on  the  measured  velocity  distribution 
at  the  rotor  inlet  plane,  a  sinusoidal  '.rake  model  (Giles,  1988) 
with  a  maodmum  velocity  defect  of  6.5%i  in  the  stator  frame  was 
chosen  for  the  simulation  of  the  flow  through  the  rotor  (case  C). 


Parameter 

Case  A 

Case  B 

Case  C 

inflow  angle 

2.1" 

48.9" 

46.7" 

outflow  angle  a 

-66.8* 

-65.5" 

-62.0* 

inflow  Mach  number 

0.14 

0.12 

0.43 

outflow  Mach  number 

0.40 

0.20 

0.65 

rotor  speed  Vr 

0.31 

0.11 

0.58 

reduced  frequency  ki 

3.45 

2.92 

4.25 

wake  velocity  defect 

4.6% 

15% 

6.5% 

Wake  width/stalor  pitch 

30% 

10% 

— 

circulation  Ti/Fo 

0.67% 

0.57% 

2.83% 

entropy  rise 

1.7*10-* 

1.2*10-* 

3.6*10-* 

efiSciency  drop  ATj,jrht,Tj> 

-0.045^0 

-0.11% 

-0.24% 

Table  1;  Results  for  single  blade  row  simulations 


The  results  are  summarized  in  table  1.  Mach  numbers  and 
flow  angles  are  mixed-out  averages  in  the  rotor  frame.  Un¬ 
listed  higher  harmonics  of  the  drculation  are  a  factor  10  or  more 
smaUer  than  the  first  harmonic. 

Single  stage  simulations 

Two  invisdd  simulations  are  presented  here;  one  for  a  turbine 
stage  and  one  for  a  compressor  stage.  In  contrast  to  the  previous 
three  cases,  the  unsteadiness  acting  on  the  downstream  blade 
row  is  not  prescribed  but  the  result  of  blade  row  interaction  in 
the  simulation.  Viscous  wakes  and  their  subsequent  interaction 
with  the  downstream  blade  row,  some  aspects  of  which  could  be 
modeled  inviscidly,  arc  not  included. 

The  turbine  case  considered  is  a  highly  loaded  transonic  high- 
pressure  turbine  stage,  often  referred  to  as  the  ACE  turbine 
stage,  which  was  tested  at  MIT  (Guenette  et  al.,  1988)  .  A  2D- 


version  of  the  rotor  profile  was  investigate'd  at  Oxford  University 
(Ashw'orth  et  al.,  1985;  Johnson  et  al..  19SS)  .  In  the  simula¬ 
tion,  stator  and  rotor  tow  are  separated  by  an  axial  gap  of  40‘’l 
rotor  axial  chord  and  the  rotor  operates  near  the  aesign  inflow 
angle  of  58".  The  outflow  of  rotor  and  stator  row  are  supersonic, 
producing  weak  oblique  shocks  at  the  trailing  edges.  A  detailed 
computational  account  of  the  inviscid  stator/rotor  interaction 
for  almost  identical  input  parameters  was  presented  by  Giles 
(1990).  The  oblique  shock  wave  originating  at  the  norile  guide 
vane  traibng  edge  first  hits  the  downstream  rotor  around  the 
crown  of  the  suction  surface  and  is  reflected  from  there.  Subse¬ 
quently,  it  moves  forward  to  the  leading  edge  with  the  reflected 
portion  impinging  on  the  pressure  surface  of  the  adjacent  blade 
and  being  reflected  to  the  suction  surface  once  mote.  Upon  lift¬ 
ing  off  the  rotor  blade,  the  primary  shock  wave  grows  in  length 
until  the  relative  motion  of  the  blade  tows  causes  it  to  strike  the 
next  rotor  blade.  In  the  mean  time  its  reflection  off  the  rotor 
leading  edge  has  impinged  on  and  reflected  from  the  nozzle  guide 
vane  suction  side.  Comparatively  large  unsteady  ampbtudes  for 
this  strong  interaction  are  expected  and  have  been  reported  for 
this  turbine  stage  (Giles,  1990).  The  increased  importance  of 
higher  harmonics  is  a  consequence  of  the  discontinuous  nature 
of  the  shock  waves. 


Parameter 

ACE 

Stage  67 

stator  inflow  angle 

0" 

-44.7" 

stator  outflow  angle 

73.5* 

3.5" 

rotor  inflow  angle 

56.7* 

59.9" 

rotor  outflow  Angle 

-67.6* 

39.7* 

stator  inflow  Mach  number 

0.14 

0.67 

stator  outflow  Mach  number 

1.13 

0.46 

rotor  inflow  Mach  number 

0.59 

l.ll 

rotor  outflow  Mach  number 

1.16 

0.62 

rotor  speed  Vr 

0.55 

0.83 

reduced  Frequency  k\ 

2.38  (rotor) 

21.5  (stator) 

circulation  I\ITo 

2.19% 

0.90% 

Fr/Fo 

ro2% 

Oo% 

1  r,/r„  »  0.20% 

0.49% 

Fs/Fo 

0.25% 

0.14% 

entropy  rise  Ym  A5m 

1.1  *  10-^ 

4.0  *  I0-‘ 

efficiency  drop 

-0.30% 

-0.03% 

Table  2:  Results  for  single  stage  simulations 


The  compressor  stage  is  known  as  NASA  stage  67.  It  is  a 
low-aspect  ratio,  single-stage,  transonic,  axial-flow  compressor 
designed  and  experimentally  investigated  at  the  NASA  Lewis 
Research  Center  (Hathaway  et  al.,  1987).  The  stage  consists 
of  22  rotor  blades  separated  from  34  stator  blades  by  an  ax¬ 
ial  gap  of  78%  rotor  axial  chord.  At  a  rotor  tip  Mach  number 
of  1.38,  strong  shocks  are  present  in  the  rotor  passage.  Shan , 
(1989)  used  UNSFLO  to  c^cnlate  the  flow  through  the  stare 
at  midspan  with  the  input  parameters  inferred  £com  the  experi¬ 
ments  at  NASA  Lewis.  The  case  considered  here  is  a  full  stage 
with  purely  potential  interaction  between  the  blade  rows.  The 
magnitude  of  the  potential  interaction  strongly  depends  on  the 
size  of  the  axial  gap  separating  the  blade  rows.  A  wide  range  of 
gap  sizes  exists  in  compressors  and  the  gap  itself  varies  hub  to 
tip.  For  the  case  presented  here,  the  gap  was  10%  of  the  rotor 
axial  chord.  The  shock  position  is  roughly  30%  rotor  axial  chord 
downstream  of  the  rotor  leading  edge  and  there  is  no  significant 
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interaction  between  the  shock  system  and  the  downstream  stator 
blade  row.  As  a  consequence,  the  unsteady  circulation  and  the 
losses  are  significantly  smaller  than  for  the  ACE  turbine  stage. 

Summary 

The  magnitude  of  the  unsteady  circulation  amplitude,  estimated 
at  5%  of  the  mean  by  Keller  (1935)  ,  tsrpically  was  of  the  order 
of  1%  with  larger  values  reserved  for  highly  loaded  blade  rows 
with  high  turning.  In  their  incompressible  considerations  within 
the  limits  of  thin-airfoil  theory,  Kemp  and  Sears  (1956)  had  con¬ 
cluded  that  the  amplitude  fluctuations  are,  generally,  substan¬ 
tially  smaller  than  5%.  The  values  found  here  suggest  that  5% 
can  be  considered  an  upper  bound  for  the  unsteady  circulation 
amplitude. 

The  amount  of  kinetic  energy  associated  with  the  shed  vor- 
tidty,  which  is  proportional  to  the  square  of  the  unsteady  circu¬ 
lation  amplitude,  is  small.  Keller  (1935)  deduced  a  magnitude 
of  0.4  to  1%  of  the  power  consumed  or  delivered  by  the  rotor 
for  circulation  amplitudes  of  5%.  The  values  reported  in  this 
paper  are  smaller  with  a  maximum  efficiency  drop  of  0.3%  for 
the  ACE  turbine  stage  which  shows  a  strong  shock  wave/  blade 
row  interaction.  The  effect  of  incoming  wake  disturbances  can 
be  of  the  same  order  of  magnitude. 

Potential  interaction  was  found  to  be  of  minor  significance 
for  the  cases  considered. 


( 


Unsteady  Viscous  Losses 
Conceptual  Model 


Unsteady  viscous  losses  are  losses  with  a  nonzero  time-mean 
due  to  viscous  dissipation  of  energy  associated  with  or  induced 
by  the  unsteady  part  of  the  flow  field,  primarily  on  the  blade 
surfaces. 

The  unsteady  velocity  distribution  induced  by  a  wall  oscil¬ 
lating  harmonically  with  firequency  w/  in  a  fluid  at  rest,  Stokes’ 
second  problem,  is  shown  in  figure  2.  On  the  left,  the  veloc¬ 
ity  distribution  seen  in  a  stationary  frame  of  reference  has  been 
sketched  for  two  particular  instants  during  a  period.  The  length 
scale  characterizing  the  extent  to  which  unsteady  viscous  effects 
spread  away  from  the  wall  during  a  period  is  fiy  =  ^vjwf.  On 
the  right,  the  unsteady  velocity  distribution  is  shown  in  a  frame 
of  reference  fixed  to  the  wall.  The  same  velocity  distribution 
at  on  the  right  it  obtained  with  a  harmonically  varying  static 
pressure  gradient,  if  the  characteristic  length  of  the  pressure 
disturbance  is  much  larger  than  th«:  thickness  of  the  unsteady 
boundary  layer  it  indnces. 

On  a  turbomachinery  blade,  figure  3,  both  a  steady  ud  an 
unsteady  boundary  layer  are  present  together  with  a  mean  flow 
pressure  gradient.  The  prMsure  disturbance  can  be  a  pressure 
wave,  the  potential  field  of  another  blade  row,  or  a  weak  shock 
wave.  Far  from  the  wall,  the  unsteady  velocity  distribution  in 
the  boundary  layer,  driven  by  the  harmonically  varying  pressure 
gradient,  differs  from  the  one  in  Stokes’  second  problem.  It  is 
influenced  by  a  nonzero,  transitional  or  turbulent,  mean  flow 
and  compressibility  effects.  In  the  immediate  neighborhood  of 
the  wall,  where  the  unsteady  dissipation  is  concentrated,  the 


Figure  2:  Unsteady  boundary  layer  on  an  oscillating  wall 


mean  flow  unsteady  flow  superposition 

atui/l=0  at  iii/t  =  0  at  w/t  =  y 


Figure  3:  Unsteady  boundary  layer  on  a  blade 


pressure  disturbance  still  induces  a  velocity  distribution  as  found 
in  Stokes’  second  problem  for  sufficiently  high  reduced  frequency. 


The  energy  associated  with  the  unsteady  flow  field  is  dissi¬ 
pated  either  in  the  free  stream  or  in  unsteady  boundary  layers  at 
the  blades,  the  hub  and  the  tip.  To  compare  the  relative  impor¬ 
tance  of  the  two  dissipation  loci  in  laminar  flow,  the  following 
order-of-magnitude  argument  is  instructive.  In  the  free  stream 
the  velocity  gradient  is  set  by  the  magnitude  of  the  unsteady 
velocity  u.,  and  the  characteristic  length  //,  of  the  disturbance. 


(15) 


For  a  pressure  wave  traveling  in  or  against  the  free  stream  di¬ 
rection,  the  characteristic  length  1/  becomes  the  wavdength  A/. 
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(16) 
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The  reduced  frequen(y  is  defined  in  terms  of  the  blade  passing 
frequency  wy,  the  local  mean  vdocity  in  the  free  stream  (at  the 
boundary  layer  edge)  F,,  and  the  axial  chord  c  as  k,  =b;/c/(/e. 
At  the  wall,  the  pressure  wave  gives  rise  to  an  unsteady  bound¬ 
ary  layer  with  a  thickness  characterized  by  and  a  vdocity 
gradient  of  order 
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The  local  average  rate  of  unsteady  viscous  dissipation  is  propor¬ 
tional  to  the  square  of  the  velocity  gradient.  The  extent  of  the 
free  stream,  characterized  by  the  axial  chord,  is  much  larger  the 
unsteady  boundary  layer  thickness.  To  obtain  the  integrated  ef¬ 
fect  of  dissipation  in  the  free  stream  and  the  unsteady  boundary 
layer,  these  different  scales  have  to  be  accounted  for.  Dissipation 
in  the  free  stream  is  negligible  if 


With  the  definition  of  Sf  and  equation  (16),  (18)  is  rewritten  as 


«  1, 


(19) 


where  the  Reynolds  number  is  defined  in  terms  of  mean  free 
stream  quantities  as  Re  =  For  the  Reynolds  numbers 

encountered  in  turbomachinery,  this  inequality  is  always  satis¬ 
fied  and  dissipation  in  unsteady  boundary  layers  dominates  over 
dissipation  in  the  free  stream.  Unless  the  level  of  the  tu- julent 
viscosity  in  the  freestream  far  exceeds  the  one  characterizing  the 
boundary  layer,  inequality  (19)  also  holds  in  turbulent  flow. 


Linearizing  the  kinematic  and  the  thermodynamic  variables 
around  their  time-mean  (and  nondimensionaiizing  them  by  ti.e 
local  mean  velocity  and  density  in  the  free  stream,  the  blade 
passing  frequency,  and  the  blade  axial  chord),  the  first-order 
perturbation  equations  of  the  streamwise  momentum  equation 
and  the  stress-strain  relations  take  the  form 


The  laminar  viscosity  ^  has  been  assumed  constant.  The  second- 
order  unsteady  dissipation  function  which,  unlike  its  firsi- 
otdei  counterpart,  has  a  nonzero  time- mean,  is 


Theory  -  Unsteady  Viscous  Losses 

The  initial  step  in  the  analytical  approach  is  to  linearize  the 
governing  viscous  equations  around  their  time-mean.  Taking 
the  firzt-order  perturbation  equations  to  their  asymptotic  high- 
frequency  limit,  the  momentum  equation  yields  the  unsteady 
streamwise  velocity  distribution.  From  the  unsteady  velocities  in 
the  boundary  layer,  the  time-mean  of  the  dominating  term  in  the 
second-order  unsteady  dissipation  function  $]  can  be  calculated 
and  related  to  a  change  in  the  time-mean  isentropic  efficiency. 

The  problem  is  described  in  a  Cartesian  coordinate  system 
(s,n)  where  the  coordinate  n  is  aligned  with  the  local  surface 
normal.  The  governing  equations  are  the  continuity  equation, 
the  streamwise  momentum  equation  in  conjunction  with  the 
stress-strain  relations, 


the  energy  equation  in  conjunction  with  the  definition  of  the 
dissipation  function  i, 


and  the  equation  of  state.  Due  to  the  boundary  layer  character 
of  the  problem,  the  normal  momentum  equation  is  replaced  by 
the  assumption  (9p/3n  =  0). 


The  high-frequency  limit  is  characterized  by  a  reduced  frequency 

k.  >  1.  (27) 

As  a  consequence, 

f//?  <  1,  (28) 

A//C  <  1,  (29) 

and  Sf/Xf  <  1  (30) 

hold  in  the  high-frequency  limit.  Inequality  (28)  implies  that 
normal  derivatives  of  unsteady  quantities  dominate  over  normal 
derivatives  of  steady  quantities,  because  the  unsteady  boundary 
layer  thickness  is  much  less  than  the  thickness  of  its  steady  coun¬ 
terpart.  Inequality  (29)  allows  one  to  neglect  streamwise  mean 
flow  derivatives  because  the  wavelength  of  the  disturbance  is 
much  less  than  the  axial  chord  which  characterizes  mean  flow 
gradients.  In  other  words,  parallel  mean  flow  is  assumed  and  the 
streamwise  evolution  of  the  steady  boundary  layer  is  neglected. 
Inequality  (30)  implies  that  the  wavelength  of  the  disturbance 
U  far  greater  than  the  unsteady  boundary  layer  thickness.  It 
allows  for  the  neglect  of  unsteady  streamwise  derivatives  in  the 
presence  of  unsteady  normal  derivatives;  it  always  holds  for  the 
large  Reynolds  numbers  found  in  turbomachinery. 

Taking  the  momentum  equation  (24)  and  the  stress-strain 
relations  (25)  to  their  asymptotic  high-frequency  limits,  the  fol¬ 
lowing  terms  drop  out: 

(27) 


(30)  (29)  (29)  (28) 
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Numbers  above  the  arrows  refer  to  inequalities  (27)  to  {30) 
and  give  the  rationale  for  the  neglect  of  individual  terms.  In¬ 
equality  (27)  leaves  the  unsteady  term  dominant  on  tc  left-hand 
side  of  the  momentum  equation.  On  the  right-hand  side,  the 
driving  unsteady  pressure  term  and  the  unsteady  shear  term 
remain.  Substituting  the  stress-strain  relation  into  the  high- 
frequency  limit  of  the  momentum  equation,  the  velocity  distri¬ 
bution  near  the  wall  is  described  by  the  differential  equation 

di  1  1  ^  „  1  ^ 

‘•at  5s' 

Since  (31)  is  a  linear  constant-coefficient  equation,  smy  distur¬ 
bance  can  be  decomposed  into  its  temporal  Fourier-modes.  For 
the  mth  mode,  the  solution  near  the  wall  is 

^  ~ vf  J})  ^ 

where  p,n  denotes  the  complex  magnitude  of  the  rn-th  harmonic 
of  the  pressure  disturbance.  Figure  4  compares  this  solution  to 
the  exact  solution  for  an  upstream  propagating  pressure  wave 
over  a  fiat  plate  with  a  laminar  boundary  layer.  In  the  Stokes 


-  exact  solution  ii  kg  =  9 

-  equation  (32)  Af*  =  0.3 

Figure  4;  Unsteady  velocity  distribution  as  a  function  of  (n/5/) 


sublayer  near  the  wall,  the  velodty  distribution  is  well  repre¬ 
sented  by  equation  (32).  The  agreement  increases  with  the  ra^ 
tio  SjSf  (which  it  proportional  to  a  laminar  flat  plate 

boundary  layer).  It  decreases  with  the  Mach  number  because 
(32)  assumes  negligible  mean  flow  at  the  edge  of  the  unsteady 
boundary  layer  and  incompressible  flow.  Doe  to  the  neglect  of 
free  stream  convection  it  asymptotes  to  the  wrong  velocity  mag¬ 
nitude  in  the  &ee  stream. 

The  dominating  term  of  the  unsteady  dissipation  function 
$2  io  4he  asymptotic  high-frequency  linut  is  obtained  from  (26) 
with  inequalities  (28)  to  (30).  Taking  its  time-mean  gives 


due  to  the  orthogonality  of  sines  and  cosines.  The  time-mean 
of  the  unsteady  dissipation  is  calculated  mode  by  mode  and 
summed  subsequently.  For  the  m-th  mode,  the  dissipation  per 
unit  volume  is 


(34) 


The  dissipation  is  largest  in  the  region  of  steep  gradients  near 
the  wall  where  the  agreement  between  high-frequency  biiiit  and 
exact  solution  is  best.  Therefore,  equation  (34)  is  integrated 
from  the  wall  into  the  free  stream  to  obtain  the  dissipation  per 
unit  surface  area. 


2mk.p„  V  2mk,'p„  Re  j  ds 


The  difference  in  the  integrated  dissipation  bet  a  eer.  tiie  ex¬ 
act  solution  and  the  high-frequency-limit  solution  of  figure  4  is 
about  10%  of  the  exact  solut.  rn.  The  effects  of  turbulent  bound¬ 
ary  layers,  mean  flow  pressure  gradients,  and  mean  flow  Mach 
number  are  addressed  later. 

Unlike  in  the  Unsteady  Circulation  Losses,  the  energy  is  not 
dissipated  downstream  of  the  last  blade  row  but  in  the  blade 
passage.  The  mean  flow  state  al  mg  the  blade  varies  and  the 
entropy  rise  per  unit  dissipated  energy  is  not  constant.  Integra- 
ticn  around  the  blade  surface  and  division  by  the  mass  flux  rh 
yiela  the  average  specific  entropy  rise  in  the  blade  passage  due 
to  unsteady  viscous  dissipation.  In  dimensional  quantities  it  is 


(3«) 

The  average  specific  entropy  rise  is  again  related  to  an  efficiency 
drop,  as  in  equation  (14). 


Results  -  Unsteady  Viscous  Losses 


The  model  has  been  applied  to  the  ACE^-turbine  stage  described 
in  the  results  section  of  the  Unsteady  Circulation  Losses.  The 
input  required  for  the  model  is  the  time  harmonics  of  the  lo¬ 
cal  pressure  gradient  along  the  blade  surfaces.  Those  can  be 
the  result  of  a  ^scous  or  an  inviscid  unsteady  simulation.  The 
dominating  unsteady  effect  of  rotor  blade  row/shock  wave  in¬ 
teraction  is  captured  in  an  inviscid  simulation.  The  parameters 
of  the  simulation  were  similar  to  the  ones  used  before  except 
that  UNSFLO  was  run  in  the  viscous  mode  (Giles  and  Haimes, 
1991a),  the  rotor  speed  was  increased,  and  the  back  pressure 
was  lowered  slightly.  This  resulted  in  stronger  stator  trailing 
edge  shock  waves  and  reduced  the  mean  rotor  incidence  angle 
to  48°,  matching  the  experiments  of  Johnson  et  al.  (1988)  .  In 
experiments  (Ashworth  et  al.,  1985)  ,  the  rotor  boundary  layer 
was  laminar  over  most  of  the  suction  snrface  for  low  levels  of  free 
stream  turbulence  and  in  the  absence  of  wakes  and  shock  waves. 
The  pressure  surface  boundsiry  layer  was  transitional.  Bar  pass¬ 
ing  experiments  to  simulate  the  stator  wakes  and  shock  waves 
showed  that  widely  spaced  wakes  and  shock  waves  temporarily 
trip  the  boundary  layer,  which  will  relax  back  to  a  laminar  state 
after  the  disturbance  has  passed.  For  a  bar  spacing  resembling 
engine  conditions,  the  boundary  layer  is  transitional  or  turbu¬ 
lent.  In  the  simulation,  the  flow  was  prescribed  turbulent  from 
the  leading  edge  on. 
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Figure  5  shows  the  specific  enttop)  rise  due  to  unsteady  vis¬ 
cous  dissipation  on  a  rotor  blade  surface.  A  value  of  j=0  cor 
responds  to  the  rotor  leadbng  edge;  positive  and  negative  values 
correspond  to  pressure  and  suction  surface,  respectively  The 
losses  are  concentrated  on  the  suction  surface  between  rotor 
leading  edge  and  crown,  where  the  pressure  gradient  fluctuates 
in  response  to  shock  wave  passing.  The  integrated  entropy  rise 
for  a  rotor  blade  is  Aj  =  this  translates  into  an  effi 

ciency  drop  of  A»j,,r/i?,,r,o=-0.09%.  These  numbers  are  subject 
to  errors  examined  in  the  following  section.  Of  the  total,  64%  is 
at  blade  passing  frequency  and  14%  in  the  first  harmonic. 


Figure  6:  Entropy  rise  in  the  steady  and  the  unsteady  nmaiaiion 

Figure  6  shows  the  time-mean  mass-averaged  entropy  rise 
through  the  turbine  stage  in  the  simulations.  The  ajdal  coor¬ 
dinate  X  if  nondim#n«onali*ed  by  the  stator  axial  chord;  the 
stator  leading  edge  is  located  at  *  =  0.  Three  observations  can 
be  made.  First,  the  steady  simulation,  in  which  case  the  con- 
servatirn  variables  are  averaged  at  the  interface,  has  a  higher 


entropy  rise  through  the  stage  ihan  the  uBs’.esCy  suuaiaUt.!.  i; 
the  steady  case,  '.he  fluxes  of  mass,  momenSum  arc  er.etgi  are 
averaged  at  the  stator  outflow  boundary  to  provide  uniform  and 
steady  rotor  inflow  conditions  The  (nonphysical)  eniiopy  rue 
at  the  interface,  accounts  lot  about  20%  of  the  total  At  a  eonse 
quec  the  iientropic  stage  efheseacy  is  0  4%  lets  m  the  steady 
simuiation  Since  the  unsteadiness  does  nol  propagate  beyond 
the  steady  interface,  the  entropy  rises  much  slower  m  the  rotor 
passage  than  in  the  unsteady  simulation  Second,  the  entropy 
rise  between  rotor  leading  edge  and  crows,  the  region  where  the 
unsteady  losses  are  concentrated,  is  about  10*’  for  the  steady 
case  and  3*10*’  lot  the  unsteady  case  The  difference  is  much 
larger  than  the  entropy  nse  that  has  been  aseribed  to  unslead.v 
VISCOUS  losses.  Third,  the  entropy  rise  between  crown  and  trail 
mg  edge  as  well  as  between  Iraiiing  edge  and  outlet  boundicy 
is  significantly  in  the  unsteady  sunuiatioc  Banag  turbulence 
modebng  in  the  boundary  layer  as  a  reason  for  the  steeper  nse. 
numerical  smoothing  acting  on  the  unsteadinesi  m  the  blade 
passage  remains  as  the  only  explanation 

Modeling  Accuracy 

Loss  model 

The  derivation  of  the  unsteady  viscous  losses  involved  (our  ap¬ 
proximations,  each  of  which  introduced  enors  The  first  two.  the 
boundary  layer  approximation  and  the  liDeaniaiion  thereof,  will 
not  be  examined.  The  errors  introduced  in  the  last  two  steps, 
the  invocation  of  the  high  freqnenry  limit  and  the  assumption 
of  laminar  flow,  dominate. 

In  laminar  flow,  the  two  main  effects  not  modeled  is  the 
high-frequency  limit  are  the  compression  of  the  boundary  layer 
and  the  free  stream  convection.  Even  larger  errors  are  due  to  the 
presence  of  a  turbulent  mean  flow.  To  estimate  the  errors,  a  lO 
code  was  written  to  solve  the  linearized  equations  for  a  flat  plate 
boundary  layer,  figure  7,  subject  to  a  pressure  wave  propagating 
at  an  arbitrary  nonzero  angle  9  relative  to  the  surface  normal. 
The  parallel  flow  assumption  was  retained  because  the  evoln- 


Figure  7:  Model  problem  for  nomericai  error  evaiaation 

tion  of  the  mean  flow  boundary  layer  takes  place  over  a  length 
scale  much  larger  than  the  steady  or  the  unsteady  boundary 
layer  thickness.  This  astumption  combined  with  a  small  steady 
and  unsteady  boundary  layer  thickness  (compared  to  the  wave¬ 
length),  reduces  the  problem  to  one  dimension.  Tbe  boundary 
was  modeled  as  an  adiabatic  wall. 
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All  mean  Telocity  proftlet  were  t&kes  to  be  idcnlickl  to  the 
tresponding  profiles  in  ste»dy  flow.  Only  very  large  dislur- 
uance  amplitudes  (and  frequencies  close  to  the  turbulent  burst¬ 
ing  frequency)  lead  the  mean  flow  profiles  to  deviate  from  the 
corresponding  profiles  in  steady  flow  (Karlsson,  1959;  Parikh  et 
al.,  1982)  .  In  laminar  flow,  the  velocity  profiles  were  obtained 
as  the  station  of  the  Falkner-Skan  equations.  In  turbulent  flow, 
a  three  layer  model  with  the  inner  law  in  Spalding’s  formulation 
(Spalding,  1961)  and  a  law  of  the  wake  in  Coles’  formulation 
(Coles,  1956)  was  used.  Both  the  laminar  and  the  turbulent 
velocity  profile  are  strictly  valid  for  incompressible  flow  only, 
but  were  used  r^ardless  of  the  free  stream  Mach  number.  The 
mean  temperature  and  density  profiles  were  obtained  from  the 
velocity  profiles  with  the  assumption  of  a  unit  Prandtl  number. 

The  laminar  pari  of  the  code  was  verified  by  comparison  to 
Lighthill’s  analytical  work  (1954)  on  oscillatory  boundary  layers 
in  an  incomprestibie  mean  flow.  The  turbulent  part  was  first 
checked  against  the  laminar  solution  for  very  high  fiequencies 
and  then  agairut  experiments  by  Cousteix  et  al.  (1977). 

The  numerical  error  evaluation  is  summarised  in  figures  8 
to  10.  Figure  8  assumes  laminar  mean  flow  and  pressure  waves 
propagating  parallel  to  the  boundary.  It  shows  the  ratio  of  the 
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Figure  8:  High-frequency-limit  dissipation  and  full  linearised 
laminar  model  dissipation 


integrated  disnpatkm  ¥](  calculated  in  the  high-frequency  limit 
(HFL)  to  the  one  calculated  with  the  foQ  linearued  equations. 
In  the  limit  (27.  — >0)  or  (k,->oo),  the  HFL  model  pves  the  cor¬ 
rect  diasipatioa.  The  lower  the  rednced  freqnency  for  a  nonsero 
free  stream  Mach  number,  the  greater  the  deviation.  This  is 
a  consequence  of  the  unsteady  boundary  layer  spreading  away 
from  the  waU  into  the  steady  boundary  layer  where  convective 
terms  in  (24)  become  important.  Depending  on  the  direction  of 
propagation,  the  HFL  yields  an  integrated  dissipation  greater  or 
less  than  the  full  linearued  model. 

The  flows  in  tnrbomachinery  are  almost  always  transitional 
r  turbnient  and  the  unsteady  boundary  layer  extents  beyond 
ne  laminar  sublayer.  In  figures  9  and  10  errors  doe  to  the  as¬ 
sumption  of  laminar  mean  flow  in  the  HFL  model  appear  super¬ 
imposed  on  errors  due  to  finite  reduced  frequencies.  They  show 


Figure  9:  High-frequency  limit  dissipation  and  full  turbulent 
model  dissipation  ■  downstream  ptopagaiiag  wave 


Figure  10:  High-frequency  Umit  dissipation  and  full  turbulent 
model  dissipatioD  -  upstream  propagating  wave 

the  ratio  of  the  integrated  dissipation  calculated  in  the  HFL  to 
the  full  turbulent  dissipation  for  a  Reynolds  number  ite,=10*. 
As  opposed  to  the  case  in  laminar  flow,  the  HFL  is  not  exact 
in  the  Emit  (XT.-rO).  ’This  is  because  the  turbulent  nature  of 
tbe  mean  flow  introduces  additional  errors  into  the  HFL  model 
which  is  an  inherently  laminar  model.  For  high  reduced  fre- 
quendet,  the  turbulent  dissipation  again  approaches  the  value 
of  the  EIFL  model  although  this  is  not  evident  within  the  range 
M  reduced  frequencies  displayed.  The  parameter  (2i/f/)  is  the 
ratio  of  the  thiekoeH  of  the  laminar  sublayer  of  the  mean  flow, 
defined  es  SOy'*'  in  wall  coordinates,  to  the  thickness  of  the  un¬ 
steady  boundary  layer.  For  frequencies  such  that  (iilS/)  >  1, 
the  unsteady  boundary  layer  is  confined  to  tbe  laminar  sublayer 
and  the  HFL  model  gives  accurate  results. 

The  ratio  (ttijiriliii)  for  a  downstream  propagating  wave 
in  turbulent  flow  behaves  in  a  dmilar  fashion  to  the  correspond¬ 
ing  ratio  in  laminar  flow.  With  higher  reduced  frequency  and 
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l^er  Ma^  number,  the  HFL  is  mote  accurate.  The  behavior  of 
for  an  upstream  propagating  wave,  however,  dif¬ 
fers  significantly  from  the  behavior  in  laminar  flow.  In  the  limit 
(M,->0),  the  ratio  is  less  than  unity  and  matches  the  value  for  a 
downstream  propagating  wave.  The  laminar  HFL  model  under- 
predicts  the  velocity  gradients  (and  uses  the  laminar  viscosity) 
and  thus  arrives  at  too  low  an  integrated  dissipation  for  finite 
reduced  frequencies  (or  h/Sf).  As  the  Mach  number  increases, 
this  ratio  rises  above  unity  because  the  HFL  overpredicts  the 
unsteady  velocity  magnitude  in  the  free  stream;  it  is  the  rela¬ 
tive  magi^ude  of  these  two  effects  that  governs  the  behavior  of 
in  figure  10. 

The  flow  through  the  rotor  of  the  ACE  stage  is  turbulent  and 
there  is  a  strong  favorable  mean  pressure  gradient  in  the  region 
of  unsteady  viscous  losses.  Calculating  the  laminar  sublayer 
thickness  from  the  wail  shear  of  the  simulation,  gives  a  ratio 
lilSfKl.2.  The  free  stream  Mach  number  in  the  rotor  frame  is 
(^,  a:  0.65)  in  the  main  region  of  unsteady  viscous  losses.  There¬ 
fore,  figures  9  and  10  suggest  that  the  unsteady  viscous  losses 
calculated  for  the  ACE  turbine  stage  =  0-09%)  are 

at  most  a  factor  1.5  too  large  and  could  be  up  to  five  times 
larger. 


Numerical  smoothing 


Consider  a  pressure  wave  propagating  through  an  unloaded  flat 
plate  cascade  with  uniform  mean  flow  17,.  The  incident  pressure 
wave  of  figure  7,  propagating  at  an  angle  8  relative  to  the  surface 
normal,  is  described  by  its  wave  amplitude 


exp  {  t  k  (s  tin  -  n  cos  6)  -  it}.  (37) 

where  k  denotes  the  wavenumber.  Upon  reflection  at  the  blade 
surface,  the  wave  will  lose  a  small  fraction  fi  of  its  amplitude  to 
dissipation  in  an  unsteady  boundary  layer.  The  time-mean  rate 
of  energy  dissipation  per  unit  surface  area  it  given  by  equation 
(35)  with  the  pressure  gradient  determined  from  (37).  The  time- 
mean  rate  of  acoustic  enetQr  incident  on  the  wall  is  given  by  the 
acoustic  intensity 


**  I  lF«i  «  j  IpI*  co*^.  (38) 

The  rate  of  acoustic  energy  reflected  is 

7,  i  (1  -  /9)’  Ipfil  «  I  (1  -  Ipl’  ST.  CO.  9.  (39) 


The  above  acoustic  intensities,  strictly  valid  only  for  ST,  «  0,  are 
used  for  simplicity.  Since  the  objective  is  an  order-of-magnitude 
argument,  these  simplifleation  do  not  compromise  the  subse¬ 
quent  condusioiM.  The  streamwise  velocity  distribution  in  the 
unsteady  boundary  layer  induced  by  the  incident  and  the  re¬ 
flected  wave  are  nearly  in  phase  because  of  inequsJity  (30).  The 
rate  of  unsteady  viscous  dissipation  per  unit  surface  area  it 


¥k  = 


(2-^)*  /TX 

2fc,  y2k.Re  d»  ’ 


(«} 


where  incompressible  mean  flow  was  assumed.  Balancing  incom¬ 
ing  intensity,  reflected  intensity,  and  unsteady  dissipation  gives 
an  estimate  for  the  attenuation  factor  P  of 


0\  =  7?,  ^2  k,IRe  tan  sin  6,  (41) 


where  the  lektioo  valid  in  the  limit  (T?,-— 0),  has 

been  used.  Alternaiively,  (41)  can  be  expressed  as 

01  =  ^2  (M./Re)  (2»/A/)  tan  6  an  tf,  (42) 


The  longest  wave  in  a  turbomachine  typieaiiy  has  a  nondimen- 
iional  Wavelength  Ay  of  order  unity.  Assuming  Reynolds  and 
Mach  numbers  typical  of  turbomacKinery  (and  tan  #1108  of  or¬ 
der  one),  0(10’)  reflections  should  occur  before  the  wave  am¬ 
plitude  is  halved.  In  an  invitcid  simulation  of  the  flow  through 
the  ACE  turbine  stage  no  more  than  three  reflections  of  weak 
shock  waves/pressute  waves  are  discernible.  The  approxima¬ 
tions  made  in  the  derivation  of  the  HFL  model  and  in  the  above 
order-of-magnitude  argument  can  not  account  for  a  difference  of 
two  orders  of  magnitude  This  suggests  that  numerical  smooth¬ 
ing  in  the  free  stream  dominates  over  dissipation  in  unsteacy 
boundary  layers. 

The  propagation  of  a  pressure  wave  in  the  uniform  free  stream 
can  be  modeled  numerically  by  the  convection-diffusion  equation 


dt 


+  (l±77.sm«>)|?  =  c,d|i-c,A*0,  (43) 


where  s  is  a  coordinate  in  the  direction  of  wave  propagation;  A 
denotes  the  mesh  spacing.  The  time  t  has  been  nondimention- 
alized  by  the  asaal  chord  and  the  speed  of  tonnd.  Numerical 
smoothing,  necessary  to  stabilite  the  code,  has  been  added  on 
the  right-hand  side.  The  time-rate  of  decay  of  a  wave  ampUtnde 
due  to  second-order  smoothing  for  shock  wave  capturing  it 


ft  ~  « j  k*  A.  (44) 

In  smooth-flow  regions,  second-order  smoothing  is  turned  off. 
There,  the  time-rate  of  decay  due  to  fourth-order  smoothing  it 

ft  ~  «4  k*  A*.  (45) 


Assuming  that  a  wave  typically  travels  one  axial  chord  between 
tnccestive  reflections,  the  time  between  reflections  is  t  a;  1  for 
moderate  Mach  numbers.  Therefore,  ft  and  ft  also  represent 
the  amplitude  attenuation  caused  by  numerical  smoothing  be¬ 
tween  reflections.  Alternatively,  they  can  be  written  as 

ft  ~  A  (2v/A/)’  and  ft  ~  €4  A*  (2v/Ay)*  (46) 


which  shows  the  strong  bias  of  fourth-order  smoothing  towards 
short  wavelengths.  Typical  values  of  the  smoothing  coeffidents 
are  (tj  =  10“*)  and  (.4  =  10"*).  A  standard  value  of  the  mesh 
spadng  is  As:3«10~*.  For  the  longest  vrave  with  a  wavelength 
of  order  unity,  fourth-order  smoothing  in  between  successive  re¬ 
flections  results  in  an  amplitude  attenuation  that  is  an  order 
of  magnitude  smaller  than  the.  attenuation  caused  by  viscous 
dissipation  in  an  unsteady  boundary  layer.  For  shorter  wave¬ 
lengths,  numerical  smoothing  starts  to  dominate.  The  attenoa- 
tion  coeffident  ft  of  the  sixth  harmonic  reproaches  nnity  for  the 
above  choice  of  A  and  £4.  This  implies,  that  numerical  smooth¬ 
ing  dissipates  waves  prematurely  in  the  free  stream  unless  the 
wavelength  is  on  the  order  of  an  axial  chord  or  more. 

The  weak  shock  waves,  whose  interaction  with  the  blade  rows 
is  the  primary  cause  of  unsteadiness  in  the  ACE  turbine  stage, 
are  subject  to  strong  numerical  dissipation,  even  in  the  absence 
of  second-order  smoothing.  The  pressure  gradients  driving  the 
unsteady  boundary  layer  '  -e  diminished  by  nnmerical  smooth¬ 
ing,  lowering  the  level  of  the  predicted  unsteady  viscous  losses. 
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Conclusions 

This  paper  has  studied  two  unsteady  loss  mechanisms  for  ge¬ 
ometries,  Mach  numbers  and  loading  realistic  for  modem  tur¬ 
bomachinery.  The  results  of  numerical  simulations  were  used 
as  input  to  the  2D,  linearized,  analytical  loss  models.  Losses 
from  both  mechanisms  are  small  compued  to  other  (nonlinear) 
sources  that  increase  the  losses  under  unsteady  operating  con¬ 
ditions,  like  a  change  in  the  transition  point,  and  compared  to 
the  uncertainty  in  measured  efficiencies.  Since  this  paper  is  2D, 
it  does  not  account  for  possibly  important  3D  effects. 

The  Unsteady  Circulation  Losses,  which  are  a  due  to  vis¬ 
cous  dissipation  of  kinetic  energy  associated  with  vorticity  shed 
in  response  to  circulation  variations,  resulted  in  an  isentropic 
efficiency  drop  of  0.3%  for  the  ACE  turbine  stage.  The  ACE 
turbine  stage  operates  at  a  comparably  high  level  of  unsteadi¬ 
ness  but  is  representative  of  current  and  future  high-pressure 
turbine  designs. 

The  Unsteady  Viscous  Losses,  i.e,  the  additional  time-mean 
dissipation  in  unsteady  boundary  layers,  that  were  calculated  for 
the  ACE  turbine  stage,  are  much  smaller  and  caused  a  drop  of 
0.09%  in  the  isentropic  efficiency.  Their  exact  magnitude  is  un¬ 
certain  due  to  the  limitations  of  the  model  used  and  the  effect  of 
numerical  smoothing  in  the  simulation.  With  the  bounds  found 
in  a  numerical  error  study,  they  can  be  of  the  same  magnitude 
as  the  Unsteady  Circulation  Losses  for  this  stage. 

In  unsteady  simulations,  numerical  smoothing  can  shift  the 
locus  of  dissipation  from  unsteady  boundary  layers  into  the 
blade  passages.  Pressure  waves,  which  are  the  input  to  the  vis¬ 
cous  loss  model,  axe  diminished  prematurely  in  the  passage.  In 
steady  simulations,  the  entropy  rise  at  the  stator/iotor  interface, 
which  is  an  numerical  artifact,  can  account  for  a  significant  part 
of  the  total  losses  in  a  stage. 
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An  Approach  for  Multi-Stage  Calculations 
Incorporating  Unsteadiness 
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Abstract 

This  paper  describes  a  mathematical  approach  to  the  calculation 
of  unsteady  flow  in  multi-stage  turbomachinery.  An  asymptotic 
formulation  is  used,  with  the  small  asymptotic  parameter  be¬ 
ing  the  level  of  unsteadiness  in  each  blade  row.  The  baseline 
flow  is  the  nonlinear  steady  flow  that  is  computed  by  many  ex¬ 
isting  multi-stage  cakulation  methods.  The  flrst  order  linear 
perturbation  is  the  unsteady  flow  field  arising  from  stator/rotor 
interactions  between  neighboring  blade  rows.  The  second  order 
correction  contains  the  information  about  the  time-averaged  ef¬ 
fect  of  unsteadiness  on  the  mean  flow.  The  advantage  of  this 
asymptotic  approach  is  that  it  leads  to  a  set  of  equations  which 
can  be  s^ved  numerically  very  much  more  cheaply  than  the  full 
nonlinear  unsteady  flow  equations,  while  still  retaining  the  key 
features  of  the  flow. 

Introduction 

During  the  last  ten  years  many  researchers  have  developed  nu¬ 
merical  methods  for  calculating  nonlinear  unsteady  flows  (Chen, 
1988,  Erdos,  1977,  Fouimauz,  1986,  Giles,  1988,1990,  He,  1991, 
HodMtt,  1984,  Rai,  1987,1989,  Siden,  1991).  The  capabilities 
of  the  codes  which  have  been  developed  vary  widely.  Some 
are  two-dimensional,  some  have  quasi-three-dimennonal  exten¬ 
sions,  and  some  are  fully  three-dimensional.  Some  assume  tn- 
visdd  flow  snd  use  the  Euler  equations,  while  others  adopt 
the  Reynolds-averaged  Navier-Stokes  equations  with  turbulence 
modeling.  Some  are  able  to  treat  arbitrary  number  of  blades  in 
each  blade  row  using  techniques  due  to  Erdos  (1977)  or  Giles 
(1988)  but  most  treat  only  very  simple  integer  ratios,  such  as 
2:3.  Most  can  analyze  stator/rotor  interaction  in  which  the  un¬ 
steadiness  is  due  to  the  relative  blade  motion.  Some  can  also 
treat  unsteadiness  due  to  prescribed  blade  vibration  (He,  1991, 
Siden,  1991,  Giles,  1991)  incoming  wakes  (Hodson,  1984,  Giles, 
1988)  or  potential  disturbances  (Giles,  1991). 


These  codes  are  now  very  useful  practical  tools  for  the  in¬ 
vestigation  of  important  unsteady  flow  phenomena  in  turboma¬ 
chinery.  This  is  clear  from  the  growing  number  of  papers  con¬ 
centrating  not  on  the  mathematical  formulation  of  the  methods 
but  on  their  use  as  ‘numerical  wind  tunnels’.  Sometimes  this 
use  is  in  conjunction  with  detailed  experimental  work,  in  which 
case  the  computation  is  used  to  understand  the  experimental 
results  more  fully  (Abhari,  1991,  Epstein,  1989,  Dorney,  1991, 
Johnson,  1990).  The  experimental  data  is  usually  limited  to  a 
few  variables  (e.g.  pressure,  total  temperature)  in  a  few  loca¬ 
tions,  whereas  the  computational  results  supply  all  variables  at 
all  locations.  In  other  cases,  the  confidence  in  numerical  meth¬ 
ods  u  sufficiently  great  that  computations  are  performed  in  the 
complete  absence  of  experimental  data  due  to  the  difficulty  or 
cost  of  performing  the  experiments  (Paulon,  1991).  This  is  par¬ 
ticularly  true  for  parametric  design  studies  (Rangwaila,  1991, 
Korakianitis,  1991). 

Despite  the  impressive  capabilities  of  these  nonlinear  time- 
marching  methods  there  is  a  major  limitation  for  industrial  en¬ 
gineering  application,  which  is  the  cost  and  compute  times  as¬ 
sociated  with  these  codes.  Invisdd,  two-dimensional,  unsteady 
computations  can  be  performed  in  minutes  on  a  CRAY,  or  less 
than  an  hour  on  today’s  new  generation  of  very  fast  worksta- 
tiona.  At  the  other  extreme  however,  viscous  three-dimensioual 
unsteady  calculations  can  require  several  days  of  CRAY  time. 
These  computational  costs  can  be  justified  when  trying  to  iden¬ 
tify  and  correct  design  problems.  They  are  also  justifiable  when 
the  codes  are  used  as  a  research  tool  to  increase  our  understand¬ 
ing  of  unsteady  flows  in  turbomachinery,  so  that  our  increased 
knowledge  can  lead  to  improved  designs  with  fewer  problems. 
However,  axial  gas  turbine  engines  are  multi-stage  in  nature, 
and  if  one  considers  multi-stage  calculations,  possibly  even  the 
calculation  of  the  flow  through  an  entire  engine,  then  nonlinear 
time-marching  methods  teem  to  be  totally  impractical  for  the 
foreseeable  future. 

It  is  not  simply  that  the  computational  work  scales  linearly 
with  the  number  of  stages.  It  may  in  fact  grow  faster  than 
linearly  because  of  the  possibility  of  slower  convergence  due  to 
the  increasing  infiow-outflow-inflow  round-trip  signal  propaga¬ 
tion  time  in  unchoked  engines.  The  main  complicating  factor 
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leading  to  extraordinarily  large  computing  times  is  the  problem 
of  the  periodic  boundary  conditions.  In  single-stage  computa¬ 
tions  there  is  a  time-lagged  spatial  periodicity  condition  which 
can  be  treated  numerically  by  one  of  two  approaches  (Erdos, 
1977,  Giles,  1988).  As  soon  as  there  is  more  than  one  stage  in 
the  calculation  there  is  an  unavoidable  problem.  Consider  a  1 1 
stage  calculation  in  which  the  first  and  second  stator  rows  have 
different  numbers  of  blades.  If  there  is  no  common  factor  in  the 
numbers  of  blades  then  each  of  the  second  stator  row  blades  will 
experience  a  different  unsteady  force  depending  on  its  circum¬ 
ferential  position  relative  to  the  blades  in  the  first  stator  tow. 
There  is  no  mathematically  correct  way  around  this  problem. 
Some  researchers  have  performed  calculations  by  modifying  ge¬ 
ometries  to  achieve  very  simple  ratios  in  the  numbers  of  blades 
(e.g.  1;1,2:3,3:4}  but  these  modifications  change  the  flow  field 
being  calculated,  and  involve  larger  and  larger  errors  as  more 
stages  are  included  in  the  calculation.  I'lrimately,  the  only  cor¬ 
rect  nonlinear  time-marching  treatment  is  to  analyze  the  entire 
annulus  in  each  blade  row.  This  is  simply  too  computationally 
expensive  to  consider  for  many  years. 

The  question  then  is  how  one  can  best  perform  multi-stage 
flow  calculations  in  a  manner  which  includes  the  prediction  of 
unsteady  effects.  Before  offering  an  answer  to  this,  it  is  appro¬ 
priate  to  flist  discuss  what  one  would  hope  to  get  out  of  such 
calculations.  The  focus  in  this  work  is  on  the  development  of  a 
tool  that  can  be  used  by  designers  as  part  of  the  regular  design 
system.  From  the  designer’s  point  of  view  one  of  the  principal 
objectives  is  an  accurate  prediction  of  mean  flow  properties  such 
as  mass  fliow,  turning,  work  and  efficiency.  This  is  the  aim  of 
research  work  by  Adaraczyk,  who  in  the  early  1980’s  pointed  out 
I  the  impracticality  of  multi-stage  unsteady  Navier-Stokes  calcu¬ 

lations.  His  point  is  that  the  unsteady  components  of  the  flow 
are  important  only  in  as  much  as  they  change  the  mean  flow. 
Therefore,  he  formulated  a  system  of  ‘passage-averaged’  equa¬ 
tions  in  which  the  effect  of  deterministic  periodic  unsteadiness 
on  the  mean  flow  is  included  through  terms  which  are  similar 
in  nature  to  the  Reynolds  stresses  in  the  Reynolds-averaging  of 
turbulent  flow  equations  (Adamczyk,  1985).  The  implementa¬ 
tion  of  these  equations  in  a  multi-stage  flow  code  has  been  very 
successful  in  predicting  the  mean  flow  in  both  turbines  (Adam- 
czyk,  1989)  and  compressors  (Mulac,  1991).  'In  particular,  his 
f  approach  has  correctly  predicted  flow  features  which  are  not  cap- 
*  toted  by  steady-state  multi-stage  flow  calculations  which  do  not 
£  include  unsteady  flow  effects.^These  ideas  of  Adamczyk  have 
^  played  a  major  role  in  stimulating  the  ideas  for  the  quadratic 
jtenns  in  the  approach  in  this  paperi'.',“  7-.  , 

The  other  major  concerns  with  flow  unsteadiness  are  the 
possibility  of  bending  or  torsional  flutter,  and  the  magnitude 
of  unsteady  loads  and  forced  response  due  to  various  kinds  of 
unsteady  interactions.  Flntter  is  by  definition  a  linear  problem 
rince  it  is  concerned  with  the  possible  growth  of  infinitesimal  un¬ 
steady  perturbations.  A  growing  body  of  evidence  suggests  that 
the  forced  response  problem  is  also  largely  linear.  Comparisons 
between  linear  and  nonlinear  calculations,  for  levels  of  unsteadi¬ 
ness  typical  of  engine  conditions,  show  that  the  linear  calcula¬ 
tions  predict  the  forced  response  to  within  5%  (Hall,  1989,  Sud- 
dhoo,  1991).  This  level  of  accuracy  is  perfectly  adequate  for  all 
practical  engineering  design  purposes  and  is  probably  within  the 
level  of  other  modeling  errors.  This  has  motivated  the  develop¬ 
ment  of  linear  Euler  methods  for  forced  response  prediction.  In 


these  methods,  the  unsteadiness  is  assumed  to  be  a  small  pcrlut- 
bation  to  a  steady  nonbnear  flow  field  obtained  from  a  solution 
of  the  Euler  equations.  Some  preliminary  research  on  a  briear 
Euler  method  was  performed  in  197!)  by  Ni  bus  the  re¬ 

cent  research  activity  began  with  the  PhD  thesis  of  Hall  (1987) 
Hall’s  research  involved  the  development  of  a  two-dimensional 
Anile  element  method  for  the  calculation  of  bnear  unsteady  flows 
in  cascades,  plus  shock-fliting  for  unsteady  flows  in  ducts  (Ball, 
1989).  By  comparing  to  results  of  unsteady,  nonlinear  Euler 
time-marching  calculations  he  demonstrated  the  validity  of  the 
bnear  approximation  up  to  large  levels  of  unsteadiness  Later, 
Giles  proposed  an  alternative  pseudo- time- marching  approach 
to  solving  the  linear  Euler  equations  (Giles,  1989)  based  on  an 
idea  of  Ni  (1975).  Lindquist  and  Giles  have  also  proved  that 
with  careful  treatment  shock  capturing  can  be  valid  for  linear 
unsteady  calculations  (Lindquist,  1991).  Pseudo-time-marching 
shock-capturing  bnear  Euler  methods  are  also  being  developed 
by  Hall  (1991)  and  Holmes  (1991).. 

This  paper  combines  the  ideas  of  Adamczyk  with  the  recent 
developments  in  bnear  unsteady  flow  modebng,  to  formulate  an 
asymptotic  approach  to  multistage  unsteady  flow  calculations. 
First,  one  calculates  a  steady  multi-stage  flow  field  in  a  manner 
similar  to  that  used  by  Denton  (1985,  1990}  and  Dawes  (1990). 
Next,  one  calculates  the  unsteady  flow  assuming  that  its  magni¬ 
tude  is  sufficiently  small  that  it  may  be  analyzed  using  the  bn- 
earized  form  of  the  unsteady  Euler  equations.  In  the  third  step, 
changes  to  the  steady-state  flow  field  due  to  the  time-averaged 
effect  of  the  unsteadiness  are  related  to  source  terms  whose  mag¬ 
nitude  is  quadratic  in  the  level  of  unsteadiness  This  is  similar 
to  the  modebng  of  Reynolds  stresses  in  turbulent  flow,  except 
that  in  this  case  the  unsteady  flow  perturbations  have  been  cal¬ 
culated  using  the  linear  Euler  analysis  and  so  the  constmetion 
of  the  quadratic  source  terms  is  a  purely  algebraic  exercise. 


1  Asymptotic  analysis 


1.1  Unsteady  Euler  equations 


The  starting  point  for  the  entire  analysis  is  the  nonlinear  un¬ 
steady  Euler  equations  describing  the  motion  of  an  invisdd, 
compressible  gas  in  three  dimensions  in  a  frame  of  reference 
rotating  with  angular  velocity  (1  about  the  z-axis. 


dU 

et 


\dz  dy  dz) 


+  5, 


where  the  state  vector  is  a  five-component  vector, 


P 

P^ 

V 

pu, 

\  P^  / 


(1) 


(2) 


the  three  flux  vectors,  F,,  and  F,,  are  given  by 


'  pv, 
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and  the  source  term  $  containing  both  centrifugal  and  coriolis 
effects  is 

/  ° 

0 

5=  pin^y  +  lUu,)  .  (4) 

p(n^z  -2nuj,) 

\  0 

The  total  energy  E  is  related  to  the  pressure,  density  and 
velocities  through  the  equation  of  state  for  am  ideal  gas. 

Integrating  over  a  static  control  volume  gives  the  integral 
form  of  the  equations  of  motion. 

s///v‘'''‘'- 

(6) 

It  is  this  equation  which  is  used  as  the  basis  for  the  discretization 
in  most  numerical  time-marching  methods. 

1.2  Asymptotic  expansions 

The  approach  which  is  proposed  uses  the  first  three  terms  in  an 
asymptotic  expansion  in  which  the  level  of  unsteadiness  is  the 
small  asymptotic  parameter.  It  is  very  closely  related  to  the 
mathematical  techniques  used  to  analyze  a  phenomenon  known 
as  ‘acoustic  streaming’,  in  which  unsteady  sinusoidal  perturba¬ 
tions  of  a  stationary  fluid  can  generate  a  small  second  order 
mean  flow  (Batchelor,  1967,  Schlichting,  1979). 

The  flow  srariables  are  represented  by  an  asymptotic  expan¬ 
sion  in  t. 

U(z,t)  =  -b  -b  0{f^).  (7) 

The  next  step  is  to  construct  the  asymptotic  form  of  all  flow 
variables  using  Taylor  series  expansion.  For  an  arbitrary  scalar 
function 

(8) 

Note  that  the  second  order  term  hu  two  pieces,  one  due  to 
the  linear  sensitivity  to  second  order  changes  in  U  (which  will 
be  referred  to  as  the  second  order  linear  perturbation)  and  the 
other  due  to  the  quadratic  sensitivity  to  the  first  order  changes 
in  (7  (which  will  be  referred  to  as  the  quadratic  perturbation). 

To  avoid  the  use  of  superscripts,  the  superscript  will  be 
dropped  unless  it  is  needed  for  clarity.  In  addition,  the  following 
notation  is  introduced. 

4  =  £A[;0> 

-  dUj  ’ 

'a  .  ^W'  (S) 


The  overlining  on  the  r.h.s  of  these  definitions  denotes  s  time- 
average.  The  reason  for  introducing  the  tune.sveraees  is  that  it 
is  only  the  mean  second-order  effects  which  arc  considered  to  be 
important  in  the  analysis  because  they  affect  quantities  such  as 
the  mean  efficiency  and  mass  flow.  The  unsteady  component  of 
the  second  order  terms  is  assumed  to  give  only  a  small  change  in 
the  unsteady  loading  on  the  blades  and  is  therefore  neglected. 

After  some  algebra  one  finds  that  the  linear  perturbations  to 
the  velocities  and  pressure  are 

-  u.p) 

“»  =  P"'(^y  -  UyP)  (10) 

u,  =  -  “.a). 


=  (7-1) 


(7  -  l)p 


and  hence  the  perturbations  to  the  flux  vectors  are 

^  u.p  -b  ' 

-b  u.pu,  u-  p 

F,  -  u.^y  +  U,pu„ 

+  u.ptt, 

,  U.{p£  -T  p)  -b  i..(pF  -e  p)  ) 

f  UyP  -b  UyP  ' 

+  “vP". 

Fy,  -  t*vP“y  +  UyffUy  +  P 

__UyPU.  -f  ttyPU, 

U,(p£  -b  p)  -b  tiy{p£  -bp)/ 

'■  u,p-bu,p  1 

U.^y  +  UrpU. 

F,  =  -+-  u.pu^  , 

_J*.Ptr,  -b  Uypu,  -b  p 
\  u.(p£ -b  p) -b  i:,(p£  +  p)  / 
and  the  perturbation  to  the  source  term  is 


S=  n’py -b  2fl^,  (15) 

n*pz  - 

0  / 

The  relationship  between  the  second  order  linear  perturba¬ 
tion  A  and  U  is  exactly  the  same  as  the  relationship  between 
the  first  order  perturbation  A  and  U. 

After  considerable  algebra  the  quadratic  velocity  perturba¬ 
tions  are  found  to  be 

5y  =  -p-'  puT 

Sy  =  -p'’pu^  (16) 

S.  =  -p-’  Jul  , 
and  the  quadratic  pressure  perturbation  is 

-5(7-l)/>(^i  +  l*J  +  «‘!)-  (17) 

The  quadratic  perturbations  to  the  flux  vectors  are 
{  _  0  _ 

+  U£^+  P 

?,  =  gypUy  -b  (18) 

UiPU,  -b  Uypu, 

\  u,(p£-bp)  -b  Uyp  -b  U,p£  -b  57p  / 
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=  °  .  _ 

_  V-y,pU,  + 

(19) 

U„/>Uj  +  dyPti 

V^{pE  +  p)  +  Uyf  +  tl,p£  +  t 

0  \ 

g^/)u„  +  (20) 

e,^u.  +  g,pu,  +  f 

5,(p£+p)  +  «,p  +  u,pE  +  ulp  / 

and  the  quadratic  perturbation  to  5  is  zero. 

1.3  Steady-state  equations 


The  advantage  of  retaining  the  time-dependent  term  is  that  nu¬ 
merical  integration  of  this  equation  until  a  ‘steady-state'  solu¬ 
tion  IS  achieved  will  result  in  the  desired  harmonic  solution  This 
allows  one  to  employ  any  of  a  set  of  standard  numerical  time¬ 
marching  methods  with  the  usual  acceleration  techmques. 

1.5  Quadratic  equations 

The  quadratic  equations  come  from  taking  the  terms  which  are 
0{<’)  and  then  time-averaging. 

f.dA,  +  ~F^dAy  +  P.dA,  -  Jjj^SdV  =  QST  (27) 

The  quadratic  source  terms  (QST)  are 


The  steady-state  equations  are  obtained  by  substituting  the 
asymptotic  expansions  of  the  fluxes  £„£,,£,  and  the  source 
term  S  into  the  nonlinear  equations  of  motion,  and  retaining 
only  those  terms  which  are  0(1).  Dropping  the  superscript 
the  equations  are  identical  to  the  standard  nonlinear  steady- 
state  Euler  equations. 

F,  dA^  +  F^  dA^  +  F.dA,=  JII^S  dV.  (21) 

As  an  iterative  numerical  procedure  to  solve  these  equations 
it  is  convenient  to  add  back  in  an  unsteady  term 

i  JJJy  JL  -  IJL  ^ 

(22) 

and  integrate  in  time  until  a  steady-state  solution  is  achieved. 
This  can  be  accelerated  by  neglecting  time  accuracy  and  using 
techniques  such  as  local  time-stepping  and  multigrid. 

1.4  Linear  equations 

The  linear  perturbation  of  the  unsteady  integral  equation  of  mo¬ 
tion  is  obtained  bom  the  terms  which  are  0(c). 

I  IL  =  IllsdV.  (23) 

The  fact  that  this  equation  is  linear  and  all  of  the  coeffi- 
dearts  depend  on  U  and  geometric  terms,  which  do  not  vary 
in  time,  means  that  there  are  harmonic  solutions  of  the  form 
&(*)exp(  The  most  general  periodic  solution  can  be  ex¬ 

pressed  as  the  real  part  of  a  linear  superposition  of  such  solutions 
for  a  number  of  different  real  positive  frequendes  ur,  all  satisfying 
the  necessary  periodicity  constraints. 

U  =  71  ^51  exp{~iwt)^  .  (24) 

A  useful  idea  due  to  Ni  ( 1975)  is  to  formulate  the  harmonic 
equation  by  retaining  a  time  dependence  in  0.  Substituting 

U  =  U{x,t)exp{-iwt)  (25) 

gives  the  unsteady,  harmonic  equation 

I  +  =jj ljS  +  u^U)dV. 

(26) 


QST  =  -  T^dA^  -f  T^dAy  T  T.dA..  (28) 


Eq.  (27)  is  a  linear  equation  in  U;  the  quadratic  source  terms 
are  the  inhomogeneity  that  drives  this  second  order  flow  pertur¬ 
bation. 

Again  using  Ni’s  idea  of  adding  a  pseudo- unsteady  term  to 
allow  the  use  of  numerical  time  integration  schemes,  one  obtains 

I,  IfL  ^  ^  *  IL  -  UP 

(29) 

In  this  form,  it  is  equivalent  to  the  linear  harmonic  equation  with 
zero  frequency  and  the  addition  of  the  quadratic  source  terms. 
Thus  it  can  be  solved  numerically  with  only  trivial  changes  to 
the  numerical  procedure. 

An  alternative  approach  is  to  define  a  time-averaged  state 

U  =  U  +  €’‘U  =  (/W-l-ct/Ol-ht’t/lri,  (30) 


and  recombine  the  steady-state  equation  with  the  second  order 
equation  to  form 

(31) 

This  effectively  introduces  into  the  nonlinear  steady-state  equa¬ 
tion  the  time-averaged  effect  of  the  linear  unsteadiness  through 
the  simple  inclusion  of  the  quadratic  source  terms. 

A  final  issue  is  the  evaluation  of  the  quadratic  source  terms 
for  linear  unsteady  solution  U  which  is  represented  as  the  real 
part  of  a  sum  of  harmonic  solutions  of  different  bequendes. 

U  =  11  ^51  &«exp(-twt)j  (32) 

It  turns  out  that  this  is  very  simple.  For  an  arbitrary  scalar 
function  A(V), 


^  -  2dUjd[/M^ 
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Thus,  the  quadratic  source  terms  can  be  expressed  as  a  sum 
of  the  quadratic  source  terms  for  each  distinct  harmonic  fre¬ 
quency 

QST  =  £QST,,  (34) 

and  each  of  the  harmonic  components  is  easily  evaluated  know¬ 
ing  U. 

2  Boundary  conditions 


2.2  Linear  unsteady  flow 

The  boundary  conditions  on  blade  surfaces  and  endwalls  are 
again  that  there  is  no  mass  flux  across  the  boundary.  When 
this  condition  is  linearized  it  produces  the  requirement  that  the 
normal  component  of  the  velocity  perturbation  is  zero. 

The  periodic  boundary  condition  is  slightly  mote  compli¬ 
cated  than  for  the  steady  state  flow  due  to  the  inclusion  of  the 
interblade  phase  angle  4>- 


2.1  Steady  flow 

The  boundary  conditions  on  blade  surfaces  and  endwalls  is  sim¬ 
ply  that  there  is  no  mass  flux  across  the  boundary,  or  equiva¬ 
lently  that  the  normal  component  of  the  velocity  is  zero. 

If  the  calculation  is  performed  for  a  single  blade-blade  pas¬ 
sage  in  a  blade  row  with  multiple  blades,  then  the  periodic 
boundary  condition  is  that 

u,(x,g,T)  =  u,{xj  +  Ae),  (35) 

where  Afl  is  the  angular  pitch  and  the  ‘primitive’  variable  vector 
Up  in  cylindrical  coordinates  is  deflned  as 


At  inflow  and  outflow  boundaries,  steady-state  non-reflecting 
boundary  conditions  can  be  constructed  (Saxer,  1991).  In  these, 
one  specifies  at  the  inflow  circumferentially  averaged  values  of 
total  temperature,  total  pressure . surd  flow  angles  as  functions 
of  radius.  At  the  outflow  dicumferentially  averaged  values  of 
pressure  are  specified,  consistent  with  radial  equilibrium.  Cir¬ 
cumferential  harmonics  of  Up  are  not  set  to  zero  but  are  allowed 
to  evolve  to  values  which  are  consistent  with  the  existence  of  an 
infinite  annular  duct  upstream  or  downstream  of  the  computa¬ 
tional  domain. 

At  stator/rotor  interfaces  a  very  similar  treatment  is  used, 
but  the  drcumferentially  averaged  values  of  Up  are  determined 
instead  by  the  constraint  that  the  total  flux  of  mass,  momentum 
and  energy  across  the  interface  must  match.  Defining  the  axial 
flux  in  cylindrical  coordinates  as 

'  pti,  ^ 

l^l  +  p 

Fm$  =  »  (37) 

k  u.(/»£+p)  , 

with  all  quantities  being  in  the  absolute  frame,  then  the  match¬ 
ing  condition  between  blade  rows  1  and  2  is  that 

■^•♦row  1  ~  ^•♦tow  2>  (38) 


with  the  notation  denoting  a  circumferential  average  over  one 
blade  passage.  This  flux-conservation  treatment  was  first  used 
by  Denton  and  csm  cause  a  small  entropy  rise  due  to  the  implicit 
mixing  process  (Denton,  1985). 


Up(x,g,r)  =  exp(i^)U,(*,^-f-AS,r)  (39) 

This  is  still  a  very  much  simpler  boundary  condition  to  imple¬ 
ment  numerically  than  the  lagged  periodic  boundary  condition 
in  single  stage  nonlinear  stator/rotor  calculations  (Erdos,  1977, 
Giles,  1988).  Also,  linearity  allows  one  to  superimpose  two  so¬ 
lutions  with  different  interblade  phase  angles  due  to  differing 
number  of  blades  in  upstream  and  downstream  blade  rows. 

The  most  interesting  boundary  conditions  are  the  inflow  and 
outflow  conditions  used  to  model  unsteady  stator/rotor  interac¬ 
tion  in  a  linear  sense.  Consider  two  neighboring  blade  rows, 
labeled  1  and  2,  rotating  with  speeds  and  fij.  The  key  idea 
is  that  a  stationary  circumferential  variation  in  one  blade  row 
becomes  a  rotating  unsteady  perturbation  in  the  other  blade 
tow  due  to  the  relative  rotation  of  the  two  tows.  The  steady- 
state  primitive  flow  variables  in  cylindrical  coordinates  in  blade 
row  I  can  be  expressed  as  the  sum  of  am  axisymmetric  part  plus 
circumferential  ha  ionic  terms. 

-lfp{x,r)  +  Tli^Y.  | 

The  angle  in  neighboring  blade  row  2  is  related  through 

+  {lit  =  9t  +  D,t.  (41) 


Hence, 

2iemBi  _  2irm(Bi  -  (fli-n,)t) 

Afii  ABi  ■ 

and  so  the  stationary  m'*  mode  in  blade  row  1  is  equivalent  to 
an  unsteady  mode  in  blade  row  2  with  the  frequency 

2xm(n,-0,)  . 


The  harmonic  boundary  terms  that  go  with -this  mode  are 

^p(*,«i.»-)  =  exp^— I/,„(z,r)  (44) 

where  the  inter-blade  phase  angle  is 


*  =  - 


2TTn  Afla 


A  slight  complication  is  the  desire  to  keep  all  frequencies  w 
positive.  If  the  above  description  gives  a  value  for  ui]  which  it 
negative  then  the  problem  is  avoided  by  noting  that  a  harmonic 
solution  with  frequency  u  and  inter-blade  phase  angle  4>  the 
complex  conjugate  of  wother  solution  with  frequency  -w  and 
interblade  phase  angle 
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2.3  Quadratic  equations 

The  sobd  wall  and  periodic  boundary  conditions  are  exactly  the 
same  as  for  the  steady  flow  analysis.  The  inflow  and  outflow 
boundary  conditions  are  also  similar  in  that  the  citcumferen* 
tial  harmonics  are  again  specified  by  steady-state  non-reflecting 
boundary  condition  theory.  At  stator/rotor  interfaces,  the  cir¬ 
cumferentially  averaged  values  of  the  second  order  flow  pertur¬ 
bations  are  set  by  the  requirement  that  the  average  fluxes  (aver¬ 
aged  both  in  time  and  circumferentially)  must  match.  Substitut¬ 
ing  the  asymptotic  expansion  of  the  flux  and  then  time-averaging 
gives  the  following  equation. 

- = - — — 0  - ^ - - — 0 

«»  +  1  =  +  ^*#)row  2- 

Canceling  the  0(1)  terms  which  match  in  the  steady-state  equa¬ 
tions,  yields 


F,0  +  ■^*srow  I  —  ■^*#  +  F**iov/  2-  (^~) 

Alternatively,  if  one  adopts  the  other  approach  to  the  second 
order  analysis  by  adding  U  to  U  to  form  the  time-average  U, 
then  one  obtains  the  matching  condition 

- ^~0  - - - 

7x0  +  <*^*»row  1  =  7,0  +  t^F,0 

row  2*  (48) 

In  this  form,  the  matching  condition  is  similar  to  a  mathematical 
constraint  derived  by  Adamczyk  to  ensure  the  correct  transfer 
>f  time-averaged,  drcumferentially  averaged  vorticity  across  the 
interface  (Adamczyk,  1991). 

3  Possible  Applications 

The  purpose  of  this  section  is  to  discuss  the  variety  of  different 
phenomena  which  can  be  analyzed  using  the  approach  outlined 
in  this  paper. 

3.1  Wake  and  hot  streak  migration 

In  a  compressor,  rotor  wakes  migrate  towards  the  pressure  sur¬ 
face  of  the  stators  downstream  (Kerrebrock,  1970).  Similarly, 
in  a  turbine  first  stage,  hot  streaks  coming  from  spatial  nonuni- 
fonnities  in  the  combustor  migrate  towards  the  pressure  surface 
of  the  rotor  leading  to  elevated  heat  transfer  (Butler,  1986).  In 
both  cases  the  mechanism  is  essentially  the  same.  The  unsteady 
velocity  of  the  wake  or  hot  streak  relative  to  the  mean  flow  drives 
the  fluid  across  the  passage  towards  the  pressure  surface.  The 
time-averaged  effect  is  quadratic  in  the  level  of  unsteadiness. 
For  example,  the  elevated  surface  temperature  in  the  turbine 
is  proportional  to  both  the  unsteady  velocity  and  the  unsteady 
perturbation  in  temperature. 

Such  phenomena  will  be  captured  by  the  second  order  anal¬ 
ysis  using  the  quadratic  source  terms.  In  the  case  of  the  hot 
streaks,  the  quadratic  term  in  the  energy  flux  which  reflects  the 
nigration  effect  is  of  the  form  UyT. 


3.2  Time-averaged  total  temperature  varia¬ 
tions 

In  fans  and  compressors  an  unsteady  rotor/stalor  interaction 
can  lead  to  a  circumferential  variation  in  time- averaged  total 
temperature  in  the  absolute  frame.  This  spatial  variation  can 
cause  significant  complications  in  the  experimental  determina¬ 
tion  of  the  efficiency  of  the  rotor.  This  phenomenon  has  been 
documented  and  explained  fully  in  a  paper  by  Epstein  e(  al 
(1989).  It  is  due  to  the  potential  field  of  the  stator  causing  an 
osdllation  in  time  of  the  circulation  around  the  rotor.  This  pro¬ 
duces  a  shed  vortex  sheet  which  in  turn,  in  the  absolute  frame  of 
the  stator,  gives  a  spatial  variation  in  total  temperature  which 
can  be  explained  through  the  time-averaging  of  Crocco’s  theo¬ 
rem.  The  phenomenon  is  linear,  with  the  amplitude  of  the  total 
temperature  nonuniformity  being  proportional  to  the  ampbtude 
of  the  oscillation  in  the  rotor  blade  circulation. 

This  phenomenon  will  be  captured  by  the  steady-state  and 
linear  perturbation  analysis  in  a  three-step  process.  This  first 
step  is  a  steady-state  combined  rotor/stator  computation.  The 
second  step  is  the  computation  of  the  linear  harmonic  unsteadi¬ 
ness  in  the  rotor  due  to  the  first  (and  possibly  second)  circum¬ 
ferential  spatial  harmonics  in  the  stator  inflow  plane.  The  un¬ 
steadiness  in  the  rotor  frame  of  reference  produces  downstream 
propagating  acoustic  and  vortical  disturbances.  A  circumferen¬ 
tial  Fourier  decomposition  shows  that  some  of  the  modes  are 
stationary  in  the  absolute  frame  of  reference.  These  modes  are 
then  input  as  boundary  conditions  for  a  steady  linear  perturba¬ 
tion  analysis  in  the  stator  which  will  produce  the  time-averaged 
spatial  variation  of  total  temperature.  In  this  application  the 
second  order  quadratic  effects  do  not  play  a  significant  role. 

3.3  Losses  due  to  unsteadiness 

One  of  the  most  important  parts  of  stator/rotor  interaction  is 
the  potential  interaction,  the  effect  of  the  pressure  field  about 
one  blade  row  on  a  neighboring  blade  row  in  relative  motion. 
The  potential  interaction  generates  acoustic  waves,  and  with 
the  exception  of  the..situation  discussed  in  the  previous  exam¬ 
ple,  most  of  the  waves  are  unsteady  in  both  frames  of  reference. 
Associated  with  these  waves  is  a  certain  level  of  acoustic  energy, 
and  the  question  which  then  arises  is  whether  this  acoustic  en¬ 
ergy  is  eventually  dissipated  (converted  into  thermal  energy  with 
an  increase  in  entropy  and  a  decrease  in  stagnation  pressure). 

This  topic  is  addressed  in  a  paper  by  Fritsch  (1992)  in  which 
it  is  argued  that  most  of  this  energy  will  be  lost.  The  physical 
m  .ichetnism  is  the  unsteady  pressure  field  generating  an  unsteady 
Stokes’  layer  at  the  bottom  of  the  boundary  layer,  leading  to  in¬ 
creased  dissipation  in  the  boundary  layer.  By  contrast,  an  order- 
of-magnitude  analysis  comes  to  the  conclusion  that  in  numerical 
computations  the  energy  is  lost  in  the  fieestream  portion  of  the 
flow  due  to  numerical  smoothing.  A  supposition  (possibly  plau¬ 
sible  or  possibly  just  hopeful)  is  that  a  numerical  computation 
may  still  correctly  calculate  the  overall  additional  lots  provided 
it  correctly  predicts  the  rate  of  generation  of  the  acoustic  energy 
due  to  the  potential  interactions. 

If  this  supposition  is  true,  and  a  nonlinear,  time-marching 
flow  calculation  including  numerical  smoothing  can  indeed  pre¬ 
dict  the  increased  losses  due  to  this  form  of  unsteadiness,  then 
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the  full  steady /linear/quadratic  analysis  in  this  paper  should  be 
equally  successful.  This  is  a  quadratic  effect  since  the  acoustic 
energy  level  is  proportional  to  the  square  of  the  unsteady  flow 
quantities. 

If  the  supposition  is  not  correct  then  it  will  be  necessary  to 
extend  the  theory  in  this  paper  to  an  asymptotic  analysis  of  the 
Navier-Stokes  equations.  It  will  sdso  be  necessary  to  use  very 
much  finer  computational  grids  to  eliminate  the  errors  due  to 
the  numerical  smoothing.  Of  course  this  will  greatly  increase 
the  computational  cost  of  such  calculations. 

4  Concluding  Remarks 

This  paper  has  presented  an  asymptotic  approach  toward  the 
calculation  of  unsteady  flows  in  multistage  turbomachinery.  The 
analysis  proceeds  in  three  phases.  First  one  performs  a  steady- 
state  multistage  analysis  in  which  the  flow  is  assumed  to  be 
steady  in  each  blade  row  in  a  frame  of  refcrehce  rotating  with  the 
blade  row.  The  stator/rotor  interface  matching  conditions  en¬ 
sure  overall  conservation  of  mSLSS  momentum  and  energy  across 
the  interface. 

In  the  second  phase  the  unsteadiness  in  each  blade  row  is 
treated  as  a  perturbation  to  the  steady  flow  field.  Assuming  the 
perturbation  is  small  it  is  represented  by  the  linearized  unsteady 
flow  equations.  The  linear  stator/rotor  interaction  is  calculated 
by  recognizing  that  a  steady  circumferential  variation  in  a  stator 
corresponds  to  a  rotating  unsteady  disturbance  in  a  rotor,  and 
ince  versa.  Each  circumferential  harmonic  produces  a  response 
of  a  different  frequency.  Due  to  linearity  these  can  be  com¬ 
puted  independently  and  superimposed  to  form  the  composite 
unsteady  flow. 

The  third  phase  calculates  the  change  to  the  mean  flow  due 
to  second  order  consequences  of  the  unsteady  flow  field.  This  is 
achieved  through  the  evaluation  of  quadratic  source  terms  which 
sre  similar  in  origin  to  the  Reynolds  stress  terms  in  the  Reynolds 
averaged  Navier-Stokes  equations. 

Through  this  process  one  can  calculate  the  key  features  of 
the  unsteady  flow  in  turbomachinery,  the  unsteady  forces  and 
moments  on  blades  due  to  interactions  with  neighboring  rows, 
as  well  as  the  average  effects  on  mass  flow  and  efficiency.  The 
efficiency  changes  come  about  as  the  result  of  wasted  kinetic  en¬ 
ergy  associated  with  the  eventual  dissipation  of  acoustic  energy 
generated  by  unsteady  stator/rotor  interaction.  A  more  com¬ 
plete  treatment  of  efficiency  prediction  would  require  the  appli¬ 
cation  of  the  entire  asymptotic  approach  to  the  Navier-Stokes 
equations.  This  is  an  algebraically  complicated  but  otherwise 
straightforward  extension  which  will  be  tackled  in  the  future. 

This  paper  does  not  present  any  computational  results,  but 
clearly  this  conceptual  approach  needs  to  be  thoroughly  vali¬ 
dated  by  comparison  to  nonlinear  time-marching  methods  and 
experimental  data.  This  will  be  addressed  in  future  work,  along 
with  the  task  of  accelerating  the  numerical  computations  using 
techniques  such  as  multigrid.  The  two  reasons  at  present  to 
beUeve  that  there  is  some  merit  to  this  idea  are  that  Adam- 
czyk  has  already  achieved  significant  improvements  using  his 
passage-averaged  flow  equations,  and  the  alternative  approach 
of  nonlinear  unsteady  multistage  calculations  requires  excessive 
computation  times. 
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TASK  III:  UNSTEADY  PHENOMENA  AND  FLOWFIELD  INSTABILITIES 

IN  MULTISTAGE  AXIAL  COMPRESSORS 

(Investigators:  E.M.  Greitzer,  C.S.  Tan,  P.P.  Silkowski,  J.P.  Longley, 

N.A.  Cumpsty,  IJ.  Day,  D.C  Wisler,  H.W.  Shin) 

Executive  Summary 

This  progress  report  describes  a  combined  experinaental  and  theoretical  program  of 
research  on  unsteady  flow  and  instability  in  multistage  axial  compressors.  The  objectives  of  the 
program  were  to  define  the  critical  fluid  dynamic  phenomena  associated  with  the  onset  of 
instability  in  modem,  multistage  machines  and  assess  those  elements  that  should  be  incorporated 
in  improved  models  of  this  process.  The  specific  research  goals  included  elucidation  of  newly 
discovered  types  of  propagating  disturbances,  having  widely  disparate  wavelengths  and  waves 
speeds,  which  occur  prior  to  rotating  stall  onset,  and  defmition  of  the  extent  to  which  current 
unsteady  flow  models  capture  the  relevant  fluid  processes  of  this  complex  dynamical  system.  The 
work  is  being  carried  out  in  collaboration  with  the  General  Electric  Aircraft  Engine  Company,  with 
the  experiments  actually  being  ctMiducted  at  the  GE  Aerodynamics  Research  Labewatory,  on  a 
facility  representative  of  advanced  aeroengine  compressor  geometries. 

Two  main  activities  are  discussed  in  the  report.  The  first  is  unsteady  measurements  of  the 
blade  passage  flo' v  field  in  the  first  rotor  of  a  multistage  compressor,  under  conditions  at  which 
stall  occurred  naturally  as  well  as  under  conditions  in  which  the  flow  was  artificially  stabilized,  so 
that  the  behavior  at,  and  even  past,  the  nominal  stall  point  could  be  addressed.  This  gives  new 
information  ctxiceming  the  intmiction  and  role  of  the  different  wavelength  unsteady  disturbances. 
The  seemd  is  a  series  of  diagnostic  experiments  with  a  forced  unsteady  inlet  distortion,  which 
allo’Mvd  detailed  assessment  of  the  ctmeepts  that  underlie  the  existing  wave  dynamic  models  of 
compressor  response  to  inlet  non-uniformity. 
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1.0  IntrodutUon 

Recent  studies  of  rotating  stall  inception  have  led  to  the  definition  of  two  different  paths 
from  the  essentially  axisymmetric  flow  in  the  pre-stall  state  to  fully  developed  rotating  stall.  The 
first  involves  the  growth  of  a  small  amplitude,  long  circumferential  length  scale  (of  the  order  of  the 
machine  radius)  traveling  wave  into  a  fully  developed  rotating  stall  cell.  The  initial  propagation 
speed  of  the  traveling  wave  is  close  to  that  of  the  developed  stall  cell,  and  small  amplitude  traveling 
waves  can  exist  for  a  long  time  (several  tens  of  rotor  revolutions)  before  any  operational  changes 
in  machine  behavior  occur.  This  process  is  well  captured  by  the  models  of  unsteady  flow  in 
multistage  machines. 

Another,  quite  different,  process  has  also  been  observed,  however,  which  involves  the 
circumferential  growth  of  a  small  circumferential  extent  disturbance  (a  few  blade  pitches)  into  a 
mature  rotating  stall  cell  in  roughly  five  to  ten  rotOT  revolutions.  These  small  extent,  or  short 
length  scale,  phenomena  will  be  referred  to  here  as  spikes,  primarily  because  their  temporal  history 
appears  considerably  more  “spiky”  than  the  long  length  scale  disturbances.  An  important  factor  in 
the  design  of  the  experiments  to  be  discussed  was  the  desire  to  focus  on  the  different  time  and 
length  scales  associated  with  the  spikes  and  the  long  wave  disturbances  (referred  to  as  LWD’s).  In 
particular,  the  axial  locations  of  probes  upstream  of  the  rotor  were  intended  to  spatially  filter  out 
short  length  scale  disturbances  and  make  the  long  wavelength  disturbances  clearer.  Also,  several 
rig  modifications  were  used  as  tools  to  focus  on  specific  time  and  length  scales,  including 
restaggering  the  rear  three  stages,  asymmetric  tip  clearance,  and  distention  screens. 

The  basic  fluid  dynamic  questions  that  were  initially  addressed  in  this  study  were: 

1 )  What  are  the  mechanisms  of  rotating  stall  inception? 

2)  Is  there  a  repeatable  pre-stall  event,  and  if  so,  is  it  a  long  wavelength  distinbance  (LWD),  a 
short  wavelength  disturbance  (spike),  or  some  combination  of  the  two? 

3)  Can  one  associate  tiie  type  of  pre-stall  events  with  con^ressor  design  parameters? 

4)  At  the  blade  passage  length  scale  how  do  the  various  fluid  dynamic  structures  (tip 


121 


voitcx/leakage  jet,  blade  wakes,  and  cndwall  boundary  layers)  behave  w*iih  changing  operating 
point  and  what  role  do  these  structures  play  in  the  stall  inception  process? 

5)  How  do  inlet  distcwtion  and  changes  in  rip  clearance  affect  the  stall  inception  process? 

6)  If  we  are  able  to  artificially  stabilize  the  first  stall  flow,  what  is  the  axisyrnrnctric  machine  and 
local  flow  field  behavior  in  the  first  stage  at  these  highly  loaded  conditions,  i.e.  what,  in  detail, 
is  the  flow  field  associated  with  "compressor  stall"? 

To  address  these  questions,  an  investigation  of  a  multi-stage  compressor  flowTicld  has 
been  conducted  involving  several  different  levels  of  experimentation.  Detailed  examination  of  the 
first  stage  rotor  has  been  carried  out  during  a  scries  of  experiments  in  which  the  flow  was 
artificially  stabilized  by  rcstaggering  the  rear  three  stages.  This  restaggering  created  a  situation  in 
which  the  behavior  of  the  rear  three  stages  suppressed  the  growth  of  long  length  scale  flow 
disturbances.  This  enabled  the  behavior  of  the  first  rotor,  at  flows  below  that  where  the  rotor 
would  experience  stall  in  a  natural  enviixMiment,  to  be  seen.  As  will  be  described,  the  experiments 
show  that  although  the  stall  inception  process  may  begin  as  a  localized  event,  the  evolution  into 
rotating  stall  is  governed  by  the  environment  defined  by  the  coupling  of  the  various  compression 
system  conqxments.  Also,  the  rotor  tip  leakage  jc*  appears  to  be  a  key  feature  in  dte  local  stall 
inception  {ffocess. 

As  a  second  facet  of  the  program,  diagiwstic  experiments  have  been  carried  out  to  examine 
the  effect  of  unsteady  inlet  distoiticms  cm  instability  onset,  i.e.  to  address  the  issue  of  the  receptivity 
of  the  compressor,  considered  as  a  multi-dimcnsional  dynamic  system,  to  inlet  distuitmnces.  This 
set  of  experiments  also  allows  assessment  of  wave  models  of  unsteady  compressor  response  that 
have  been  developed. 

The  first  experiments  showed  clearly,  perhaps  for  the  first  time,  that  the  path  to  stall  can 
involve  the  interaction  of  disturbances  of  both  long  and  short  wavelength  pmuibations. 
Undmtanding  the  connection  between  these  types  of  disturbances  is  a  critical  new  conceptual  step 
in  attacking  this  impe^nant  problem. 
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2.0  EApflcimental  facilities 

The  investigation  included  theoretical  and  experimental  phases  and  was  a  collaborative 
university-industry  effort,  involving  groups  at  MIT  and  Genwal  Electric  Aircraft  Engines.  There 
was  also  participation  from  the  Whittle  Lab,  Camteidge,  England.  The  experiments  were  carried 
out  at  the  General  Electric  Aerodynamics  Research  Laboratory,  under  the  aegis  of  Dr.  D.C. 

Wisler.  The  facility  used  was  a  large  scale,  low  speed  research  compressor  that  presents  a  model 
of  many  aspects  of  the  actual  high  speed  compressor  flow  field  (Fig.  1).  The  large  size  minimizes 
probe  blockage  and  allows  for  easier  access  and  instrumentation  whilst  g»\  ing  correct  Reynolds 
number  at  low  operating  speed.  A  four-stage  compressor  was  used  with  the  NASAAjE  Energy 
Efficient  Engine  blading.  Further  information  on  this  facility  and  the  blading  can  be  obtained  from 
Ref[l]. 

In  addition  to  standard  steady  state  instrumentation,  high  response  pressure  transducers  and 
hot  wires  were  used.  A  radially  traversable  slant  wire,  to  give  three  components  of  velocity,  and  a 
total  pressure  probe  were  also  used  at  rotor  inlet  and  rotor  exit.  The  bulk  of  the  experiments  were 
carried  out  with  fixed  arrays  of  40  probes,  consisting  of  20  hot  wires  and  20  pressure  transducers. 
At  inlet  guide  vane  exit,  there  were  ten  reference  hot  wires  at  midspan  and  ten  corresponding 
pressure  transducers  in  the  casing.  The  ten  refnence  locations  were  circumferentially  distributed  at 
eight  equally  spaced  locations  to  facilitate  spatial  Fourier  analysis.  Two  closely  spaced  sensors 
were  also  used  to  allow  cross-correlations  over  various  distances,  to  combat  spatial  aliasing.  Most 
of  the  unsteady  data  was  sampled  at  10  kHz  and  low  pass  filtered  at  3  kHz,  with  several  cases 
sampled  at  1  kHz  and  low  pass  filtered  at  300  Hz. 

Ccnnplementing  the  reference  probes  were  ten  additional  hot  wires  and  ten  casing  pressure 
transducers  with  circumferential  distribution  i^ntical  to  the  reference  probes  refored  to  above. 

The  reference  probes  were  fixed  at  their  locations  for  all  runs  as  a  control,  but  the  test  probes  were 
moved  to  various  axial-radial  positions.  The  test  probes  were  used  at  four  axial  locations:  one  half 
of  a  radius  and  one  quarter  of  a  radius  upstream  of  rotw  inlet,  one  at  rotor  inlet,  and  one  at  rotor 
exit,  at  locations  roughly  0.2  chord  upstream  and  downstream  of  the  rotor  respectively.  At  R/2  or 
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R/4  upstream,  the  hot  wires  were  at  midspan,  but  at  rotor  inlet  and  rotor  exit,  the  hot  wires  were 
used  at  20%  and  80%  immersion  (tip  and  hub),  as  well  as  at  other  immersions,  as  shown  in  Fig. 

2.  To  focus  on  small  length  scale  events,  experiments  were  also  performed  with  eight  probes 
equally  spaced  over  three  blade  pitches. 

The  baseline  build  overall  compressor  characteristic  (pressure  rise  versus  axial  velocity 
parameter  at  constant  rotor  speed)  is  displayed  in  Fig.  3.  Figure  4  is  the  performance  of  the  first 
stage  only.  Several  points  are  labeled  for  future  reference.  Point  A  is  design  and  point  E  is  just 
before  rotating  stall  occurs. 

For  the  restagger,  the  rear  three  rotors  were  closed  by  10°.  Figure  5  she  vs  the  result  of  the 
restagger  on  the  overall  pressure  rise  characteristic.  The  stalling  flow  coefficient  was  changed 
from  0.345  to  0.308.  Thus,  the  purpose  of  the  restagger,  which  was  to  extend  the  stable  operating 
region  of  the  first  stage  and  provide  a  clearer  view  of  the  spikes  by  suppressing  the  long  length 
scale  disturbances,  has  been  achieved.  There  is  an  alnnjpt  change  in  slope  of  the  overall 
mismatched  characteristic  at  the  stalling  ^  of  the  baseline  build. 

Figure  6  displays  the  first  stage  characteristic  for  the  baseline  and  restaggered  builds  with 
several  operating  points  labeled  for  future  reference.  In  the  range  of  flows  comnaon  to  the  two 
builds,  the  first  stage  pressure  rise  characteristic  is  the  same  for  both.  The  extended  operating 
region  of  the  first  stage  characteristic  (achieved  with  the  restagger),  however,  is  positively  sloped. 
Stability  analysis  predicts  that  a  positively  slc^ied  compressor  stage  would  exhibit  flow  instability; 
here  the  rear  three  stages  are  negatively  slcq)ed  and  the  overall  compressor  has  a  negative  slope  and 
hence  is  stable. 

A  seccHid  variation  from  the  baseline  build  was  to  open  the  first  stage  rotor  clearance  from 
gap/span  =  1.3%  to  3.1%.  This  was  done  both  symmetrically  as  well  as  by  opening  the  clearance 
over  a  120°  section  of  the  annulus. 
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3.0  Results 

3.1  Steady  State  Blade  Row  Flow  Field  Data 

The  steady  state  blade  row  data  are  presented  in  two  ways:  radial  profiles,  i.e.  mass 
averages  over  one  circumferential  pitch  at  a  given  radius,  and  contour  or  vectOT  format  over 
roughly  l.S  blade  pitches. 

Figure  7  shows  the  radial  velocity  distribution  at  inlet  to  the  first  rotor  at  various  operating 
points.  This  should  be  close  to  zero  if  the  fiow  is  unstalled.  Although  the  flow  angles  in  the  tip 
region  for  the  extended  operating  region  of  the  restaggered  build  were  out  of  calibration  range  for 
the  slant  wire,  useful  information  can  still  be  obtained  from  the  data.  In  particular,  the  plot  shows 
that  the  tip  flow  field  does  not  deviate  frcnn  the  expected  direction  until  operating  point  J  of  the 
restagger  build,  which  is  the  baseline  stalling  flow  coefficient.  As  the  flow  coefficient  decreases 
further,  an  increasingly  large  portion  of  the  tip  flow  field  is  out  of  range  for  the  slant  wire. 

Figure  8  is  a  display  of  the  slant  wire  data  at  rottx-  exit  at  the  stable  operating  point  closest 
to  stall  in  the  baseline  build.  The  data  has  been  converted  to  the  relative  (rotcn-)  fiame  and  is 
viewed  at  a  yaw  angle  of  55°  fiom  axial,  i.e.  an  angle  aligned  with  the  bulk  through  flow  direction, 
so  the  figure  represents  the  cross-plane  or  secondary  flow.  The  scale  vector  at  the  top  shows  the 
direction  of  rotor  rotation  and  the  magnitude  of  blade  tip  q>eed,  Ut^.  Blade  wakes,  which  have 
strong  radial  transport,  and  the  tip  vortex/leakage  jet  flow  is  clearly  visible. 

Hgures  9  and  10  show  the  throughflow  velocity  at  different  flow  coefficimts  with  the 
contour  of  velocityAJtip  <  0.45  shaded.  A  change  in  the  tip  region  between  points  D  and  K  is 
seen;  diere  appears  to  be  an  almost  axisymmetiic  region  of  low  velocity  near  the  tip  at  point  K 
Figure  1 1  shows  the  axial  velocity  distribution  at  50%  immersion  (mid-span)  at  rotor  exit  for 
roughly  1.5  pitch  at  three  diiferent  flow  coefficients,  and  it  is  seen  that  the  rotor  wakes  do  not  grow 
much  circumferentially  in  0  as  0  decreases.  Figure  12  displays  axial  velocity  profiles  fen-  two 
operating  points  at  rotor  exit  for  the  baseline  and  for  a  larger  first  stage  tip  clearance.  An  increase 
of  axial  flow  deficit  can  be  seen  near  the  tip  with  the  larger  clearance. 

To  summarize  the  steady  state  data,  as  the  flow  ((|>)  decreases  the  tip  vortex/jet  structure 
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and  the  rotor  wakes  do  not  change  qualitatively,  but  close  to  stall  the  rotor  leakage  jet  becomes 
larger.  Furthermore,  for  the  restaggered  case,  starting  at  the  baseline  stall  flow  coefficient,  a  region 
of  appreciable  radial  velocity  at  rotor  inlet  grows  in  radial  extent  with  decreasing  flow  coefficient. 
Finally,  although  the  data  is  not  shown  here,  the  larger  tip  clearance  build  had  a  larger  axial  velocity 
deficit  in  the  tip  region,  than  the  baseline  build. 

4.0  TimfirResotved  Data 
4.1  Baseline  Build 

Figure  13  shows  time  traces  for  eight  equally-spaced  casing  pressure  transducers  at  rotor 
inlet.  Probe  numbers  increase  in  the  direction  of  increasing  6  in  the  rotor  rotation  direction,  and 
time  is  in  units  of  rotor  revolutions.  'Hie  figure  displays  a  situation  in  which  the  compressor  was 
set  at  a  stable  operating  point  close  to  stall  and  the  throtde  slowly  closed  until  the  machine  went 
into  stall.  A  spike  is  first  visible  at  t  ~  76.5  on  pressure  transducer  number  5.  The  growth  of  the 
spike,  in  amplitude  and  circumferential  extent,  is  seen  over  the  next  few  rotor  revolutions,  resulting 
in  a  developed  stall  cell  by  t  ~  82.  There  are  no  apparent  propagating  disturbances  prior  to  time  ~ 
76.5. 

Guidelines  have  been  drawn  in  Fig.  14,  with  the  slope  of  these  lines  representing  the  speed 
of  propagation  of  the  disturbance  around  the  annulus.  The  change  in  slc^  of  these  lines  indicates 
the  change  in  the  speed  of  propagation  of  the  disturbance  as  it  decelerates  from  70%  at  inception  to 
45%  as  a  mature  stall  cell.  No  evidence  of  any  long  wavelength  disturbances  above  the 
perturbation  levels  of  (perturbation  velocity)/U  =  0.5%,  p’/pU^  =  1%  could  be  found.  The  data 
analysis  techniques  used  confirmed  the  spike’s  70%  initial  speed  and  subsequent  deceleration  to 
45%  with  growth  in  magnitude  and  circumferential  extent 

An  example  of  the  data  analysis  technique  used  is  shown  in  Fig.  15.  At  each  discrete  time, 
data  fiom  eight  equally  spaced  probes  were  Fourier  analyzed  and  the  resulting  magnitudes  and 
phases  for  the  various  harmonics  displayed.  The  size  of  the  symbols  are  scaled  by  the  amplitude 
of  the  disturbance.  No  propagating  disturbances  are  visible  prior  to  time  =  77,  at  which  point  the 
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disturbance  grows  rapidly  to  final  size  within  roughly  five  rotor  revolutions.  The  continuous 
change  in  slope  of  the  phase  of  the  first  Fourier  harmonic  indicates  the  deceleration  of  the  spike. 
(The  regular  pattern  seen  for  times  prior  to  t  ~  77  in  harmonic  one  is  attributed  to  a  slight  rotor 
geometry  asymmetry.) 

Some  experiments  were  conducted  with  die  eight  equally  spaced  test  wires  alternating 
between  hub  and  tip,  80%  and  20%  immersion.  Figure  16  displays  the  results  from  one  of  these 
tests  conducted  at  rotor  inlet,  which  shows  the  three-dimensional  nature  of  the  spike.  Between  t  -7 
7  and  t  ~  78,  the  spike  is  seen  as  an  axial  deficit  traveling  in  the  tip  region,  with  no  disturbance 
visible  in  the  hub  wires.  Only  after  the  tip  deficit  has  grown  appreciably  is  any  disturbance  visible 
at  the  hub  at  t  ~  78.5.  The  hub  disturbance  is  first  seen  as  a  flow  increase  (partially  compensating 
for  the  tip  blockage),  but  as  the  spike  grows  into  a  stall  cell,  a  full  span  deficit  eventually  occurs. 

Time-resolved  measurements  thus  show  the  baseline  build  stalled  via  radial  and 
circumferential  growth  of  an  initially  small  length  scale  disturbance.  The  growth  occurred  in 
roughly  five  rotw  revolutions  and  no  modal  signals  were  found.  The  disturbance  initial  speed  of 
70%  decreased  to  45%  as  it  grew  into  a  mature  stall  cell. 

5.0  Restagyered  Build 

As  the  throttle  was  closed  for  the  mismatched  build,  the  following  sequence  of  events 
occutred.  When  the  flow  coefficient  decreased  to  the  level  at  which  the  baseline  encountered 
rotating  stall,  spikes  would  appear  and  disappear.  As  the  flow  was  decreased  further,  to  J,  mac 
spikes  would  appear  per  unit  time  with  the  spike  width  increasing  from  roughly  2-3  blade  pitches 
to  3-4  blade  pitches. 

As  tile  flow  coefficient  was  decreased  from  J,  the  spike  disturbances  became  more  regular, 
they  did  not  decay  after  a  short  appearance,  but  persisted.  Time  traces,  FFTs,  and  ctMrrelations  all 
showed  this.  The  disturbances  traveled  at  roughly  70%.  At  operating  point  N,  there  were  12 
spikes  traveling  about  the  annulus  at  70%.  For  any  further  decrease  in  flow,  one  of  these  would 
grow  into  rotating  stall  in  the  same  fashion  as  the  baseline  case.  As  in  the  baseline  case,  the  spike 
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disturbance  gives  an  initial  tip  axial  velocity  deficit  and  a  hub  velocity  increase  until  the  disturbance 
evolves  to  cover  the  full  span.  The  spike  disturbance  circuntferential  size  at  the  hub  is  on  the  order 
of  one  or  two  pitches. 

The  disturbances  with  the  restaggeied  build  have  similar  features  to  the  spike  transients  for 
the  baseline  build.  Both  appear  first  at  operating  point  J  and  travel  at  a  speed  of  roughly  70%. 

Both  have  a  radial  velocity  distribution  with  a  tip  deficit  and  a  hub  increment  in  flow.  The 
circumferential  size  (few  pitches)  of  the  disturbances  of  the  restaggered  build  is  comparable  to  that 
with  the  baseline  build,  early  in  its  growth.  The  stabilized  disturbance  with  the  restaggered  build 
thus  resembles  the  spikes  of  the  baseline  build  at  scmie  time  during  the  five  rotor  revolution 
growth  period  of  the  latter.  It  therefore  appears  that  the  stabilized  disturbances  of  the  restaggered 
build  and  the  transient  disturbance  of  the  baseline  build  may  be  the  same  phenomena,  although  this 
has  not  been  proven  ccmclusively.  This  is  important  because  it  suggests  that  one  can  use 
restaggering  to  enable  more  accurate  flow  diagnosis  of  these  events. 

6.0  Response  Qf  Multistage  Compressors  to  Rotating  Inlet  Flow  Distortions 

The  second  main  thrust  in  this  task  has  been  the  assessment  of  the  effect  of  unsteady  inlet 
disturbances  on  multistage  COTipressor  performance  and  stability.  The  motivation  for  the  work 
arises  out  of  two  compelling  facmrs  concerning  operalnlity  of  gas  turbine  engines.  First  is  that 
stability  limitations  often  arise  not  with  uniform  flow,  but  as  a  result  of  an  inlet  ncMi-uniformity,  or 
distOTticHi.  Further,  distortions  encountered  in  practice  usually  have  significant  asymmetry  to 
ttem,  so  that  the  problon  of  compressor  response  to  circumferential  inlet  distortion  is  a  teal  one 
indeed. 

The  seccmd  aspect  is  that  the  instability  process  is  an  inherently  unsteady  one,  and  the 
behavior  of  multistage  machines,  considered  as  dynamic  systems,  is  not  well  understood.  The 
ctHisequence  is  that  it  is  not  scientifically  correct  to  break  up  the  response  into  smaller,  much 
simpler  modules,  and  study  each  one  separately  in  the  hope  that  one  can  smnehow  learn  about  the 
overall  behavior.  To  do  so  is  to  ignore  the  real  issues  of  the  dynamic  coupling  between  blade 
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passages,  blade  rows,  and  compressor  stages  which  give  rise  to  the  response  of  interest. 

These  background  facts  led  us  to  the  set  of  experiments  to  be  described,  which  made  use  of 
a  rotating  screen  to  create  an  unsteady  inlet  non-uniformity,  which  served  as  a  unique  diagnostic 
tool.  Although  the  use  of  the  rotating  screen  arose  primarily  for  this  reason,  it  can  be  mentioned 
that  circumferential  non-uniformities  which  can  rotate  at  fixed  speed  around  the  compressor 
annulus  can  occur  in  two-spool  compression  systems,  where  one  of  the  compressors  or  the  fan 
could  be  in  rotating  stall,  subjecting  the  other  to  a  non-uniform  rotating  flow  pattern. 

There  is  little  detailed  data  concerning  the  behavior  of  compressors  when  subjected  to  such 
types  of  distortions,  but  experimental  results  [2],  [3]  suggest  that  disturbances  which  rotate  at 
speeds  similar  to  stall  cell  propagation  speed  have  the  most  pronounced  effect.  This  behavior  can 
also  be  seen  in  the  computations  of  [4]  which  linked  the  phenomena  to  the  resonant  response  of 
the  compressor,  i.e.  to  a  correspondence  between  the  forcing  and  the  natural  eigennnodes  of  the 
compression  system. 

The  investigation  of  the  effect  of  rotating  flow  non-uniformities  on  multistage  compressor 
performance  described  here  had  three  central  objectives:  to  obtain  information  on  the  overall 
behavior  (surge  line  changes)  as  a  functitm  of  disttmion  rotation  speed;  to  characterize  the  unsteady 
response  of  the  compressor,  and  to  assess  present  thecnretical  approaches  to  congnessor  stability. 
The  last  is  particularly  important  since  the  existing  models  make  use  of  rudimentary  descriptions 
of  the  unsteady  flow  in  the  compressor  blades,  and  it  is  not  known  to  what  extent  this  affects  the 
predictions  of  ovoall  behavior. 

For  the  present  experiments,  the  General  Electric  large  scale  research  conpressor  described 
earlier  was  equiiq)ed  with  a  drive  mechanism  which  could  rotate  a  distortitm  screen  at  speeds 
between  ±100%  rotor  speed  [5].  The  advantage  of  a  screen-generated  rotating  distortiem,  rather 
than  one  generated  by  another  compressor  in  rotating  stall,  is  that  it  allows  independent  control  of 
the  severity  of  the  non-uniformity,  the  rotation  direction  and  speed  and  the  compressor  operating 
point.  The  screen  was  located  1 .5  radii  upstream  of  the  inlet  guide  vanes  and  produced  a  roughly 
square  wave  stagnation  pressure  distortion  of  120°  circumferential  extent  with  amplitude 
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corresponding  to  1.2  x  [Vi  pCCxmean^^J-  non-uniformity  was  chosen  to  give  a  large  change  in 
the  compressor  stall  point  with  no  rotation,  based  upon  predictions  of  the  theoretical  model. 

Time  mean  and  fast  response  instrumentation  were  used  to  measure  the  mean  overall  static 
pressure  rise  of  the  compressor,  the  unsteady  static  pressure  distribution  (both  upstream  and 
downstream  of  the  compressor)  and  the  time  varying  velocity  field  at  various  locations.  These 
measurements  were  processed  to  analyze  both  the  steady  flow  redistribution,  as  seen  in  the  fiame 
of  reference  of  the  distortion  screen,  and  the  behavior  of  perturbations  which  were  unsteady  in  this 
coordinate  system. 

In  what  follows,  we  shall  present  the  experimental  results  in  close  conjunction  with  theory; 
this  will  serve  to  address  the  capabilities  of  the  flow  models  that  have  been  developed. 

7.0  Uisfio 

The  theoretical  model  used  for  the  predictions  and  the  design  of  the  experiment  was  an 
extension  of  that  reported  in  [4],  [6],  and  [7].  It  is  a  general  fluid  dynamic  stability  analysis  for  the 
compressor  and  the  associated  flowfields.  The  con:q>ressor  performance  is  specified  by  an 
axisymmetric  pressure  rise  characteristic  with  a  correction  for  non-steady  effects  based  on  a  simple 
model  ftn*  the  ineitia  of  the  fluid  within  the  blade  passages.  The  theory  is  well  documented  in  [4], 
[6],  and  [7],  so  only  the  changes  necessary  for  this  investigation  are  commented  upon. 

For  the  present  situadtm,  the  mean  (background)  flow  equations  were  solved  in  a 
coordinate  system  rotating  with  the  inlet  distortion,  so  the  flowfield  redistribution  may  be 
ctmsidered  as  time  invariant  In  the  engiite  frame,  however,  tlte  potential  field  which  causes  the 
upstream  redistributicm  is  unsteady,  and  the  stagnation  pressure  of  a  fluid  particle  is  not  conserved 
even  for  an  inviscid  flow.  This  introduces  an  additional  term  into  the  equations,  but  solution  of  the 
mean  flow  can  be  carried  out  in  the  same  manner.  Once  the  flow  redistribution  has  been  calculated 
(in  the  rotating  coordinates),  the  stability  of  the  flowfleld  is  assessed,  as  in  the  (viginal  form  of  the 
theory,  by  calculatitm  of  the  growth  rates  associated  with  general  unsteady  perturbations. 
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8.0  Experimental  Results  and  Comparison  With  Theory 

The  behavior  of  the  small-amplitude  perturbations,  which  is  important  in  determining  the 
stability  of  the  non-uniform  compressor  flowfield,  is  set  by  the  background  flow,  which  is  a 
measure  of  the  compressor's  response  to  the  imposed  inlet  distortion.  The  theory  shows  that  the 
velocity  field  near  stall  is  primarily  determined  by  unsteady  flow  effects  and  that  these  change 
significantly  as  the  rotation  speed  of  the  distortion  varies.  A  study  of  how  screen  rotation  rate 
affects  the  stability  boundary  of  the  compression  system  is  thus  useful  in  assessing  the  importance 
of  the  unsteady  fluid  dynamics  and  the  degree  to  which  the  modelling  is  correct. 

Figures  17  through  20  show  overall  effects  of  propagation  speed  of  the  inlet  distortion 
(rotation  rate  of  the  screen)  on  the  onset  of  flow  instability  in  the  multistage  compressor.  Tlie 
horizontal  axis  in  the  figure  is  the  non-dimensional  rotation  rate,  rotation  speed  divided  by  rotor 
rotation,  and  the  vertical  axis  is  the  value  of  flow  coefficient,  Cx/U,  at  which  the  instability,  i.e. 
"stall",  is  encountered.  A  higher  value  of  this  quantity  means  a  decreased  stable  flow  range;  the 
compressor  stalls  at  a  higher  flow.  In  the  figures,  the  experimental  data  are  shown  by  the  solid 
lines,  and  the  calculations  by  the  symbols.  There  is  mOTe  than  one  line  in  some  of  the  figures 
because  there  was  a  limitation  on  the  mechanical  speed  of  the  screen.  Data  were  thus  taken  at 
several  rotor  RPM  in  order  to  be  able  to  achieve  a  higher  non-dimensional  value  of  this  quantity. 

Figure  17a  shows  the  compressor  tested  by  Plumley,  as  described  in  [8],  while  Figure  17b 
presents  this  information  in  more  operational  terms.  In  either  of  the  figures,  it  is  seen  that  counter- 
rotation  has  a  markedly  less  severe  effect  on  compressor  instability.  What  is  of  more  interest  is 
the  effect  on  stall  margin,  for  example  as  measured  by  the  distance  between  design  flow  and  stall 
flow.  Figure  17b  shows  that  a  steady  inlet  distortion  decreases  this  margin  by  about  twenty-five 
percent  fiom  the  uniform  flow  value.  With  the  rotating  distortion,  however,  this  decrease  is  more 
than  doubled,  and  the  margin  decreases  to  roughly  one-third  of  the  uniform  flow  value. 

A  scnnewhat  over-simplified  explanation  for  the  decrease  is  that,  for  the  rotation  rates  at 
which  the  large  effects  are  seen,  the  system  is  being  subjected  to  disturbances  that  have  a  temporal 
and  spatial  structure  similar  to  its  natural  eigenmodes.  Because  of  this  "resonance",  in  the  mean 


131 


flow  as  well  as  in  the  unsteady  perturbations,  the  system  (the  compressor)  thus  responds  much 
more  strongly. 

The  theoretical  stability  boundaries  predicted  by  the  model  are  also  shown  for  the  highest 
compressor  speed.  Trends  in  the  predictions  can  be  said  to  be  qualitatively  in  agreement  with  the 
data  but  they  do  not  indicate  the  second  peak  with  frequency,  nw  do  they  show  the  difference 
between  counter-rotational  and  co-rcxational  distortions.  Although  a  definitive  explanation  for  the 
differences  between  these  aspects  of  the  theory  and  experiment  has  not  yet  been  made,  there  is  an 
aspect  of  the  model  behavior  that  may  be  important  and  w.  be  discussed. 

The  calculation  of  the  theoretical  stability  boundary  involved  a  general  flowfield 
disturbance,  with  no  limitations  were  placed  on  characteristic  length  scale  nor  reduced  frequency. 
However,  the  waveform  predicted  for  the  most  unstable  disturbance  is  strongly  dependent  on  the 
screen  rotation  speed.  For  normalized  rotation  rates  of  zero  or  less  (counter-rotation),  the  predicted 
instabilities  have  long  characteristic  length  scales  and  propagate  around  the  annulus  at  20-30%  of 
rotor  speed  in  the  absolute  frame.  In  the  region  where  the  predicted  change  in  the  staWlity  limit  is 
greatest,  calculated  instability  waves  have  a  much  shorter  characteristic  length  scale,  and  do  not 
propagate  relative  to  the  screen  position.  These  localized  disturbances  have  a  high  reduced 
frequency,  may  not  be  adequately  resolved  by  the  simple  model  for  unsteady  flow  in  the  blade 
passages,  and  hence  may  be  more  stable  than  the  model  suggests.  If  so,  this  would  result  in  a 
larger  predicted  change  in  the  stability  boundary  than  that  measured,  as  is  the  case. 

Figure  18  shows  the  data  for  the  compressor  discussed  previously.  It  can  be  seen  that  the 
variation  is  considerably  less  than  with  the  machine  of  Fig.  17,  but  the  two-peaked  structure  and 
the  decrease  in  sensitivity  with  counter-rotation  are  still  present.  (The  confutations  have  not  yet 
been  done  for  this  compressor  build.) 

Figures  19  and  20  give  different  results  for  two  other  compressors.  Again,  the  overall 
picture  is  one  of  a  significant  decrease  in  stable  flow  range  over  a  band  of  rotation  rates,  but  the 
plots  show  a  qualitatively  different  type  of  behavior  than  Figs,  17  and  18.  There  is  a  single  hump 
in  the  overall  plot  of  stall  point  versus  screen  rotating  rate,  in  good  agreement  with  the  theoretical 


description.  In  addition,  although  not  shown  in  the  figures,  the  time-resolved  data  (with  uniform 
flow)  also  show  qualitatively  different  behavior  in  these  last  two  compressors.  Lx)ng  wavelength 
disturbances  can  be  seen  to  propagate  round  the  circumference  prior  to  stall. 

To  summarize  the  overall  experimental  results  as  well  as  the  interaction  between  theory  and 
experiment,  there  appears  to  be  a  direct  link  between  the  overall  behavior  as  regards  the  distcHtion 
response,  and  the  details  of  the  unsteady  stall  inception  process.  Current  mathematical  models  of 
the  flow  capture  the  behavior  very  well  in  cases  in  which  the  long  wavelength  disturbances  exist; 
for  situations  in  which  the  shorter  length  scale  disturbances  are  dominant,  the  models  do  show 
some  features  of  the  flow,  but  the  precise  degree  of  applicability  is  not  clear. 

9.0  Summary  and  Conclusions 

New  fluid  dynamic  phenomena  associated  with  the  stall  inception  process  and  the  dynamic 
response  of  multistage  compressors  have  been  identified  in  the  present  experimental  investigation. 
TThe  results  are  thought-provoking  in  that  they  impact  directly  on  the  development  of  improved 
theoretical  predictive  tools  for  self-excited  as  well  as  forced  unsteady  motions  in  compressors. 
More  specifically,  we  cite  the  following  points  tiiat  have  been  made: 

1)  There  arc  two  routes  to  ctwnpressor  instability,  one  via  the  growth  of  the  amplitude  of  long 
wavelength  disturbances  and  one  in  which  a  shcnt  length  scale  perturbaticm  grows 
drcumfetentially,  as  well  as  perhaps  in  amplitude. 

2)  An  experiment  has  been  carried  out  to  examine  separately  the  role  that  each  of  these 
disturbances  play  in  multistage  machines,  artificially  stabilizing  the  long  wavelength 
disturbances. 

3)  Undo*  conditions  in  which  the  long  wavelengtii  distuitances  were  stable,  die  short  wavelength 
perturtiations  would  appear  and  propagate,  but  would  not  lead  to  rotating  stall. 

4)  The  interaction  of  the  two  types  of  perturbaticms  thus  appears  to  be  cracial  in  die  stall  cmset 
process. 

5)  There  is  a  strong  effect  of  distortion  rotation  speed  and  direction  on  compressor  stability. 
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6)  A  simple  theoretical  model  of  unsteady  compressor  response  to  rotating  distonions  shows 
predictions  that  are  qualitatively  similar  to  the  experimental  results,  for  several  compressors 
and  quantitatively  similar  for  others,  the  difference  between  the  two  classes  appears  to  be  the 
existence  of  small  amplitude,  long  wavelength  disturbances  prior  to  instability. 

7)  Both  the  detailed  rotoi  passage  results  and  the  information  with  the  rotating  inlet  distortion 
point  to  study  of  the  fluid  dynamics  of  the  short  wavelength  disturbances,  both  structure  and 
evolution,  as  an  item  of  extremely  high  priority  for  unsteady  flows  in  multistage  compressors. 
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Fig.  12:  Axial  velocity  at  rotor  exit  (pitchwise  average). 
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Fig.  14:  Pressure  at  rotor  inlet  as  the  throttle  is  closed  (8  equally  spaced  casing  transducers). 
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Fig.  17b:  Calculated  and  experimental  neutral  stability  flow  coefficients 
for  standard  distonion 
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Fig,  19:  Effect  of  distortion  rotation  rate  on  multistage  compressor  stability  (Compressor  3). 
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of  distortion  rotation  rate  on  multistage  compressor  stability  (Compressor  4). 
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TASK  IV:  VORTEX  WAKE-COMPRESSOR  BLADE  INTERACTION  IN 
CASCADES:  A  NEW  RAPID  METHOD  FOR  UNSTEADY 
SEPARATION  AND  VORTICITY  FLUX  CALCULATIONS 

(Investigators:  J.E.  McCune,  A.  Gioulekas) 

1.0  Introduction 

In  previous  progress  reports  under  this  grant,  we  have  described  the  motivation,  as  applied 
to  bluff  airfoils  and  cascaded  blades,  for  the  study  of  unsteady  boundary  layer  separation  in  the 
aerodynamics  of  turbomachineiy.  We  have  also  reported  the  gradual  evolution  under  the  program 
of  a  useful,  much-simplified  method  for  predicting  die  location  and  motion  of  the  boundary  layer 
separation  points  on  such  blades.  This  is  the  key  information  required  to  correctly  predict  the 
interaction  of  an  airfoil  or  a  blade  with  its  own  unsteady  vortical  wake  at  reduced  frequencies  found 
to  be  typical  of  unsteadiness  in  compressor  blade  rows,  and  the  consequent  prediction  of  fluid 
dynamic  forces. 

This  work  has  now  been  completed  and  repotted  cominehensively  in  the  Ph.D.  Thesis  of 
Dr.  Alexandros  Gioulekas,  submitted  to  die  faculQr  of  MIT’s  Department  of  Aeronautics  and 
Astronautics,  and  accepted  in  May  1992.  The  tide  of  Gioulekas’  thesis  suggests  its  application(s) 
in  the  framewoik  of  modem  aerodynamics,  bodi  internal  and  external:  “An  Alternative  to  the 
Kutta  Condition  for  high  Frequency,  Separated  Flows”  (May  1992).  A  copy  of  the  thesis  abstract 
is  atuiched  as  an  appendix,  and  the  complete  thesis  is  available  upon  request. 

A  reproduction  (see  next  section)  of  the  Table  of  (Contents  of  the  thesis  furtha*  provides  a 
useful  overview  of  die  topics  covered  in  this  Task  over  die  past  few  years.  It  also  provides 
(indirectly)  a  listing,  through  Chapter  6  approximately,  of  the  topics  addressed  in  summary  form 
in  the  earlier  progress  reports  to  the  AFOSR  under  this  grant 

In  this  final  rqioil,  in  the  sections  following  the  Gioulekas  Table  of  Contents,  we  focus  on 
our  new  and  more  complete  results,  obtained  since  the  last  progress  repart  The  focus  during  that 
discussion  is  on  the  trends  discovered,  as  a  result  of  this  research,  for  both  the  mean  part  (time 
average)  of  the  blade  response  and  for  the  variations  with  reduced  frequency  of  the  unsteady  part 
Some  of  these  trends  may  appear  surprising  at  first  but  have  been  inter|Hcted  successfully  by 
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Gioulekas  in  his  thesis.  In  addition,  some  very  interesting  “rules  of  thumb”  have  emerged  for  the 
steady  (time  average)  circulation,  or  turning,  that  a  loaded  blade  can  create  in  an  unsteady  flow 
environment. 

2.0  OverYiew  of  Completed  Research 

A  brief  summary  of  the  research  accomplished  under  this  Task  is  effectively  provided  by 
the  Table  of  Contents  of  Gioulekas’  thesis; 

Table  of  Contents* 

1  Introduction  14 

1 . 1  Survey  of  previous  work  and  connection  to  the  present  work  14 

1.2  Synopsis  of  the  thesis  29 

1.3  Overview  32 

2  Flow  in  a  boundary  layer  with  a  rapidly  oscillating  free-stream  velocity  33 

2.1  Assumptions  33 

2.2  The  division  of  the  flow-field  34 

2.3  The  non-dimensional  form  of  the  equations  41 

2.4  The  solution  48 

2.5  !Zero-th  order  approximation  49 

2.6  Second  order  approximation  51 

3  Unsteady  separation  59 

3. 1  Conditions  for  unsteady  separation  59 

4  A  criterion  for  predicting  unsteady  separation  71 

4.1  Derivation  of  the  criterion  71 

4.2  Nondimensional  form  of  the  unsteady  sq>aration  criterion  79 

5  How  the  interaction  betweoi  tiie  airfoil  and  its  wake  detomines  the  airfoil 

circulation  and  tile  force  and  moment  acting  (HI  the  airfoil  83 

5.1  The  circulation  84 

5.2  Airfoil  in  oscillating  stream  vs.  oscillating  airfoil  in  steady  stream:  udiat 

is  the  diffnence  in  the  aerodynamic  force  and  moment?  91 

5.3  Inertial  and  airfoil  fiames  of  reference  92 

5.4  Calculation  of  the  force  using  the  unsteady  B^ouUi  equation  93 

5.4.1  The  pressure  coefficient  93 

5.4.2  The  mapping  of  the  physical  to  the  circle  plane  95 

5.4.3  The  force  and  moment  coefficient  96 

5.5  Calculation  of  the  force  using  the  impulse  96 

5.6  Calculation  of  the  moment  using  the  moment  of  impulse  1(X) 

5.7  The  fiee  wake  convection  104 


*  Extracted  from  Gioulekas’  Ph.D.  Thesis.  The  page  numbers  refer  to  tiic  thesis,  not  to  this  report 
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6  The  influence  of  reduced  firequency,  strength  of  flow  unsteadiness,  and 

ellipse  slenderness  on  unsteady  separation  1 06 

6. 1  Comparison  between  th^retical  predictions  and  experiment  106 
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The  results  discussed  in  the  following  sections  are  selected  primarily  from  Chapters  6  and 
7  of  the  above  thesis. 


3.0  Trajectories  of  the  Unner-Siirface  Sena 

wi  Bluff  Airfoiis  or  Blades  at  Angle  of  Attack 
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Application  of  the  generally  accepted  separation  criteria  for  unsteady  boundary  layor  flow 
for  practical  problems  involving  blade-wake  interaction  has  long  been  recognized  as  quite 
coirplex.  Direct  numerical  study  has  been  awkward,  time  consuming  and  expensive,  with  results 
diat  are  both  difBcult  to  interpret  and  to  generalize. 

But  in  die  research  reported  in  Gioulekas'  diesis,  and  in  our  earlier  progress  reports  to 
AFOSR,  we  have  been  able  to  achieve  a  substantial  simplifleation  of  this  conqiutational  task 
without  sacrificing  the  important  physics  involved.  This  is  pardculariy  so  when  the  unsteadiness 
occurs  at  moderately  high  ‘^reduced  frequencies”,  taefU  >  1,  for  reasons  originally  pointed  out 
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separately  by  C.C.  Lin  and  by  Lighthill.  In  addition  to  the  simplifications  suggested  by  Lin,  we 
have  been  able  to  add  an  extension  of  the  ideas  of  Stratford*  to  provide  an  algebraically  simple 
computational  criterion  for  separation,  involving  only  external  par^eters  of  the  flow,  steady  or 
unsteady. 


The  result  is  that  we  can  now  readily  calculate  the  instantaneous  location  and  motion  of  the 
boundary  layer  sqiaration  point(s),  without  having  to  go  through  a  complete  boundary  layer 
calculation  at  each  time  step.  This  means  that  when  a  blade  or  an  airfoil  is  interacting  with  its  own 
unste^y  wake,  the  effect  of  the  wake  itself,  back  on  the  boundary  layer  that  is  “feeding”  it,  can  be 
readily  included,  and  this  is  a  key  feature  of  the  results  reported  here. 

The  examples  included  below  illustrate  the  kind  of  separation  trajectories  one  finds  above  a 
slender,  bluff  airfoil  in  unsteady  flow,  and  how  these  depend  on  some  of  the  key  parameters 
involved.  The  first  cases  shown  are  for  an  ellipse  of  slenderness  ratio  1  ;20  at  zero  mean  angle-of- 
attack.  In  that  case,  of  course,  the  separation  trajectories  on  the  lower  surface  of  the  blade  are  the 
mirror  image  of  those  shown  here.  Results  for  finite  mean  angle-of-attack  are  discussed  in  the 
following  section,  in  terms  of  the  blades’  aerodynamic  response,  both  for  the  steady  and  for  the 
unsteady  component 

.  The  value  of  the  reduced  frequency,  ox^oo,  in  the  first  case  shown  for  this  section,  is 


^2-  (OC  _n 

A,  -y  -y  , 


and  results  are  shown  as  a  function  of  the  strength  of  the  unsteady  oscillations,  as  measured  by 


H«e,  c  is  the  airfoil  chord  and  the  external  flow  is  U„, = U,,  +  0«,  e'^“* 

Figure  1,  extracted  from  Chapter  6  of  the  thesis,  shows  separation  point  trajectories  near 
the  ellipse’s  trailing  edge  fOT  e  =  0.01, 0.OOS,  0.001,  and  0.  The  steady  flow  separation  point,  as 


*  An  extensive  bibliogn^hy  is  provided  in  Gioulekas’  thesis,  aivi  has  also  been  made  available  in 
earlier  AFOSR  progress  reports.  A  feature  of  Gioulekas’  thesis,  in  fact,  is  the  usefulness  of  his 
bibliogrtqrhy. 
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calculated  by  Howarth  (sec  discussion  in  previous  reports,  or  in  the  thesis),  is  also  shown  for 
reference.  Each  trajectory  shown  represents  the  path  traversed  by  the  separation  point  over  one 
period  of  oscillation. 

The  Hgure  shows  that  as  the  strength  or  amplitude  of  unsteadiness  decreases  from  1%  to  0, 
the  trajectories  shrink  in  range  and  converge  on  the  steady-flow  separation  point,  as  would  be 
expected.  However,  one  also  readily  sees  that  not  all  effects  are  simpiy  “linear”  in  e. 

Figure  2  gives  corresponding  results  for  fixed  e,  but  with  varying  reduced  frequencies, 
ranging  from  9  to  36.  At  first,  these  results  may  seem  surprising,  because  one  sees  that  as  the 
frequency  increases  at  fixed  amplitude,  the  trajectories  again  converge  toward  the  steady-flow 
separation  location. 

This  phenomenon  is  explained  by  Gioulekas  in  his  thesis.  In  simple  terms,  at  high 
frequency  the  response  of  the  unsteady  part  of  the  flow  to  the  wall's  resistance  is  confined  to  an 
ever-thinner  sheath,  analogous  to  Stokes  flow  for  an  oscillating  infinite  wall.  In  the  limit,  this  inner 
sheath  becomes  decoupled  from  the  “Prandtl-typc”  boundary  layer,  which  is  the  part  of  the  flow 
that  responds  primarily  to  the  steady  part  of  the  wall's  resistance,  that  is,  to  the  mean  (time- 
averaged)  flow.  Thus,  at  veiy  high  reduced  frequencies,  the  “Stokes  layer”  and  the  “Prandtl  layer” 
become  decoupled,  and  in  the  limit,  the  boundary  layer  separation  properties  are  determined 
completely  by  the  mean  flow,  precisely  because  the  frequency  is  high. 

4.0  Acrodvnimiic  Response  for  Loaded  Separated  Blades 

With  the  new  method  made  available  by  his  analysis  of  the  problem,  and  by  his 
corresponding  extension  of  the  “Stratford  Criterion”  to  unsteady  flow,  Gioulekas  was  able  to 
systematically  investigate  the  aerodynamic  response  of  bluff-edged  airfoils  over  a  range  of  blade 
slenderness  and  of  mean  angle-of-attack. 

As  we  have  pointed  out  in  our  earlier  reports,  die  new  mediod  of  rapid  location  of  the 
separation  trajectories  plays  the  same  role  in  determining  the  blade  cm.  ulation  as  does  the  classical 
Kutta  condition  for  unseparated,  sharp-edged  blades  or  airfoils.  It  is  natural,  therefore,  to  compare 
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the  results  Gioulekas  obtained  with  those  that  would  be  olrtained  if  the  Kutta  condition  were 
assumed  to  apply. 

A  second  benchmark  is  the  circulation  that  would  be  obtained  in  steady  flow  over  the 
separated  bluff  airfoil  at  angle-of-attack,  as  first  obtained  by  Howarth. 

With  these  limits  in  mind.  Tables  1  and  2  provide  useful  and  revealing  summary 
information  to  guide  our  understanding  of  blade  performance  in  unsteady  flow. 

In  Table  1,  the  results  for  an  elliptic  airfoil  of  fineness  1:20,  at  angle-of-attack  =  5**,  and  for 
an  amplitude  of  unsteadiness  e  =  0.04,  are  shown  for  a  range  of  reduced  frequencies.  Tabulated 
are  the  mean  (time-averaged)  circulation  T  developed  by  the  blade,  the  time-varying  circulation  F, 
and  the  phase-shift,  to  be  compared  with  classical  unsteady  airfoil  theory.  The  values  on  the  left 
are  tiiose  that  are  obtained  if  the  Kutta  condition  is  (arbitrarily)  assumed.  Those  on  the  right  are  the 
responses  obtained  using  the  present  theory,  and  are  thus  based  on  the  more  realistic  prediction  of 
the  instantaneous  location  and  movement  of  the  boundary  layer  separation  point(s). 

We  note  that  the  values  of  T  actually  achieved  by  the  blade  are  systematically  Iowct  than 
tiiose  it  “could”  achieve  if  Kutta  a4)plied. 

We  also  note  that  FHowaitht  the  steady  flow  value  of  circulation  obtained  by  the  same 
airfoil,  at  the  same  angle-of-attack,  taking  into  account  the  separation  of  the  upper-surface 
boundary  layer,  would  be 

^Howarth  =  0.267  , 

about  8%  less  than  FKutta-  Based  on  our  observations  of  the  separation  trajectories  described  in 
Secti(Hi  3.0,  it  is  clear  why,  as  the  fiequency  increases,  the  actual  F  approaches  the  lower, 
separated  steady-flow  value,  Fifowaith^  for  a  bluff-edged  airfoil. 

In  any  case,  the  results  provide  us  with  useful  bounds  for  tire  mean  circulation  obtainable 
by  a  bluff-edged  blade  at  angle-of-attack  in  any  unready  environment: 

^Howarth  ^  ^  ^KutUt  • 

These  bounds  can  be  used,  for  example,  to  estimate  the  maximum  loss  of  mean  turning  angle  that 
might  occur  in  a  blade  row  as  a  result  of  unsteadiness.  The  higher  the  frequency,  the  closer  the 
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lower  value  is  approached.  And,  at  fixed  ftequency,  the  thinner  the  blade,  the  more  closely  PRutta 
can  be  obtained. 

Table  2  serves  as  one  check  on  the  theory,  and  verifies  that,  as  the  blade  is  chosen  thinner 
and  thinner,  conditions  approach  those  appropriate  to  the  classical  Kutta  condition,  other 
parameters  being  held  constant  This  is  true  not  only  for  the  mean  circulation.  F,  but  also  for  the 
properties  of  the  unsteady  component  of  the  blade  response.  The  amplitudes  F,  and  the  phase 
angles  0r,  obtained  from  the  new  theory,  agree  remarkably  well ,  for  very  thin  elliptical  blades, 
with  the  results  of  classical  thin  airfoil  theory.  Table  2  shows  this  agreement  for  the  example  of  an 
elliptical  blade  of  fineness  ratio  1 :50,  at  reduced  frequency  =  9. 

Summary  conclusions  of  the  kind  given  in  the  above  two  tables  were  obtained  by 
systematic  study  of  the  results  obtained  for  many  cases  (each  for  a  given  set  of  the  above- 
mentioned  parameters)  by  runs  such  as  are  illustrated  in  Figs.  3  dirough  5.  These  figures  show 
how  the  aerodynamic  responses  were  actually  obtained  and  analyzed,  and  illustrate  only  a  small 
sample  of  the  results  reported  by  Gioulekas. 

5.0  Future  and  Ongoing  AppMcations 

The  extensive  results  available  in  Gioulekas’  thesis  rq>resent  a  substantial  study  of  the 
blade-wake  interaction  effect  for  driven,  or  imposed,  unsteady  oscillations,  including  realistic 
treatment  of  the  boundary  layer  separation  properties. 

Future  work  includes  qrplication  of  this  technique  to  the  determination  of  nahiral- 
osdllation  regimes  (analogous  to  die  Kaiman  vortex  street  behind  a  cylinder)  for  thin  elliptical 
blades  or  airfoils. 

The  method  also  has  direct  application  to  the  dynamic  stall  problem  for  airfoils  in  external 


flow. 
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(a)  Kutta  Condition 

(b)  Separation 

A* 

r 

f 

r 

r 

St 

i 

0.288 

0.0015 

-50® 

0.279 

0.003 

-46® 

1 

0.288 

0.001 

-52® 

0.272 

0.002 

-59® 

m 

0.288 

0.0005 

-54® 

0.271 

O.OOl 

-65® 

TABLE  1:  The  circulation  r(t)  =  F  +  f  cos(o)t  +  Gp)  around  an  ellipse  calculated  by:  (a) 

assuming  smooth  flow  at  the  trailing  edge,  and  (b)  using  t^e  new  separation  criterion, 
^  the  reduced  frequency  increases  (b :  a  =  1 : 20,  a  =  5°,  e  =  0.04).  Note  that 
r  ->  FHowarth  =  0.267  as  the  reduced  frequency  increases. 


( 


Slenderness 

(a)  Kutta  Condition 

(b)  Separation 

Ratio 

r 

A 

r 

St 

r 

f 

1:20 

0.288 

0.0015 

-50® 

0.279 

OJ0O3 

-46® 

1:50 

0.279 

0.002 

-46® 

0.279 

0.002 

-46® 

TABLE  2:  The  circulation  F(t)  =  F  +  f  cos(oi)t  +  Op)  around  an  ellipse  calculated  by:  (a) 

assuming  smooth  flow  at  the  trailing  edge,  and  (b)  using  die  new  separation  criterion, 
as  the  ellipse  becomes  slimmer  (a  =  5®,  =  9,  e  =  0.04). 


( 


*  All  “F  values”  presented  here  are  made  dimensionless  with  U«oC,  i.c.,  each  number  is 
(cjrculation)AJoaC. 
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,  Fig.  1 :  The  separation  trajectory  for  varying  strength  of  flow  unsteadiness,  e  =  0.00 1 , 0.005, 

1  0.01,  and  constant  reduced  frequency,  =  9.  The  separation  trajectory  shrinks  and 

approaches  the  location  of  steady  separation  as  e  approaches  zero. 

I 

I 
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Fig.  3:  The  stationary  points  t  =  T/4,  /;  t  =  3T/4  of  the  suction-side  separation  trajectory  roughly 
correspond  to  the  stationary  points  in  r*(t).  The  ellipse  (b/a  =  1 : 20)  is  placed  at  an 
angle  of  attack  a  =  5°  in  a  stream  oscillating  with  strength  e  =  0.04  at  reduced  frequency 
X2  =  9. 


1 


— Q  0  —  ^ 
-0  0-  > 


Fig.  4:  Time  evolution  of  the  non-dimensional  circulation  for  varying  reduced  frequency  =  9, 
18, 36)  and  constant  strength  of  flow  unsteadiness  (e  =  0.04).  The  circulation  is 
bracketed  by  rKutia  =  0.283  and  FHowarth  =  0.267.  The  ellipse  (b/a  =  1 : 20)  is  at  an 
angle  of  attack  a  =  5®. 
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I 

I 
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( 


Fig.  5:  Time  evolution  of  the  lift  coefficient  for  varying  reduced  frequency  =  9, 1 8, 36)  and 

constant  strength  of  flow  unsteadiness  (e  =  0.04).  The  ellipse  (b/a  =  1 : 20)  is  at  an  angle 
of  attack  a  =  5®. 
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APPENDIX  IV.l 


An  alternative  to  the  Kutta  condition 
for  high  frequency,  separated  flows 

by 

Alexandros  Gioulekas 

Submitted  to  the  Department  of 
Aeronautics  and  Astronautics 
in  partial  fulfillment  of  the 
requirements  for  the  Degree  of 

Doctor  of  Philosophy 

Abstract 


An  alternative  to  the  Kutta  condition  for  determining  the  circulation  around  a  bluff 
airfoil  in  unsteady,  separated  flow  is  presented.  For  such  flows,  there  is  a  need  for  a 
practical  criterion  which  would  avoid  the  detailed  boundary  layer  calculations  and  would 
predict  the  time  evolution  of  the  airfoil  circulation  based  on  the  external  potential  flow 
only. 

This  criterion  would  play  for  unsteady,  separated  flows  the  role  that  the  Kutta 
condition  plays  for  flows  past  thin  airfoils.  It  turns  out  to  be  a  criterion  for  predicting 
the  location  and  movement  of  the  “separation  points”,  because  they  determine  the  net 
vorticity  flux  shed  into  the  wake  and  thus  the  rate  of  change  of  the  airfoil  circulation. 

The  laminar,  two-dimensional  flow  about  a  bluff  airfoil  at  angle  of  attack,  when  the 
external  flow  oscillates  at  a  high  reduced  frequency  is  considered. 

At  high  frequencies,  the  vorticity  generated  as  the  wail  resists  the  imposed  unsteadi* 
ness  is  confined  to  a  thin  layer  near  the  blade  surface  (“Stokes  layer”)  and  its  contri* 
bution  to  the  displacement  thickness  is  proportional  to  an  inverse  power  of  the  reduced 
flrequency  and  thus  small.  Outride  this  region,  the  unsteady  part  of  the  boundary  layer 
velocity  is  approximately  the  external  potential  oscillation.  Based  on  this  observation 
and  following  C.C.  Lin,  the  boundary  layer  velocity  can  be  divided  into  two  coupled 
velocity  distributions,  one  predominantly  oscillatory  (“Stokes  velocity”)  and  another 
predominantly  steady  (“Prandtl  velocity”).  The  main  contributor  to  the  diq>lacement 
thickness  is  the  latter.  Therefore,  separation,  identified  by  a  dramatic  increase  in  the 
displacement  thickness,  can  be  located  by  calculating  the  evolution  of  the  “Prandtl 
velocity” and  finding  where  the  latter  bifurcates. 

.  / 

Stratford’s  ideas,  modified  to  account  for  unsteadiness  in  the  “Prandtl  velocity' 
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(arising  through  the  coupling  of  the  “Prandtl”  to  the  'Stokes”  Sow  by  the  no^slip 
condition  on  the  wall  and  Reynolds-streas  terms  in  the  momentum  equation),  lead  to 
a  criterion  for  unsteady  separation  that  uses  as  only  inputs  parameters  of  the  external 
Sow  and  avoids  a  detailed  boundary  layer  calculation. 

The  airfoil  circulation  is  calculated  by  an  iterative  method  which  calculates  how  the 
interaction  between  the  airfoil  and  its  wake  affects  separation. 

An  ellipse  is  adopted  as  a  study  case,  and  results  are  presented  for  varying  angle 
of  attack,  ellipse  slenderness,  reduced  frequency,  and  strength  of  unsteadiness.  In  the 
limit  of  a  very  slender  ellipse,  the  theory  recovers  the  results  from  the  classical  unsteady 
wing  theory,  which  assumes  the  Kutta  condition. 

The  theory  predicts  that  there  exist  two  limits  for  the  mean  value  of  the  circulation. 
The  upper  limit  is  the  value  of  the  circulation  for  which  the  trailing  edge  becomes  a 
stagnation  point  (FKutta)-  The  lower  limit  is  the  value  of  the  circulation  in  steady  flow 

(FHowtrth)- 

The  pressure>side  "separation  point”  for  all  practical  purposes  can  be  considered 
fixed,  even  at  small  angles  of  attack.  On  the  other  hand,  the  "separation  point”  on  the 
suction-side  oscillates  with  amplitude  proportional  to  the  strength  of  the  flow  unsteadi¬ 
ness,  and  inversely  proportional  to  the  reduced  frequency.  When  the  reduced  frequency 
increases  or  the  strength  of  flow  unsteadiness  decreases,  the  trajectory  of  this  "separa¬ 
tion  point”  shrinks  and  tends  toward  the  position  of  steady  separation.  Since  the  mean 
I  location  of  the  suction-side  "separation  point”  controls  the  mean  value  of  the  circu¬ 

lation,  and  the  amplitude  of  its  excursion  determines  the  amplitude  of  the  oecillatory 
component  of  the  circulation,  the  above  trends  explain  how  the  circulation  responds  to 
changes  in  the  flow  unsteadiness. 


Thesb  Supervisor:  James  E.  McCune 


Title:  Professor  of  Aeronautics  and  Astronautics 
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3.  AIR  FORCE  RESEARCH  IN  AERO  PROPULSION  TECHNOLOGY 
(AFRAPT)  PROGRAM 


The  research  at  MIT  is  strongly  tied  in  with  the  AFRAPT  Program.  The  students  who 
have  participated  in  this  program  during  the  contract  period,  the  advisors,  and  the  research  projects 


are: 

Trainee: 

Advisor: 

Project: 

Jabin  Bell 

Prof.  A.H.  Epstein 

Active  Control  of  Jet  Engine  Surge 

Trainee: 

Advisor: 

Project: 

Aaron  Gleixner  (fmished  M.S.  Thesis,  8/92) 

Prof.  A.H.  Epstein 

Unsteady  Phenomena  in  Transonic  Axial  Flow  Turbines  (RR) 

Trainee: 

Advisor: 

Project: 

Martin  Graf 

Prof.  A.H.  Epstein 

Characterization  of  Fully  Scaled  Turbine  Blade  Boundary  Layers 

Trainee: 

Advisor; 

Project: 

Daniel  Gysling  (finished  M.S.  Thesis,  8/89;  Ph.D.  expected,  8/93) 

Prof.  J.  Dugundji/Prof.  E.M.  Greitzer/Prof.  A.H.  Epstein/Prof.  J.  Kerrebrock 
Active  and  Structural  Control  of  Aeromechanical  Systems  (AFOSR) 

Trainee: 

Advisor: 

Project: 

Dana  Lindquist  (finished  Ph.D.  Thesis,  12/91) 

Prof.  M.B.  Giles 

Development  of  Adaptive  Procedures  for  Turbomachinery  Flow  Computations 

Trainee: 

Advisor: 

Project: 

Knox  Millsaps  (finished  Ph.D.  Thesis,  5/92) 

Prof.  M.  Martinez-Sanchez 

Rotor  Dynamic  Instability  Due  to  Alford  Forces  (NASA  MSFC) 

Trainee: 

Advisor: 

Project: 

G.  Scott  McNulty 

Prof.  E.M.  Greitzer 

Self-Adaptive  Aeropropubion  Conqx>nents 

Trainee: 

Advisor: 

Project: 

Taras  Palczynsld  (firtished  M.S.  Thesis,  8/92) 

Prof.  M.  Martinez-Sanchez 

Tip  Clearance  Excited  Turbomachine  Rotordynamic  Instabilities  (NASA  MSFC) 

Trainee: 

Advisor: 

Project: 

George  Pappas  (firtished  M.S.  Thesis,  8/90) 

Prof.  A.H.  Epstein 

Temperature  Distortion  Effects  on  Turbine  Heat  Transfer  (RRI) 
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Trainee:  Peter  Silkowski  (finished  M.S.  Thesis,  2/90;  left  MIT,  8/92) 

Advisor:  Prof.  E.M.  Greitzer/Dr.  C.S.  Tan/Prof.  N.A.  Cumpsty  (Cambridge  University) 

Project:  Unsteady  Phenomena  and  Flow  Instabilities  in  Multistage  Turbomachines 

(AFOSR) 


A  more  detailed  description  of  the  work  carried  out  on  these  projects  is  given  in  Refs.  3. 1 

3.2,  and  3.3. 
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AFOSR-85-0288,  September  1990. 
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4.  PUBLICATIONS  AND  PRESENTATIONS  (Period  10/91  - 10/92) 


Epstein,  A.H.,  “Unsteady  Effects  in  Turbine  Aerodynamics  and  Heat  Transfer,”  ASEA  Brown 
Boveri  Research  Laboratory,  Baden,  Switzerland,  October  1992. 

Fritsch,  G.,  Giles,  M.,  “Second-Order  Effects  of  Unsteadiness  on  the  Performance  of 
Turbomachines,”  ASME  Paper  92-GT-389,  presented  at  the  International  Gas  Turbine  and 
Aeroengine  Congress  and  Exposition,  Cologne,  Germany,  1992. 

Giles,  M.B.,  “An  Approach  for  Multi-Stage  Calculations  Incorporating  Unsteadiness,”  ASME 
Paper  92-GT-282,  presented  at  the  International  Gas  Turbine  and  Aeroengine  Congress  and 
Exposition,  Cologne,  Germany,  1992. 

Giles,  M.B.,  “Current  Status  of  CFD  Capability  for  Unsteady  Flow  Analysis,”  presented  at 
NASA  Lewis  Workshop  on  “Unsteady  Row  in  Turbomachines;  Potential  Benefits,  Present 
Status,  Future  Research”,  May  1992. 

Greitzer,  E.M.,  Tan,  C.S.,  “Unsteady  Row  in  Turbomachines:  Where’s  the  Beef?”,  presented  at 
NASA  Lewis  Workshop  on  “Unsteady  Row  in  Turbomachines:  Potential  Benefits,  Present 
Status,  Future  Research”,  May  1992. 

Guenette,  G.R.,  Pappas,  G.,  Epstein,  A.H.,  “The  Influence  of  Non-Uniform  Spanwise  Inlet 
Temperature  on  Turbine  Rotor  Heat  Transfer,”  presented  at  80th  PEP/AGARD  Symposium, 
Antalya,  Turkey,  October  1992. 

Longley,  J.P.,  Greitzer,  E.M.,  “Inlet  Distortion  Effects  in  Aircraft  Prqiulsion  System  Integration,” 
AGARD  Lecture  Series  183  on  “Steady  and  Transient  Performance  Rediction  of  Gas  Turbine 
Engines,”  June  1992. 
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5.  PROGRAM  PERSONNEL 

PrincipallnYfiSiigator; 

Edward  M.  Greitzer 

H.N.  Slater  Professor  of  Aeronautics  and  Astronautics 
Director,  Gas  Turbine  Laboratory 

Co-lnvesti^tors: 

Alan  H.  Epstein 

Professor  of  Aeronautics  and  Astronautics 
Associate  Director,  Gas  Turbine  Laboratory 

Michael  B.  Giles 

Associate  Professor  of  Aeronautics  and  Astronautics 
James  E.  McCune 

Professor  of  Aeronautics  and  Astronautics 
Choon  S.  Tan 

Principal  Research  Engineer 


Graduate  Research  Assistants: 


9/89  -  8/92 
7/88  -  5/92 
6/88  -  5/92 
9/90  -  present 
1/92  -  present 


Peter  Silkowski  (AFRAPT  student) 
Gerd  Fritsch 
Alexandros  Gioulekas 
Martin  Graf  (AFRAPT  student) 
David  Sujudi 
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6.  INTERACTIONS 

As  we  have  noted  previously  and  as  may  be  inferred  from  the  presentations  that  were 
given,  there  are  considerable  interactions  between  Gas  Turbine  Laboratory  personnel  and  industry 
and  government,  as  well  as  with  other  universities.  This  occurs  in  several  forms.  One  of  these  is 
collaborative  research  projects  with  industry  and  government  laboratories.  This  creates  a  direct 
link  with  those  persons  who  are  actively  working  in  similar  areas  in  other  establishments,  and  aids 
considerably  in  helping  transfer  technology. 

A  second  type  of  interaction  is  through  presentations  and  discussions  with  industry  or 
government  agencies.  These  were  listed  in  Section  4  covering  the  contract  period.  It  should  be 
noted  that  only  those  which  involved  discussion  of  the  work  supported  by  AFOSR  under  this 
grant  are  given. 

The  Gas  Turbine  Laboratory  also  has  an  active  seminar  program  which  brings  speakers 
from  industry  and  government.  The  speakers  for  the  past  ''ear  are  listed  below;  those  for  previous 
years  have  been  listed  in  Refs.  [6.1]  and  [6.2]. 

Mr.  W.  Culhane,  American  Airlines 

‘Today’s  Centralized  Fault  Detection  System,  Maximum  Information 

Yielding  Minimum  Value” 

Mr.  B.  Koff,  Pratt  &  Whitney  Aircraft 

“Spanning  the  Globe  with  Jet  Propulsion” 

Mr.  J.  Kuhlberg,  Pratt  &  Whitney  Aircraft 

“Jet  Engine  Controls:  Past,  Present,  and  Future” 

Dr.  D.  McCormick,  United  Technologies  Research  Center 

“Vortical  and  Turbulent  Structure  of  Lobed  Mixer  Free  Shear  Layers” 

Dr.  L.  Nichols,  NASA  Lewis  Research  Center 

“Numerical  Propulsion  System  Simulation” 

As  mentioned  in  the  description  of  Task  ID,  there  is  strong  coupling  with  the  Whittle 
Laboratory  at  Cambridge  University.  This  laboratory  is  of  comparable  size  to  the  Gas  Turbine 
Laboratory  and  has  an  excellent  and  very  active  group  of  faculty  and  students  focussed  on 
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turbomachinery  research.  Another  close  interaction  which  spans  several  of  the  tasks  is  with 
Professor  F.E.  Marble  of  Caltech.  Marble  is  the  former  Head  of  the  Jet  Propulsion  Center  in  the 
Division  of  Engineering  and  Applied  Sciences,  and  has  worked  on  a  wide  range  of  fluids 
problems  connected  to  propulsion.  He  has  participated  in  discussions  in  the  Gas  Turt>ine 
Laboratory  for  over  a  decade,  and  his  physical  insight  and  ability  to  penetrate  to  the  heart  of 
complex  issues  have  been  enormously  helpful  to  us  in  that  time. 

E.  M.  Greitzer  spent  a  month  at  the  Ohio  Aerospace  Institute  as  a  Lecturer  in  Residence  on 
Turbomachinery.  This  included  a  series  of  lectures  on  turbomachinery  aerodynamics  at  NASA 
Lewis  Research  Center,  as  well  as  lectures  and  discussions  with  industry  and  government 
laboratories  and  universities  in  the  Ohio  region:  General  Electric  Aircraft  Engines,  Wright 
Laboratories,  Allison  Gas  Turbine  Division,  University  of  Cincinnati,  and  Case  Western  Reserve 
University. 

E.M.  Greitzer,  C.S.  Tan,  and  M.B.  Giles  participated  in  a  NASA/industry/university 
workshop  on  unsteady  flow  in  turbomachines.  The  purpose  of  the  workshop  was  to  assess 
available  information  on  the  importance  of  unsteady  flow  in  practical  turbomachines  and  see 
whether  a  consensus  existed  on  the  magnitude  of  the  gains  to  be  achieved  through  better 
understanding  of  unsteady  flow,  and  the  issues  that  were  central  to  the  potential  gains.  A  further 
goal  was  to  define  and  prioritize  those  fluid  dynamic  topics  in  unsteady  flow  on  which  attention 
should  be  focused.  The  participants  had  strong  consensus  that  critical  unsteady  flow  problems, 
with  substantial  payoff,  exist  in  the  areas  of  gas  turbine  engine  and  compressor  stability,  and  in 
aeroelasticity  in  turbomachines.  Much  of  the  discussion,  therefore,  was  aimed  at  the  area  of 
greatest  uncertainty,  the  unsteadiness  associated  with  rotor  stator  interaction,  i.e.  unsteady  flow 
effects  which  have  a  length  scale  on  the  order  of  the  blade  pitch. 

A  report  on  the  woikshop  has  been  written  by  J.J.  Adamczyk  (NASA),  A.J.  Strazisar 
(NASA)  and  E.M.  Greitzer  (MIT).  A  summary  of  this  report  is  given  below: 
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Summary 

NASA  Government! Industry lUniversity  Workshop  on 

Unsteady  Flow  In  Turbomachines:  Potential  Benefits,  Present  Status,  Future  Research 

On  May  20-21, 1992  a  workshop  was  held  at  NASA  Lewis  Research  Center  on  the  topic 
of  “Unsteady  How  in  Turbomachines”.  The  purpose  of  the  meeting  was  to  bring  together  active 
workers  in  the  field,  from  university,  industry,  and  government,  to  address  the  following 
questions: 

1.  What  areas  show  potential  for  gas  turbine  engine  performance  increase  through  improved 
understanding  of  unsteady  flow? 

2.  What  are  the  magnitude  of  the  gains  to  be  expected? 

Associated  with  the  answers  to  (1)  and  (2)  are  several  follow-on  issues; 

3.  What  is  the  current  state  of  the  art  in  understanding  in  the  areas  defined  to  be  of  importance? 

4.  What  topics  should  be  explored  further? 

5.  How  does  one  do  this  in  an  appropriate  manner  (What  is  the  balance  of  basic  and  applied 
research,  near  term  versus  long  term,  etc.)? 

6.  If  C.F.D.  codes  are  to  be  used  to  quantify  the  effects  of  unsteady  flows  on  aero-performance 
and  surface  heat  transfer,  how  should  these  codes  be  calibrated? 


A.  H.  Epstein  participated  in  a  study  of  wide  breadth  carried  out  by  the  Aeronautics  and 
Space  Engineering  Board  of  the  National  Research  Council.  The  purpose  was  to  assess  the  current 
state  of  aeronautics  in  the  United  States  and  to  help  identify  the  technology  advances  necessary  to 
meet  the  challenges  of  the  next  several  decades.  The  conclusion  of  the  study  are  summarized  in 
the  report  "Aeronautical  Technologies  for  the  Twenty-First  Century"  issued  by  the  ASEB. 

E.  M.  Greitzer  participated  in  an  Air  Force  Scientific  Advisory  Board  Committee  to 
examine  failures  in  titanium  components  of  Air  Force  aircraft  engines. 
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7.  DISCOVERIES,  INVENTIONS,  AND  SCIENTIFIC  APPLICATIONS 


During  the  present  contract  period,  there  have  been  no  inventions. 
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8.  CONCLUDING  REMARKS 

Several  aspects  of  the  multi-investigator  program  deserve  comment.  First,  one  can  see 
explicitly  in  the  research  conducted  strong  technical  links  with  industry  and  government 
laboratories.  A  consequent  benefit  is  the  connection  with  problem  areas  that  are  of  real  interest; 
this  is  useful  and  stimulating  to  faculty  and  students  alike.  A  second  positive  feature  is  that  there 
are  personnel  at  these  organizations  who  are  directly  involved  in  the  programs.  These  people  are, 
within  the  organizations,  the  “owners”  of  the  scientific  and  engineering  results  that  emerge  from 
the  research.  In  this  sense,  there  is  much  less  need  for  other  formal  mechanisms  for  technology 
transition,  since  this  occurs  at  the  grass  roots  level. 

Another  item  is  the  strong  coupling  between  investigators,  which  has  been  fostered  by  this 
program.  We  believe  this  is  extremely  useful  in  “cutting  edge”  aeropropulsion  research;  the 
problems  are  not  only  conceptually  difficult,  but  arc  also  complex.  Having  several  different  types 
of  talents  focussed  on  the  different  areas  provides  a  substantial  synergism  to  the  overall  program. 
The  multi-investigator  aspect  of  the  research  program  is  a  powerful  resource. 


